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SUMMARY 


The performance of the Space Shuttle Orbiter's Ku-Band Integrated 
Radar and Communications Equipment is analyzed for the radar mode of 
operation. The block diagram of the rendezvous radar subsystem is 
described. Power budgets for passive target detection are calculated, 
based on the estimated values of system losses. Requirements for process- 
ing of radar signals in the search and track modes are examined. Time- 
multiplexed, single-rchannel , angle tracking of passive scintillating 
targets is analyzed. Radar performance in the presence of main lobe 
ground clutter is considered and candidate techniques for clutter sup- 
pression are discussed. Principal system parameter drivers are examined 
for the case of stationkeeping at ranges comparable to target dimension. 
Candidate ranging waveforms for short range operation are analyzed and 
compared. The logarithmic error discriminant utilized for range, range 
rate and angle tracking is formulated and applied to the quantitative 
analysis of radar subsystem tracking loops. 
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1.0 


INTRODUCTION AND OVERVIEW 


This is the first report for the "Study to Investigate and Evaluate 
Means of Optimizing the Ku-Band Combined Radar/Communication Functions 
for Space Shuttle" (Contract NAS 9-14614, Exhibit D). The study consists 
of analytical tasks directed towards optimizing the design of a baseline 
configuration for the Ku-Band Integrated Radar and Communication Equipment. 
The Ku-band equipment described and analyzed in this report is being 
developed by Hughes Aircraft Company for the Rockwell International /Space 
Division, Downey, California. The equipment contractor is now proceeding 
with its proposed baseline design, but a number of areas of equipment 
function require analytical support to make the performance compatible 
with the requirements [1]. 

Although, as the description implies, the Integrated Equipment 
performs both the radar and the communication functions, the tasks 
described in this report pertain primarily to the radar mode of opera- 
tion. Analysis of various facets of the optimized communication mode 
operation has been described in one of the recent Axiomatix reports [2]. 
Additional tasks pertaining to both the radar and to the communication 
modes of the equipment’s operation, as well as the tasks dealing with the 
integration aspects of the equipment design v/ill be covered during the 
remaining phase of this contract. 

Briefly, the Ku-Band Integrated Radar and Communication Equipment 
time shares the two respective functions during the various phases of the 
Space Shuttle orbital flights. The radar is used for rendezvous with 
orbiting passive targets, as well as with cooperative (transponder-equipped) 
targets. Rendezvous requirements include initial detection of passive 
targets at the maximum range of 12 nmi. The active targets are to be 
detected at a 300 nmi maximum range. 

Subsequent to detection, the radar must track the targets to a 
range as short as 100 feet, where the stationkeeping takes place. Target 
parameters measured by the radar during the tracking phase include: 

(1) range, (2) range rate, (3) angle, and (4) angle rate. The radar is 
of a coherent pulse doppler type with a peak transmitter power of 50 watts 
nominal. A 36- inch diameter center- fed parabolic dish is used as an 
antenna, whose terminals also provide the tv/o-axis monopulse difference 
signals used for angle tracking. Both the transmitter, which is a TWT 
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amplifier, and the antenna with its associated monopulse electronics are 
shared with the communication mode. 

In the communication mode, the Integrated Equipment is used as a 
spaceborne terminal to provide a Ku-band two-way link betv^een the Shuttle 
Orbiter and the ground. This Ku-band link is via the TORS. The forward 
link to the Shuttle consists of a 216 kbps data stream spread by a 
3.03 Mcps PN code to reduce the power density incident upon the earth. 

The return link consists of three independent data streams phase- 
multiplexed for transmission on a single carrier. The data rates of the 
three channels range from about 200 kbps for the lov/est rate channel to 
as high as 50 Mbps for the wideband channel. 

The forward link communication signal is, received at the Shuttle 
on the nominal frequency of 13,775 GHz. The return link communication 
transmission to the TORS is on a 15,003 MHz carrie'^. The radar operates 
in the approximate vicinity of the conmunication receive frequency. This 
permits time sharing of the receiver and the monopulse electronics betv/een 
the two functions. The radar operating frequency is hopped over approxi- 
mately 200 MHz to reduce target scintillation. 

The baseline design proposed for the Ku-Band Integrated Radar and 
Communication Equipment by Hughes [3] is based upon a subdivision of the 
equipment into five subunits which are: 

(1) Deployed Mechanism Assembly (DMA) 

(2) Deployed Electronics Assembly (DEA) 

(3) Electronics Assembly 1 (EA-1) 

(4) Electronics Assembly 2 (EA-2) 

(5) Signal Processor Assembly (SPA). 

The deployed mechanism assembly (DMA) includes the antenna reflector, 
feed, antenna gimbals, drive motors, gyros, digital shaft encoders, rotary 
joints and connecting v/aveguides, and associated wires and cables. 

The deployed electronics assembly (DEA) consists of a common 
receiver and transmitter, used for both radar and communication, and a 
common crystal-controlled exciter/LO generator, also shared by the two 
functions. The interconnection between the DEA and the DMA is via two 
rotary joints and cable wraps. 

Both the DMA and the DEA are located within the Shuttle payload 
bay and thus have to be deployed for their operation upon the opening 
of the payload bay doors. 
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The communications electronics assembly 1 (EA-1) includes the PN 
despreader, the Costas loop phase-lock demodulator, the antenna servo 
control electronics and motor drive amplifiers, the data multiplexer/ 
demultiplexer, and input/output buffers. As a subunit of the communica- 
tion receiving equipment, the EA-1 also contributes to the TORS signal 
acquisition and forward link data detection. 

The electronics assembly 2 (EA-2) consists of the radar second 
mixer, second IF amplifier, I/Q detectors, digital Fourier transform 
(DFT) doppler processor, radar angle tracking circuits, radar range track- 
ing circuits, radar range rate tracking circuits, radar timing and control, 
and the radar AGC detector. 

The signal processing assembly (SPA) includes management logic, 
fomard data processing, and return link data multiplexing. 

Both of the electronics assemblies, as well as the signal processor 
assembly, are located in the avionics bay of the Shuttle. The interconnec- 
tion with the two deployed assemblies is via wires and coaxial cables. 

Of the five assemblies listed above, the functions performed by 
the EA-2 are given main consideration in this report. The functions of 
other assemblies are described only where necessary to explain, equipment 
operation in the radar mode. 

Section 2.0 presents the description of the radar system as con- 
figured by the appropriate assemblies. The summary of power budgets for 
passive radar target detection is given in Section 3.0, with the support- 
ing analysis contained in Appendix A. 

Section 4.0 summarizes radar search and track mode operations, with 
detailed description given in Appendix B. 

Time-multiplexed single-channel angle tracking of passive targets 
is described briefly in Section 5.0 and in detail in Appendix C. Section 
6.0 highlights the radar performance in the presence of main lobe ground 
clutter with supporting details given in Appendix D. The problems asso- 
ciated with stationkeeping and short range operation are addressed in 
Section 7.0. Ranging waveforms for short range operation are summarized 
in Section 8.0, with details in Appendix E. Section 9.0 describes the 
maximum likelihood estimator used for the range, range rate, angle, and 
angle rate estimation. Supporting analysis is given in Appendix F. 

Conclusions and recommendations are presented in Section 10.0. 
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2.0 RADAR SYSTEM DESCRIPTION 

The functional block diagram for the Integrated Ku-Band Radar/ 
Communication Equipment operating in the "Radar" mode is shown in Figure 1. 
The diagram is presented in a manner which reveals the signal flow inter- 
action between the assemblies utilized for the radar operation. Subunits 
utilized for the communication function are shown in this -diagram only if 
they are time-shared with the radar function. 

The four subassemblies of the Ku-Band Integrated Equipment shown 
in Figure 1 are; Deployed Electronics Assembly (DEA), {2) Deployed 
Mechanism Assembly (DMA), (3) Electronics Assembly 2 (EA-2), and (4) 
Electronics Assembly 1 (EA-1). The function of each assembly, as well 
as that of the overall system, can be understood best by following the 
signal flow. 

The coherent frequency synthesizer located within the DEA provides 
all the RF, LO, and IF signals required to generate, transmit, receive 
and demodulate the radar signals. The radar signal originates in the 
Radar Modulator located within the DEA. The modulator accepts a 651 MHz 
CW signal from the synthesizer and keys it in accordance with a selected 
transmit/receive T/R pulse modulation format. The latter is developed 
by the timing and control unit within the EA-2. The pulse modulated 
651 MHz signal is then up-converted to the 13.8 GHz (approximate) radar 
operating band by mixing with the five hopped frequencies of the trans- 
lation oscillator (TO) signal. This upper sideband translation is shov/n 
in Table 1 * 


Table 1. Radar Transmit Frequencies Synthesis 


T.O. 

Frequencies 

(MHz) 

Radar Modulator 
Center Frequency 
(MHz) 

Up-Converter 

Output 

(Tx Frequencies) 
(MHz) 

13,128 N 


/ 13,779 

13,180 / 


\ 13,831 

13,232 > 

Add 651 

^ 13,883 

13,284 \ 


/ 13,935 

13,336 ^ 


^ 13,987 


Radar Frequencies 
(Upper Sideband) 


In the communication mode, only the 13,128 MHz TO signal is mixed 
with the 1,875 MHz communication modulator carrier, resulting in a single 
transmit center frequency of 15,003 MHz (15.003 GHz). 
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As shown in Table 1 , the radar frequencies are spaced 52 MHz 
apart, resulting in a total diversity bandwidth AF of 208 MHz. 

The output of the up-converter is applied to a preamplifier, 
which consists of a high frequency FET device. The preamplifier, in 
addition to supplying the gain, provides limiting prior to the TWT. 

This limiting minimizes the AM-to-PM conversion of the transmitted signal 
passing through the TWT. The power output P-j of the FET limiting pre- 
amplifier unit is about +7 dBm. 

In the search mode and the "long range" (over 3000 feet) track 
mode, the radar signal is amplified by the TWT. This amplification 
brings the drive power to the nominal 50 v/atts peak output power Pg- 
The tube gain assumed for this case is aq dB. The function of the vari- 
able attenuator is to reduce the TWT output pov/er for the short-range 
operation. 

The output of the TWT is also gated in accordance with the T/R 
signal. Such gating prevents the tube noise from interfering with the 
low level signals received during the search and the long range tracking 
operations. Note that the TWT output is delivered to the High Power 
Pulser by means of Diplexer 1. This diplexer restricts the gating func- 
tion only to radar region of the TWT output signal. Diplexer 2 re-inserts 
the gated radar signal into the main transmitter output RF line. The 
latter connects the transmitted signals (both radar and comm), via a 
circulator and a rotary joint, to the main antenna. Switch 3 is used 
only in the communications mode. For the initial acquisition of the 
TORS, the TWT output is radiated via the auxiliary horn to facilitate 
spatial search for the TDRS signal by the TORS. 

In the short range radar mode, the TWT amplifier is bypassed 
and the preamplifier output is applied directly to Diplexer 2. Switch 1 
is used to gate the low level signal and Svn'tch 2 is in the SR position. 
The gated low level signal, like the TWT output, is applied to the main 
antenna via the circulator and the rotary joints. The position of 
Sv/itch 2 therefore determines whether the TWT or the preamplifier output 
is applied directly to the main antenna. 

The main antenna, which is a 36- inch diameter, prime- fed parabola, 
is a part of the deployed mechanism assembly (DMA). The antenna uses a 
five-element monopulse feed. Four elements are used to generate the aAZ 
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and AEL errors in the monopulse bridge. The fifth element is used for 
transmission and reception of either the radar or the communication signals. 
During the reception, the fifth antenna element provides the sum (s) 
channel signal. The polarization of the main antenna is selectable. For 
tne communication mode, it is circular (RHC) for both transmit and receive. 
During the radar mode, the polarization is linear and transmission is 
only via the main antenna. 

The reception of the signal returned by either a passive or an 
active target is performed by both the main parabolic and auxiliary horn 
antennas. The purpose of the horn reception is to guard against the side- 
lobe lock-ons. The main antenna provides the sum (e) signal to the 
receiver. The monopulse bridge provides the two difference (a) signals, 
i.e., aAZ and AEL. The sum channel output R (s) passes through rotary 
joint 1, the circulator, the transmit/ receive (T/R) sv/itch, and terminates 
at the input of the E-channel first mixer, M(e). 

The two aAZ and AEL signals are time-multiplexed into a cotmion 
channel, henceforth referred to simply as the A-channel. During the 
spatial search for the radar target, the receiver channel, which is used 
exclusively for A-signal amplification during the track mode, amplifies 
the output of the auxiliary horn. This is accomplished with Switch 4 in 
the "Guard" position. The guard signal passes through rotary joint 2, 

Switch 5, the T/R switch, and terminates at the input of the A-channel 
receiver first mixer M(a). 

The first mixers translate their respective Ku-band input signals 
to the first IF frequency. In the radar mode, the first LO signal 
applied to these mixers is the same frequency-hopped signal which is used 
for the up-conversion in the radar transmitter. The incoming frequency- 
hopped signal is therefore de-hopped in the manner indicated in Table 2. 

For monopulse tracki ng in both the radar and communication modes, 
a portion of the E-channel signal is applied to the A-channel, thus form- 
ing a a±aAZ±aEL signal. As shown in Figure 1, this summation is per- 
formed at the first IF. Because the guard channel is not required in the 
tracking mode. Switch 6 remains in the e + a position during tracking. 

During the search phase, the guard channel has to be sampled 
periodically to prevent the lock-up onto the sidelobes. The sampled 
guard channel output is then compared to the main lobe, i.e., the sum 
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channel signal. Alternating the position of Switch 6 between the s±A 
and the GD (Guard) terminals permits time-sharing of the single second 
IF channel of the radar receiver. It must also be noted that, during 
search, the A-signals are not used and, consequently, the signal v;hich 
is 2±a during tracking is actually only Z during acquisition. Further- 
more, as explained in Section 4.0, the tracking of salient radar measure- 
ment parameters such as range, velocity (range rate), and angle are per- 
formed on the e + aAZ ± aEL signal, which is the only signal supplied to 
the third IF of electronics assembly 2 (EA-2).* 

Table 2. First Conversion Frequency ne-Hopping 
(Low Side LO Injection) 

Radar Frequency First LO Frequency 

(MHz) (MHz) 

13,779 13,128 

13,831 13,180 

13,883 13,232 

13,935 13,284 

13,987 13,336 

'if 

For communication mode only, the 13,128 MHz first LO fre- 
quency is used. The translation is then: 13,775- 13,128 

= 647 MHz for the first communication IF frequency. 


First Radar 
IF Frequency* 
(MHz) 


651 


The second IF radar signal is amplified and applied to a coherent 
processor which consists of in-phase (I) and quadrature (Q) detectors. 

The analog outputs of these detectors are low-pass filtered and then 
converted to digital data which contains both the amplitude and phase 
information of the received signal. 

The subsequent processing of the radar signal is all digital. 

The processing consists of presum, digital Fourier transform (OFT), post- 
detection integration (PDI), and logarithmic discriminant formatter. 
During the search mode, the output of the PDI is applied to the threshold 
detector which stops the antenna search pattern when the target is 


For the communication mode only, the E-channel is translated to 
the third IF located in EA-1. The recovery of the angle tracking infor- 
mation is performed by envelope detection of the e±aAZ±aEL signal at 
the first IF and by subsequent demultiplexing. 
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detected. Subsequently, the angle and range tracking units take over 
and the system is sv/itched from the "track" mode to the "search" mode. 
Descriptions of the search and track modes are presented in Section 4.0 
and Appendix B of this report. 
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3.0 POWER BUDGET FOR RADAR PASSIVE TARGET DETECTION 

2 

The initial detection of the passive 1 m target at the 12 nmi 
range is the driving requirement for the radar system's power budget. 
Appendix A of this report presents a detailed description of the pertinent 
factors contributing to the power budget for the Ku-band radar in the 
passive target mode. Based on these factors, power budgets are calculated. 
The following two specific cases are considered: (1) the improvement due 

to the use of frequency diversity is negligible (slow fading), and (2) 
the improvement due to the use of frequency diversity is ideal (fast 
fading). 

Table 3 shov/s the power budgets for these tv/o cases. As can be 
seen from this table, the design margin for the slow fading case is nega- 
tive, specifically -6 dB. In comparison, the fast fading case assumption 
results in a positive design margin of +2 dB. The latter is a more real- 
istic case. 

An alternate set of power budgets, which incorporate recent changes 
in system parameters, is also considered. Although these changes were 
incorporated into the system design to decrease implementation complexity, 
the penalty paid turns out to be the performance degradation for the 
passive search mode. Table 4 shows the updated power budgets. As is 
shov;n in Table 4, the design margin for the slow fading case is lowered 
to -8.5 dB and the margin for the fast fading case is reduced to -0.54 dB. 

Conservative assumptions used for the updated povjer budgets may 
account for the negative margins. It is therefore anticipated that per- 
formance in some cases will be improved and the design margins will be 
increased accordingly over those presented here. There are various alter- 
natives to improve the design margin to a more satisfactory value, namely, 

6 to 10 dB for the passive mode. 

In all cases considered, the gain of the main antenna v/as assumed 
to be G=38.5 dB, which is the most recently known value [4]. Further- 
more, all of the budget estimates presented in this report assume that 
two scans of the volume to be illuminated can be accomplished in 1 minute. 
At 12 nmi, the two-sided cone angle that must be searched is 40 degrees. 

It has been shown in [5,6] that, at these values of cone angle and range, 
two scans can be completed in 1 minute with an antenna spatial dv/ell time 
of 36 msec. This value of spatial dv/ell time is used in this report for 
pov/er budget estimation. 
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Tabic 3. Power Budget for Passive Search at 12 rmii llUliout Range Designated 


ltf3m 

Term 

Section 
Humber 
in This 
Report 

Value 

Budget /^1 
Case 

(Slov; Fade) 
(dB) 

Budget 112 
Case i<2 
(Fast Fade) 
(dB) 

1 

% 

3.1 

Required Peal: SRR 

18.2 

11.7 

2 


3.4 

G=38.5dB 

-77.0 

-77.0 

3 


3.3 

X*= 0,0218 m 

33.24 

33.24 

A 

0 

3.2 

1 

0 

0 

5 




33 

33 

6 


3.6 

R= 12nmi =22,2 km 

173.9 

173.9 

7 

k 

3.7 

1. 38x10"^^ 

-228.6 

-228.6 


Ts 

3.7 

1500®K 

31.76 

31.76 

9 

Coherent Addition 
Time 

Losses 

3.5 

1.195 msec 

29,23 

29.23 

10 

2.0 

3.8 

2 

-3 

-3 

n 

Transnnt Losses 

3.8.1 

RFP 

3.5 

! 

3.5 

12 

1 

Beam Sliape 

3.8,2 

Scan Alignment- 1 

Lateral Scan 1 

3.2 

3.2 

13 1 

Threshold 

3.8.3 

Constant FAP 

1.0 

0,5 


Processor Loss 

3.8.4 

1 

1.7 

1.0 

15 

Range Gate 
Straddling Loss 

3.8.5 


0 

0 

16 

Doppler 

3.8.6 . 

Kismatch 

1.0 

0,7 

17 

PDl 


Postdetection Inte- 
gration BEF [5] 

1.25 

T.25 

18 

Pre-Sum Mismatch 

3.8.7 

Due to Doppler 

0.6 

0,6 

19 

% 

3.5 

Required Peak Power 

22.98 dB 
{198.6 w) 

14.98 dB 
(31,5 w) 

20 

't 

3.5 

Duty Factor =0.2 

-7.0 

-7.0 

21 

’’avg 


Required Average 
Power 

15.99 dB 
(39.7 w) 

7.99 dB 
(6.3 w) 

22 

Pjj Design 


Reference [ 5 ] 

50 w 

50 w 

23 

Pavg 


Reference [5l 

10 w 

10 w 

2^ 

Design Margin 

■ 


-6.0 dB 

+2.0 dB 
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Table Alternate Pov/er budget for Passive Search at IZ nmi VMtliout Range Designated 


Item 

Term 

Section 
Number 
in This 
Report 

Value 

Budget //3 
Case frl 
(Slow rode) 
(dll) 

Budget //4 
Case 112 
(fast Fode) 
(dB) 

1 

•"p 

3.1 

Required Peak StIR 

18.2 

11.7 

2 

g2 

3.4 

G=38.5 dB 

-77.0 

-77.0 

3 

x2 

3.3 

X» 0.02161 m 

33.31 

33.31 •• 

4 

5 

3*2 

ln.2 

0 

0 

5 

(4i.)^ 



33 

33 

6 


3.6 

R= 12 nmi = 22.2 Lm 

173.9 

173.9 

7 

k 

3.7 

1.38x10"22 

-228.6 

-228.6 

8 


3.7 

1500“K 

31.76 

31.76 

9 

Coherent Addition 
Time 

3,5 

1.062 msec 

29.74 

29.74 

10 


2,0 

2 

-3 

-3 


Losses 

3*0 




11 

Transmit Losses 

3.8.1 

Rrr 

3.5 

3.5 

12 

Beam Shape 

3.8,2 

Scan Alignment^ 
Lateral Scan 

3.2 

1 

3.2 

13 

Threshold 

3.8.3 

Constant FAP 

1.0 

0,5 

14 

Processor Loss 

3.8.4 


1.7 

1.0 

15 

Range Gate 
Straddling Loss 

3.8.5 


2.0 

2.0 

16 

Doppler 

3.8.6 

Mismatch 

1.0 

0.7 

17 

PDI 


Postdotection Inte- 
gration Ref. Cs] 

' 1.25 

1.25 

18 

Pre-Sum I'llsmatch 

3.8.7 

Due to Doppler 

0.57 

0.57 

19 

'’P 

3.5 

Required Peak Power 

25.53 dB 
(357.3 w) 

17.53 dB 
(56.G w) 

20 


3.5 

■Duty Fcctor=0.2 

-7.0 

-7.0 

21 

’’avg 


Required Average 
Power 

18.54 dB 
(71.5 w) 

10.54 dB 
(11.3 w) 

22 

Pp Design 


Reference [s] 

50 W 

50 W 

23 

Pavg Design 


Reference [5] 

10 Vf 

10 w 

24 

Design Margin 



-8.5 dB 

-0.54 dB 
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4.0 RADAR SEARCH AND TRACK MODES 

The objective of the Ku-band radar system on the Shuttle Orbiter 
is to detect the presence of a target and then obtain continous accurate 
estimates of the various target parameters, namely, its range, range rate 
(velocity), azimuth angle, elevation angle, and angle rates. The target 
is first detected when the radar is in the search mode; then the accurate 
estimates are obtained in the track mode. 

The radar employs coherent processing over each RF frequency 
transmitted and noncoherent postdetection integration (PDI). Follov/ing 
the intermediate amplifiers, in-phase and quadrature phase signals are 
sampled and digitized (A/D converted), providing target amplitude and 
phase information. For longer range search and acquisition (R>3000 ft), 

16 doppler filters are fonned, covering the doppler interval defined by 
the repetition rate. This is an ambiguous doppler interval; this ambi- 
guity is removed via differentiation of the range estimates. 

Following the filters, the target magnitude is determined, summed 
over the number of RF frequencies (noncoherent postdetection integration) 
and compared to a threshold. The threshold is set from noise measure- 
ments, so as to maintain a constant false alarm probability (CFAP). This 
can be implemented by any of several methods. 

In the search mode, only the sum channel data is processed. V/hen 
the threshold is exceeded and target detection is declared, the auxiliary 
antenna signal is processed and its magnitude is compared with that of 
the main antenna to eliminate sidelobe detected targets. The auxiliary 
antenna has a peak gain which is approximately 20 dB less than that of 
the main antenna. Also, the sidelobes of the main antenna are approxi- 
mately 20 dB down from the main lobe. This provides roughly a 20 dB 
main lobe/guard antenna ratio to detect and eliminate sidelobe targets. 

In the modes for which angle tracking is required, namely, GPC 
acquisition and autotrack, the sum-plus-angle error channels are processed. 
The azimuth and elevation error signals are time-multiplexed, thereby 
eliminating the need for a second and third matched processing channel. 

The Shuttle Ku-band radar signal flow block diagram is shown in 
Figure 2. The antenna assembly consists of a main antenna and an auxiliary 
antenna, which are subunits of the deployed mechanism assembly (DMA). 
Further processing of antenna signals takes place in the deployed 
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Figure 2. Radar Block Diagram 
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electronics assembly (DEA), where various switching and filtering opera- 
tions take place. The antenna assembly is followed by several IF mixing 
stages to bring the carrier frequency down to a stable IF of fp= 78.143 
f-IHz. The output of the I-Q stage that follov/s is a pair of baseband sig- 
nals which are in phase and in quadrature with the input signal. 

The available information is converted by an A/D converter to 
digital data which samples both channels at a rate of fg = 480 kHz, which 
is slightly larger than the Nyquist rate of the baseband signals at the 
output of the I-Q stage. (The IFF bandwidth in both the I and Q channels 
is equal to 237 kHz.) The output of the A/D converter is a set of complex 
numbers which undergo various stages of digital processing before esti- 
mates of the target parameters (range, velocity, azimuth angle, elevation 
angle) are obtained. The first processing stage is a presumming stage, 
followed by a set of range gates and doppler filters consecutively. The 
arrangement of the gates is varied for the two modes of operation. In 
the search mode, vjhere range designation is not present, there are four 
range gates per pulse, while in the track mode there are only tv/o non- 
overlapping gates: an early gate and a late gate. 

For long ranges, the magnitudes of the outputs of the doppler 
filters are summed over the RF frequencies and compared to a threshold 
to detect the existence of a target. For short ranges (R< 3000 ft), the 
detection takes place v/ithout using any doppler filtering. 

Once a target is detected, the radar is sv/itched to a track mode 
and the digital signals undergo postdetection integration, followed by 
an error estimation unit which is initiated by a crude estimate obtained 
by the search mode operation. The final estimates of the various param- 
eters are: 

(1) Fed as inputs to the servo motors for angle track. 

(2) Used to adjust the position of the early and late gates 
for range track. 

(3) Displayed for the astronauts for visual reading (all 
variables being tracked). 

(4) Utilized to update the state of the error estimators to 
obtain new estimates. 

Frequency diversity, using five RF frequencies, is employed in 
both the search and track modes. Table 1 (page 4) lists these frequencies. 
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The radar modulation parameters are summarized In Table 5. 


Table 5. Radar Modulation Parameters as a Function of Range 


Search Mode 


Designated 

Range 

(nmi) 

T 

(ysec) 

PRF 

(Hz) 

Frequency 
Dwell Time 
(msec) 

Spatial 
Dwell Time 
(msec) 

8-22 

66.3 

2987 

5.3 

31.8 

4-8 

33.2 

7177 

2.2 

13.4 

2-4 

16.6 

7177 

2.2 

13.4 

1-2 

8.3 

7177 

■ 2.2 

13.4 

0.5-1 

4.1 

7177 

2.2 

13.4 


In explaining the search mode operation, low range signal format 
Is chosen for demonstration. All other ranges follow identical processing 
except for the short range (<0.5 nmi) v/hich does not utilize the doppler 
filters. 

Sixteen pulses are transmitted at each RF frequency. The duration 
of each pulse is 66.3 ysec, while the time betv/een consecutive pulses is 
335 ysec, as shov/n in Figure 3. Four range gates (R^ through R^) are ■ 
used to cover the designated range. It is important to note that no dop- 
pler correction is used at the presuminer in the search mode. For each 
RF frequency, 64 complex numbers are stored in a memory designated at 
M(l), as shown in Figure 4. The output of each range gate at every RF 
frequency is passed through a bank of 16 doppler filters, implemented as 
discrete Fourier transforms (OFT) or fast Fourier transforms (FFT), and 
the magnitudes of the outputs of these filters are calculated and compared 
to precalculated thresholds. The target is detected when the magnitude 
of one or two outputs of adjacent doppler filters pertaining to one or 
two range gates is exceeded. An initial estimate of the range and the 
range rate of the target are calculated from the knowledge of the doppler 
filters and range gates whose outputs have exceeded the threshold. These 
initial estimates are fed into the error estimators as initial conditions 
to start the successive estimation process in the track mode. Thresholds 
are calculated to produce an error alarm rate of one false alarm per hour. 
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Figure 5 Illustrates the signal format in the track mode. As 
shown in the figure, the signal in the track mode consists of 64 pulses 
for each RF frequency, instead of the 16 pulses of the search mode. For 
long range tracking, the frequency dwell time is equal to 21.44 msec. 

Each frequency dv/ell time is divided into four time slots by the sv/itch- 
ing combination of SWl and the phase encoder switch in the antenna assembly. 
The signals in the succeeding time slots are E + AAZ, i-aAZ, e + aAZ, and 
E-AEL, respectively. The time slots are arranged so that the first two 
are used to measure the azimuth angle, while the last two are used for 
the elevation angle. Each one of the previous time slots consists of 
16 pulses. The pulses are passed through an early gate and a late gate, 
which are located around a predicted estimate of the range, Rp(n), from 
the processor timing unit. The width of each gate is equal to the pulse 
vndth. The output of each range gate is a complex number z(I,J,L,K), 
which is a function of four parameters (I, J, L, K): 

I denotes the RF frequency being transmitted; I =1,2,..., 5. 

J denotes the time slot in each frequency such that 
J = 1 corresponds to s + aAZ 
J = 2 corresponds to x - aAZ 
J = 3 corresponds to 2 + AEL 
J = 4 corresponds to 2 -AEL. 

K denotes a particular pulse in each time slot; K=l,...,16. 

L denotes the range gate under consideration: 

L=-l corresponds to Early Gate 
L = +l corresponds to Late Gate. 

The detailed description of signal processing involved in the 
search and track modes is presented in Appendix B. 
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Figure 5. Track Mode Signal Format 
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5.0 TIME-MULTIPLEXED SINGLE-CHANNEL ANGLE TRACKING 

OF PASSIVE TARGETS 

In the Ku-band integrated radar and communications for the Orbiter, 
the angle tracking subsystem, when operating in the communications mode, 
has received extensive attention [7]. In Appendix C of this report, the 
radar angle tracking of passive targets is studied. Both thermal noise 
and target scintillation are taken into account. It is shown that the 
performance requirements for the angle tracking of passive targets meets 
the specification requirements with a design margin of greater than 10 dB 
at 10 nmi and closer. The angle rate estimation requirements were not 
included in this report. 

A block diagram of the single channel* monopulse angle tracking 
loop is shown in Figure 6. All of the operations pertinent to the angle 
tracking functions are as shown in this figure. In particular, the pre- 
sum operation assumes that the range gate tracking loop is also tracking 
satisfactorily. Degradations in signal-to-noise ratio due to errors in 
the range gate tracking loop would not seriously affect the operation 
of the angle tracking loop unless these degradations were very signifi- 
cant. This is because of the ample design margin of the angle tracking 
loops. 

The same comments can also be made about the doppler tracking 
loop, which sets the centers of the doppler filters. 

Th.e performance of the angle tracking loop is shov/n in Figure 7, 
where the RMS tracking error due to thermal noise and target scintilla- 
tion is plotted against range in nautical miles. For all ranges less 
than 18 nmi, the RMS error is below the specification value of 0.153 deg, 
and naturally improves as the range decreases. 

In summary, there appears to be no difficulty in angle tracking 
capability of the radar for passive targets from the standpoint of the 
thermal noise. The passive mode was chosen for this detailed study since 
it imposes the most severe requirements on signal-to-noise ratio. 


* 

Single channels implies multiplexing of three information 
channels (e, AAZ and AEL) into a single channel. 
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Figure 6- Block Diagram of Single Channel Monopulse Angle Tracking Loop 
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6.0 RADAR PERFORMANCE IN PRESENCE OF MAIN LOBE 

GROUND CLUTTER 

The possibility of main lobe ground clutter interference arises 
when the Shuttle radar is working with a target which is in a position 
requiring the radar antenna to be pointed in the direction of the earth 
and close to vertical. The geometry characteristic of this case is shown 
in Figure 8. In addition to the geometrical causes for ground interfer- 
ence, the signal structure utilized by the radar may be such that clutter 
discrimination is not easily implemented. For example, a radar system 
whose parameters are optimized for operation at ranges of 1'2 nmi or less, 
such as the Shuttle Ku-band system, will typically employ PRF in the 
range of 3000 to 8000 kHz, i.e., just low enough to provide nonambiguous 
acquisition of range. Such high PRFs, however, will cause the ground 
return at orbit altitudes ranging from 100 nmi to 500 nmi to appear as 
a target “blip" within the same range interval as the desired target. 

In other words, with PRFs matched to target range (specifically passive 
target) clutter discrimination in time is not possible. 

Furthermore, the magnitude of the clutter return itself, at a 
minimum orbit altitude of 100 nmi, may be such that the passive target 
signal is masked by clutter at ranges far shorter than the maximum speci- 
fied acquisition range of 12 nmi. 

The detailed quantitative analysis of the ground clutter return 
is given in Appendix D. It is shown there that, at vertical incidence, 
the ground return is of a beamwidth limited type for all pulse widths 
greater than 1 ysec and for orbit altitudes ranging from 100 nmi to 
500 nmi. Figure 9 illustrates this fact. The nature of beamwidth limi- 
tation is evident by the fact that the ground return, even for h= 500 nmi, 
approaches its asymptotic maximum value vn'th “response times" less than 
1 ysec,* Because the long range tracking and acquisition mode pulse 
length .of the system are above 2 ysec beamwidth, limited operation is 
assumed for the case of main lobe nadir return.) 

Utilizing the analytical techniques described in Appendix 0, we 
have calculated the relative power (in dBm) of the desired passive target 


Response time is arbitrarily defined here as the time required 
for the main lobe nadir clutter return to reach 9D% of its maximum power 
value. 


Shuttle 

Orbiter 



Figure 8. Rendezvous Geometry* With Ground Clutter Return iti Main Lobe 
Not to scale 
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return(s), Pg, an?i the power of the nadir clutter return at various alti- 
tudes. Figure 10 presents the pertinent information. Note that all pov/er 
levels shown are referred to presum output. This reference point is 
chosen because beyond it analytical comparison becomes nearly impossible, 
due to extremely complex time-frequency relationship of the ground clutter 
return. 

From Figure 10, it is evident that, if one accepts the criteria 

of Pg/Pj,^20 dB (at presum output) and all radar modes functional, the 

system range at h =100 nmi is limited to 1 Ion for a 1 m^ target and to 

2 

about 1.8 km for a 10 m target. This is a significant reduction from 

2 

the specified maximum range of 12 nmi (22.2 km) for a 1 m targe.. Note, 
however, that at 0.5 nmi (0.927 km) and with a =10 m , parameters typical 
of a Skylab reboost mission, P 5 /PQ = 30 dB, even at h=100 nmi. This 
makes ‘the radar system under consideration acceptable for the Skylab 
reboost mission without requiring special modifications or parameter 
changes. 

From the standpoint of improving the performance of the Ku-band 
radar system for targets other than the Skylab, one should consider 
"fixes" available to the system. One such fix, which is actually a 
built-in one, is the frequency agility employed by the radar for com- 
bating target fluctuations. Figure 11 depicts the altitude versus inci- 
dence angle diagram of the clutter interference region for the two PRFs 
currently employed by the Ku-band radar for the search and track modes. 
According to Figure 11, PRFs of 7177 kps and 2987 kps result in “clutter- 
free" orbit altitudes of 192 nmi and 460 nmi, respectively. The 192 nmi 
orbit may be acceptable for the Skylab reboost mission, but the 460 nmi 
orbit altitude, although within the 500 nmi maximum specification, may 
be too high to be associated v/ith clutter-related problems. Nevertheless, 
the intrinsic ability of the Ku-band system to reject clutter by frequency 
agility may be of benefit to some of the Shuttle missions. 

Suppose now that one considers a frequency agility pattern which 
eliminates first-time ground returns up to 500 nmi. The region of direct 
(i.e,, first-time-around) interference would then be as shown in Figure 12. 
This is a highly desirable result, because it provides a "clear" operation 
at all altitudes of up to 500 nmi. At 500 nmi, however, the ground return 
will be at least 14 dB down from the 100 nmi return. The penalty paid 
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Power Received From a :“issive Target Relative to 
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power refers to pre-suir output) 
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for achieving the clutter-free operation at altitudes between 100 nmi 
and 500 nmi is the shortening of the coherent frequency dwell time. For 
a PRF of 2987 kps, this means that only about 4 pulses are available for 
coherent integration at any one frequency. 

Thus, a 16- point Fourier transform cannot be implemented with 
"fast" hops {1.24 msec) commensurate with the 100 nmi ground return time. 

On the other hand, if one considers a possibility of coherently summing 
consecutive 4-pulse groups to form a 16-point transform, the clutter- 
free region operation depicted by Figure 12 may be achieved. The factors 
to cor.ider are: (1) transmitter coherency over several discrete frequency 

dwells, (2) target coherency over these dwells, and (3) the effect of 
time-multiplexed F-transformation on the memory storage requirements. 
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7.0 STATIONKEEPING CONSIDERATIONS FOR THE KU-BAND RADAR 

One of the most difficult requirements for the Ku-band radar is 
the capability of establishing a stationkeeping environment for the 
target at approximately 35 feet from the Orbiter. Because of the opera- 
tional planning of the Skylab reboost, the Ku-band radar will have to be 
able to maintain stationkeeping when looking into the earth clutter 
background. 

As the operation is presently envisioned [8], there will be a 
ground-controlled sequence of maneuvers to get the Orbiter within 300 nmi 
of an active target and within 10-20 nmi of a passive target. The Ku-band 
radar will then acquire and track the target down to a braking or initial 
stationkeeping region which will be at a range of 1000-3000 feet. The 
final approach will then be carried out to a range of 35 feet, where it 
can be retrieved by the 50-foot boom. 

The geometry of the braking, initial stationkeeping, and final 
approach is shovm in Figure 13, which is similar to that shov/n in [8]. 

Four positions are shown. Although it is not pertinent to the radar con- 
siderations, from an energy consumption point of view, stationkeeping 
will be significantly easier in positions V-j and Vg than in positions Rq 
and R2* What is important from the radar point of view is that some 
targets (primarily military) will be at low altitudes, thereby possibly 
requiring an approach from above, with an initial stationkeeping position 
at R2. In this position, the earth clutter problem must be addressed. 

It is envisioned also that the requirements on range and range rate esti- 
mates v/ill be more accurate in the initial stationkeeping positions 
and 82* 

The radar operation at these close ranges have been broken down 
into three modes; 

(a) Short mode 300 to 3000 feet 

(b) Short-short mode 100 to 300 feet 

(c) Short- short- short mode 35 to 100 feet 

Whether these modes can be consolidated v/ill depend on which implementa- 
tion is chosen. Preliminary comments concerning the choice of implemen- 
tation are given below. The approach taken is as follows: 

(1) What needs to be done to modify the existing radar structure 
so as to meet all. requirements in the Short mode of operation without an 
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Figure 13. Geometry of Approach and Braking 
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undue amount of additional hardware? 

(2) If a separate unit is employed for the Short modes, what 
are the candidates and which positions of the existing RF sections could 
also be used by this unit? 

It is our understanding that search and track capability is 
required for the Short-Short mode and that it is highly desirable for 
the Short-Short-Short mode. One of the first problems is to evaluate 
exactly what it means to track in angle, range and range rate at ranges 
as short as 35 feet. Figure 14 illustrates the point. With a 36-inch 
(D=0.91 meter) diameter antenna, the near field region extends to about 
76 meters, or about 250 feet. This means that, at the distance of 35 feet, 
the main antenna beams are not well defined and the monopulse capability 
must be questioned. The target reference point to vfhich the radar param- 
eters (range, range rate, angle and angle rate) are to be measured also 
needs to be established. 

On the positive side, signal-to-noise ratio should not be a prob- 
lem. In fact, to protect targets from high levels of electromagnetic 
radiation, the TWT may be bypassed, thereby obtaining a 44 dB reduction 
in transmitted peak power. In addition, to protect the receiver from 
saturation, additional attenuators (ferrite attenuators have been proposed) 
can be inserted in the receiver front end. Implementation difficulties 
center around the rise and fall times of the RF switches, particularly 
since the pulse widths will be between 30 and 122 msec. 

For the existing 36- inch Ku-band antenna, as it applies to the 
Short range modes of operation, it has already been pointed out that the 
near field extends to approximately 76 meters. Therefore, the Short- 
Short and Short-Short-Short modes are entirely within the near field of 
the existing antenna. This probably precludes any monopulse capability 
in these modes. 

One fix is to defocus the antenna at short ranges so that the 
target would be entirely within the main beam. The present antenna cannot 
be defocused sufficiently to permit an 18-meter target to be entirely 
vnthin the main beam. A dipole antenna \i»ould provide a sufficiently vade 
beam, but vrauld not provide any monopulse capability. 

A candidate vfhich could provide monopulse capability is to use 
the monopulse feed of the existing antenna, either by physically rotating 



Figure 14. Antenna Directivities for Radar Tracking of a Passive 
Target at Stationkeeping Ranges 
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the horn by 180®, as shovm in Figure 15a, or to switch to a second horn 
mounted behind the main horn, as shown in Figure 15b. The antenna pat- 
terns of these feeds are ±30° or, equivalently, a 60° beamwidth. In 
Figure 15a, the concept is conveyed of using the same horn for all modes 
of operation, but rotating it 180°, thus eliminating the dish for short 
range operation. Figure 15b illustrates an equivalent fix, only using * 
two horns, switching from one to the other for short range operation. 
This, of course, also bypasses the dish. With a 60° beamwidth, even an 
18-meter target can fit within the main beam out to a range of 18 meters. 
The gain drops to 3 dB, but that appears to be adequate at these short 
ranges. 

It is recommended that this and other modifications of the exist- 
ing antenna subsystem be examined as candidate solutions for the short 
range modes. 

At these short ranges {as low as 35 feet), in order to transmit 
and receive a complete pulse, the pulse must be less than 70 nsec. In 
order to compensate for switching times and finite rise and fall times, 
the pulse v/idth must be substantially less than 70 nsec. If the trans- 
mitted pulse is kept larger than that, then the receiver gate can be 
opened only for the last part of the pulse. This suggests that the 
trailing edge be used for range tracking. The problem is more intricate 
than this, however, since for an 18-meter target, a pulse of any width 
substantially less than 120 nsec will be spread in time to at least 
120 nsec, due to time dispersion by the target. For a point target, 
of course, the pulse will return with little or no time dispersion. 

For large targets, therefore, the concept of range from a tracking point 
of view is not so immediate. The most satisfying range from an opera- 
tional point of view is to have the closest range to the target. 

It is recommended that the various candidate modifications of 
the existing Ku-band radar to accommodate the specifications at short 
range be examined carefully with particular emphasis on those suggested 
in [9]. It is apparent from this discussion that the modifications will 
involve all aspects of the system. 


(a) One Horn - Rotated Mechanically 


Position of Horn for 



Position of Horn for 
Short Range Operation 




Figure 15. Candidate Short Mode Antenna Configurations 
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8.0 RANGING WAVEFORM CANDIDATES FOR SHORT RANGE OPERATION 

. The latest trend tov'/ards defining the short range operation of 
the KU'-band system is to consider two intervals: 3000 to 300 feet and 
300 to 30 feet. Because the target return power at these ranges is more 
than adequate, signal-to-noise ratio is not a limiting factor. On the 
other hand, accuracy of the range and range rate (velocity) measurements 
becomes of paramount importance. In view of these requirements, one 
should consider possibilities of using waveforms other than conventional 
pulse for short range operation. 

In Appendix E, several alternate (to pulse) candidate waveforms 
for short range (R<3000 ft) radar are introduced. The basic require- 
ment in choosing any' one of these waveforms is the high range resolution 
accuracy, preferably combined with high range rate estimation accuracy. 

The latter can be achieved by the use of a bank of doppler filters or 
via the a-g estimator at these short ranges. 

The ranging waveforms under consideration are: (1) linear FM 

signal, (2) short pulse signal, (3) tone ranging signal, and (4) pulse 
coded signal. The resolution capabilities, as well as the ambiguity 
functions, of these waveforms are presented and compared in Appendix E. 

However, in order to choose the best ranging v/aveform for close 
range tracking, further rigorous analysis of all the candidates is 
required. This analysis should take into consideration all resulting 
implementation losses, as well as the added complexity (if any) of the 
transmitter and the receiver. These factors should be weighed against 
the performance capability of the existing pulse system. 
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9.0 MAXIMUM LIKELIHOOD ESTIMATOR DESCRIPTION 

In the radar tracking mode, a logarithmic error discriminant is 
used [10] to obtain error signals for the various parameters of interest, 
namely, range, range rate, azimuth angle, and elevation angle. All of 
these parameters are mathematical duals of each other and a single model 
can apply to all of them. In each case, the information is imbedded in 
a pair of voltages ]J and _V, whose difference is approximately proportional 
to the target parameter under consideration. 

Both IJ and )I_ can be assumed to be narrov/band processes and a maxi- 
mum likelihood estimate of any desired parameter can be deduced by taking 
the logarithm of the powers of the two vectors; thus, 

2 = log (uVv\) . 

In Appendix F, the maximum likelihood (M.L.) estimate of the ratio of the 
power in two narrowband processes is considered. These processes are 
assumed to be Gaussian and, in general, will be correlated. The in-phase 
and quadrature phase components of each process are sampled N times, where 
the time between samples is assumed to be large enough so that the samples 
are statistically independent. Numerical computation of the standard 
deviation formula (31) in Appendix F is shown in Figure 16 for N= 2,5,10. 
In the figure, p is the normalized correlation coefficient between the 
two processes. The results agree closely with those derived from an 
empirical formula in [10]. 
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p = Normalized Correlation 

Figure 16. standard Deviation of the Random Variable z as a 

Function of Correlation 
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10.0 CONCLUSIONS AND RECOMMENDATIONS 

10.1 Power Budgets 

It has been shown that, with the present set of system parameters 
(as of April 1977), v/hen the radar is searching for a 1 m passive target 
in the "Autotrack" mode at 12 nmi, the design margin is -0.5 dB. This 
is a reduction of 2.5 dB in design margin from the previous set of system 
parameters. Host of the loss is due to the significant increase of range 
gate straddling loss for the new set of system parameters. This may not 
be a serious situation in that the power budgets were performed at 12 nmi, 
and all radar tracking loops and estimators must be operating by the 
time range is reduced to 10 nmi. This is a time period corresponding 
to 2 nmi to initiate track- For a coherent radar, vnth digital signal 
processing, it is expected that this can be carried out in much less 
time. A gain of 2.3 dB in design margin is obtained by carrying out the 
power budget at 10.5 nmi. This is due to range reduction alone. There 
is also a 1+ dB gain in range gate straddling loss at 10.5 nmi. 

Therefore, even though the design margin is still not very large, 
it is at least a couple of dB on the positive side. 

It is recommended that further power budgets be carried out for 
the GPC-ACQ mode at 10-12 nmi to correctly assess the effect of range 
designate errors on search and track performance. Also, at short ranges 
(R<2 nmi), similar pov/er budgets are recommended in both modes with 
particular emphasis on the effect of errors in range designate. 

10*2 Angle Track 

A thorough analysis of the Ku-band radar which is angle tracking 
a passive target is presented in Appendix C. The analysis is applicable 
to all modes of operation which require angle tracking, namely, GPC-ACQ, 
GPC-DES, and Autotrack. At all ranges of interest (R< 10 nmi) for passive 
targets, the angle tracking system has a design margin of at least 70 dB. 
The crosscoupling errors into the angle tracking loops from the range 
and doppler tracking loops v/ere described but not evaluated. With the 
present design margin and the fact that range and doppler crosstalk track- 
ing errors may not be significant, it is expected that the design margin 
for the angle tracking loop will remain adequate. 
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10.3 Range and Doppler Tracking 

It is recommended that the range tracking loop and doppler tracking 
loop be analytically developed in a manner similar to that of the angle 
tracking loop presented In this report. It is anticipated that both the 
range and doppler tracking loops v/111 not have a design margin as satis- 
factory as the angle tracking loop. 

In addition, it is recommended that the a- 6 tracking loop in the 
range tracking loop and the angle rate estimator be similarly developed. 
This would then clarify the performance capability of all aspects of the 
Ku-band radar when tracking a passive target in all modes in which tracking 
is required. 

10.4 Stationkeeping 

Preliminary considerations concerning the stationkeeping require- 
ments of the Ku-band radar are presented in Section 7,0, Some suggestions 
are given as to how the presently configured system could be altered to 
meet the stationkeeping requirements, both in the presence and absence 
of clutter. The geometry of the Skylab reboost and the military satel- 
lites are such that it will be necessary to provide stationkeeping in 
the presence of earth clutter. 

There are several areas that are deserving of significant effort. 

In the Short-Short mode (100-300 feet) and the Short-Short-Short mode 
(35-100 feet), the present antenna is operating in the near field. This 
has significant implications on the received signal and in monopulse 
tracking. One potential solution is the use of the feed-only part of 
the Ku-band antenna when in these modes. Preliminary studies show that 
this candidate configuration has promise. It is recommended that this 
potential modification, as well as other candidates, be carefully 
evaluated. 

The other major area is the radar signal processing in the Short 
mode, Short-Short mode, and Short-Short-Short mode. One particular 
implementation is presently being evaluated, which totally eliminates 
clutter for all ranges to the earth up to 500 nmi. Other implementations 
have been proposed for these modes and are also under investigation. 
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10.5 


Ground Clutter 


Analysis of radar operation in the presence of main lobe ground 

clutter indicates that both acquisition and tracking requirements can be 

met for Skylab reboost missions with the present system configuration. 

This applies even to the case of nadir return. Meeting the general spe- 

2 

cifications, that is, of acquiring and tracking aim passive target at 
12 nmi may not be possible unless either the radar modulation formats 
are modified from their present configuration or angles of incidence 
close to vertical are avoided. The investigation of methods for altering 
the modulation formats without impacting considerably the existing design 
is recommended. The goal of such an investigation will be to provide 
the required. radar mode capability regardless of incidence angle and of 
the magnitude of the ground return. 

The areas requiring further analysis, as outlined above, will be 
included with the remaining tasks of the program. 
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APPENDIX A 


POWER BUDGETS FOR THE KU-BAND RADAR FOR THE SPACE SHUTTLE 
ORBITER VEHICLE OPERATING IN THE AUTOTRACK MODE 

by 

Charles L. Weber 
SUMMARY 

Power budgets are developed for the Ku-band radar in the "Auto- 
track" search mode for a passive target at 12 nmi. Two cases are con- 
sidered: (1) the improvement due to the use of frequency diversity is 

negligible (slow fading), and (2) the improvement due to the use of 
frequency diversity is ideal (fast fading). 

For the slow fading case, the design margin is -6.0 dB, while 
for the fast fading case, the design margin is +2.0 dB. The latter is 
the more realistic case. 

An alternate set of power budgets is also presented v/hich incor- 
porates recent changes in system parameters. Although these changes 
were instituted to decrease implementation complexity, there is a degra- 
dation in radar search performance in the passive mode. For the slow 
fading case, the design margin is lowered to -8.5 dB and in the fast 
fading case, the design margin is reduced to -0.54 dB. 

Some of the cases presented are rather conservative. It is there- 
fore anticipated that performance in some cases will be improved and 
that the design margins v/ill be improved over those presented. 

There are various alternatives to improve the design margin to a 
more satisfactory value, namely, 6 to 10 dB for passive radar. 

In all cases considered, the gain of the main antenna was assumed 
to be G = 38.5 dB, which is the most recently knov/n value [9] 

Furthermore, all of the results in this report assume that tv/o scans 
of the column to be illuminated can be accomplished in 1 minute. At 
12 nmi, the two-sided cone angle that must be illuminated is 40 degrees. 
It has been shown [1,2] that, at these values of cone angle and range, 
two scans can be completed in 1 minute with a spatial dwell time of 
36 msec. This is the spatial dwell time used in this report. 
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POWER BUDGETS FOR THE KU-BAHD RADAR FOR THE 
SPACE SHUTTLE ORBITER VEHICLE 

by 

Charles L. Weber 


1,0 INTRODUCTION 

Sample power budgets are presented for the Orbiter Ku-band radar 
described in [1], Most of the hardware losses predicted in [1] are used 
in the budget computations in this report. The passive automatic search- 
and-detect mode only is considered. This mode is the most critical from 
the required transmitted power point of view. Each of the parameters is 
considered individually before the composite pov/er budgets are assembled. 


2.0 THE RADAR EQUATION 


The radar equation that is appropriate for passive mode detection 
of a pulse doppler radar which is essentially using an optimal receiver [2-4] 
is given by 


R = = 

V Nn 


[(4.)^R'*][kT5][L] 


where R = ensemble averaged peak SNR required per pulse for a given 
P probability of detection, P^, and false alarm probability, 
Pfg. In the case of an optimal radar which employs fre- 
quency diversity, the time of the pulse corresponds to 
the total time during v/hich coherent addition takes place 
during one frequency dwell time. 

E = peak signal energy received during the time in which 
coherent addition takes place 

Nq = equivalent one-sided noise pov/er spectral density of the 
entire receiver system = kTs 

a = average target cross section 

X = radar RF wavelength 

G = peak antenna power gain 

Pp - peak transmit power = P^^g/d^ 

d^ = transmitter duty factor 

- average transmit power 
avg 

= total time of coherent addition at one RF = tM 
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T = transmitted pulse width 

H = number of transmitted pulses per RF frequency 
R = range 

k = Boltzmann constant = 1.38x10”^^ Joules/*K 
T_ = system noise temperature, ®K 

L = total system losses, 

3.0 SEARCH RADAR SYSTEM PARAMETERS 

Each of the above parameters is now considered In regard to the 
Ku-band radar for the Space Shuttle Orbiter Vehicle described In [1], 
The reasoning and/or source of the choice of each parameter value is 
Indicated. 


3.1 Required Peak Siqnal-to-Hoise Ratio Per Pulse 

R. 


= ^ 
■p N, 


’0 


This value is dependent on the assumptions made about the mode 
of detection. The assumptions to be used in these power budgets are 
among the follovnng cases. The overall probability of detection is 0.99 
and the false alarm probability is 10” . The 0.99 value is a specifica- 
tion requirement, and the a =10”^ value is approximately equivalent to 
one false alarm per hour, which is the specified value. The value of 
Rp is relatively insensitive to the choice of a. Hence, the choice is 
not critical. The cases shown in Table 1 all assume two scans, so that 
the probability of detection per scan is 0.9. The follov/ing definitions 
are used in Table 1: 

N = number of pulses, noncoherently added 
Pj = probability of detection per scan 

n' = false alarm number 
a = 0,693/n' . 

The table takes into account target fluctuation loss in both the slow 
and fast fading cases- 
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Table 1. Required Value of Signal-to-Noise Ratio to Meet the Indi- 
cated Performance Levels Under the Indicated Assumptions 



Requirements 

SNR Needed 

Case 

Description 

n 

p ** 

•^d 

n' 

Rp (dB) 

Ep/Ng (dB) 

1 

Sv/erling I 
(Slow Fading)* 

6 

0.9 

10® 

18.2 

15.2 

2 

Swerl i ng ' I I 
(Fast Fading)* 

6 

0.9 

10® 

11.7 

8.7 

3 

Swerling I 

1 

0.9 

10® 

24.2 

21.2 

4 

Nonfluctuating 

Target 

1 

0.9 

10® 

16.3 

13.3 

5 

Nonfluctuating 

Target 


0.9 

10® 

10.1 

7.1 

6 

Swerling II 


0.9 

10® 

12.6 

9.6 


The phraseology being adopted (somev/hat incorrectly) is that the 
Sv/erling I model corresponds to slow fading and Swerling II to fast 
fading. The interpretation is based on the following assumptions, namely, 
that frequency diversity is being employed and that the target is initi- 
ally modeled by Swerling I fluctuations. If the frequency diversity is 
not effective so that the Swerling I model still applies, then the term 
"slow fading" is used. If the frequency diversity is completely effective, 
and independent samples are being received upon changes in RF, then 
Swerling II applies, and it is called "fast fading." 

The cumulative probability of detection is specified at 0.99. 

With two scans of the column in 60 seconds, each scan then has a proba- 
bility of detection of P^ = 0.9. 


1 I I -J:. 
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3.2 Average Target Cross Section 

2 

a = 1 m , as per specification 
= 0 dB. 

3.3 Radar RF Wavelength 

The RF for the radar is 13:775 GHz. This corresponds to 
A = 0.0218 meters 
= -33.24 dB. 

3.4 Peak Antenna Power Gain 

The antenna diameter has been specified as d= 36" = 0.914 m. 

Based on the rule of thumb, 

0g (deg) - ( 180 /i-),d/x) , 

the 3 dB beamwidth is 0g = 1.32 deg. The 3 dB beamwidth quoted in [1] 
is 1.6 deg. Based on the rule of thumb [8], 

G = (e„ In deg) , 

(V 

the peak power gain of the antenna is 
G = 41 dB . 

The peak power gain in [1] is .38.9 dB, so that 
6^ = 77.8 dB . 

More recently, however, the peak power gain is given as 38.5 dB [9], so 
that 

G^ = 77.0 dB . 

The latter value is used in the power budgets in this renort. 

3 . 5 T. ^ansml tted Signal Characteristics 

Pp = peak transmitter pov/£r 

r 

Given all other parameters and assumptions in the radar equation, 
the required peak transmitter power will then be obtained. 


Nl .1 I j 1 1 . ? J I 
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For the automatic search mode of a passive target, the PRF is 
presently set at 2666, so that the pulse repetition interval (PRI) is 

PRI - Tp = 375 psec . 

Also, the pulse vndth is 

T = 74.7 psec , 

so that the transmitter duty factor is 

d^ = 0.20 

and the maximum unambiguous range is 30.4 nmi. The spatial dwell time , 

Tjj, based on the algorithm described in [1] is 

= 36 msec . 

This also agrees with that developed by Udalov [5] for the Ku-band radar 
search mode of a passive target. 

The timing of the RF signal during one spatial dwell time is as 
shown in Table 2, where the frequency dwell time is seen to be 6 msec 
over each of the five RFs. 

Table 2. One Spatial Dwell Time 


RF 

fr 

^2 

^3 

^4 

>^5 

*^6 

Total 

Number 
of RF 
Pulses 

16 

16 

16 

16 

16 

16 

96 

msec 

6 

6 

6 

6 

6 

6 

36 


There are 16 RF pulses during each RF frequency dwell time. 

These are coherently added in the digital implementation, since the 
bandwidth of each of the doppler filters is matched to that associated 
with 16 pulses at one RF. 

This is not a centerline radar, since there is no centerline 
filter, and since wideband sampling is performed at IF. This allows the 
coherent addition of all coherent lines in the receiver RF spectrum over 
the entire IF bandwidth. 
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The total time of coherent addition at one RF is 

= T M , 

where t = 74.7 sec 
M = 16 pulses . 

Therefore, 

= total time of coherent addition 
= 1195.2 ysec 
= -29.2 dB . 

3.6 Range 

The specification calls for meeting the above requirements for 

_ p 

a (a=m ) passive target at 12 nmi = 22,236 km. Hence, 

R = 22. 236 X 10^ meters 
R^ = 173.9 dB. 

3.7 Receiver Noise Parameters 

The equivalent one-sided power spectral density of the receiver 
noi se i s 

Nq = kTg watts/Hz, 

-23 

where k = Boltzmann constant = 1.38x10 Joules/°K 
= -228.6 dB 

Tg = overall system noise temperature 

The components in the overall system noise temperature can be 
computed from the diagram shown in Figure 1. 



Figure 1. Simplified Diagram of Losses for Ku-Band Radar Receiver 
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The receiver is assumed to have a noise figure of 5 dB. All component 
losses [1] between the antenna and the receiver add to 2.8 dB, 

When the first stage of a receiver is an attenuator, as is the 
case here, it can be shown [6] that the systein noise. temperature, com- 
puted at the antenna, is 

h = (L-DTg + LT^. 

where = antenna temperature 

= additional antenna temperature due to clutter 

L = total attenuation from all lossy components between the 
antenna and receiver input 

Tg = receiver noise temperature = l)Tg 

= receiver noise figure 

Tg = standard room temperature = 290“ K . 

The design in [1] assumes the receiver noise figure is F^=5.0 dB or, 
equivalently, T_ = 627“K, and the attenuation from all lossy components 
between the antenna and receiver front end is L = 2.8 dB. The antenna 
temperature is assumed to be T^=5.5°K and the additional antenna temper- 
ature due to earth clutter when looking horizontally is T =36.5“K. Then 

= ^A-^'^c-^^L-DTgH-LT, 

= 5.5 + 36.5 + 263 + 1195 
= 1500“K 
= 31.76 dB . 

This is increased somev/hat as the look angle increases. When looking 
directly tov/ards the earth, the clutter temperature increases to 263. 5°K, 
and Tg increases to 

Tg = 1727“K = 32.37 dB. 

3.8 Losses 

In this section, the pertinent losses are listed, along with the 
references where the values v/ere derived or where the losses were specified. 
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3.8.1 Transmit Losses 

The RFP and the response [1] have specified the transmit losses at 
3.5 dB. At this point, no attempt is being made to dispute this value, 
so it is accepted and used in the computations. 

3.8.2 Scan Alignment Loss and Lateral Scan Loss 

The specification calls for a total beam shape loss of 2 dB. In 
[3], it is shov/n that the scan alignment loss is 1.03 dB each v^ay and the 
lateral scan loss is 0.58 dB each way, when the antenna movement from one 
circle to the next in the spiral scan is 0.6 6g. The spiral scan described 
in [13 and [5] calls for an antenna movement of 0.7 9g. The difference 
in scan alignment loss and lateral scan loss between adjacent antenna 
paths of 0.6 0g and 0.7 0g is negligible. Thus, the two-v/ay loss due to 
beam shape which includes scan alignment loss and lateral scan loss is 
set at 2(1.03+ 0.58) = 3.2 dB, 

3.8.3 Threshold Loss (Constant FAP) 

The specification calls for a threshold loss of 1 dB due to main- 
taining a constant false-alarm probability by monitoring the system noise 
level. This value is used in the computation, realizing that it should 
be possible to reduce this value to a few tenths of a dB. 

3.8.4 Processor Loss 

It is estimated that the algorithms in the radar digital signal 
processor will result in approximately 1.7 dB loss [1]. This value is 
also used in the power budget, while anticipating that it may be possible 
to reduce it. 

3.8.5 Range Gate Straddling 

The present implementation calls for four range gates when the 
initial range to the target is not designated by GPC [1], as shov/n in 
Figure 2. Each range gate is of the same width as the transmitted pulse. 
The corresponding values of range are shown at the beginning of each 
range gate which corresponds to zero straddling loss. For example, at 
12 nmi, the received pulse will show up exactly in Range Bin 2, As the 
range decreases, straddling loss increases, but not as fast as the R^ fac- 
tor increases signal povier. The net effect is that, at 12 nmi, there is 
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Tp= 375 usec 

Transmit Pulse 



74.9 vsec 





75 vsec 

75 psec 


Range 


Bin 3 


Bin 1 

Bin 2 

6 nmi 1 

2 

18 


75 ysec 
Bin 4 


Figure 2. Range Gate Diagram 



Figure 3. Description of Bank of Doppler Filters 
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no range gate straddling loss and, at all shorter ranges, the SNR increases. 
As a result, the range gate straddling loss is set to zero. Straddling 
loss and eclipsing loss are analyzed in detail in [7]. 

3.8.6 Doppler Filter Loss 

Following each of the four range bins is a bank of 16 doppler 
filters, as shown in the block diagram of Figure 3. These filters are 
uniformly spaced over the frequency range ^PRF^* passive 

search and detection at 12 nmi, the PRF is 2666 Hz, so that the frequency 
difference between the center of adjancent filters is 166 Hz, as shown 
in Figure 3. If the 3 dB (one-sided) is set at FpRp/16, then at the points 
where the adjacent filter gains cross (at ±83 Hz about the center of the 
filter), the gain is down approximately 1 dB for a first order Butterworth 
filter, and down approximately 0.7 dB for a [{sin x)/x]^ filter. This loss 
is classified as mismatch due to doppler. In our pessimistic power budget, 
we use -1 dB, and, in the more realistic budget, -0.7 dB is used. 

There may be additional loss due to the lack of a proper match 
between the pulse bandwidth, the PRF, and the bandwidth of the FFT. 

3.8.7 Pre-Sum Mismatch 


There is an additional loss due to doppler in the pre-sum. The 
pre-sum adds all samples over the width of one range bin from a trans- 
mitted pulse. There is no doppler compensation in the pre-sum. If there 
is no doppler, then this coherent addition is ideal and there is no loss. 
As the doppler frequency increases, the pre-sum loss increases. This 
loss is in both the I and Q channels and can be described in terms of 
the following integral and is shown in Figure 4. 


^PRE-SUM 



jZirfot 

e 


dt 


2 


where Tp is the transmitted pulse width and f^ is the doppler frequency. 
This is a very close approximation to the actual pre-sum loss, since the 
sampling rate is high. The pre-sum loss becomes 


•-PRE-SUH 


I sin (irfoTp) 
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which is plotted in Figure 5. For the radar in [1], Tp = 75 usee. The 
range rate varies over -75 fps to +148 fps. At f = 13.775 GHz, each 
fps is equal to approximately 28 Hz, so that the doppler frequency is 
v/ithin the range 

F[j s [-2100, +4144 Hz] . 

Ambiguity clearly exists, since the PRF= 2666 Hz, This ambiguity 
is being resolved by filtering and differentiating the range measurements. 
The maximum pre-sum loss therefore occurs at Fp= 2666 Hz. Then, 

TpFp = (75 usec)(2666 Hz) = 0.2. 

This is shown as the dot in Figure 5, and the resulting maximum pre-sum 
loss is 0.6 dB. 

4.0 POWER BUDGETS FOR PASSIVE SEARCH 

In Table 3, the results of the previous sections are summarized 
for the slow and fast fading cases. The assumptions that the target 
range is not designated and that detection must be accomplished at 
R=12 nmi are maintained. 

Budget #1, the slow fading case, assumes frequency diversity is 
present but totally ineffective and that independent samples are not 
obtained v/hen the RF frequency is switched. As a result, we have a 
Swerling I target. The postdetection integration is still noncoherent 
from frequency to frequency. 

With all of the assumptions listed in Table 3, the required 
average power in Budget #1 is 39.7 watts, which is higher than the design 
value of 10 watts [1] by 6.0 dB. The corresponding peak power is 198.6 
v/atts when the duty factor is d^=0.2. 

With Budget #2, the fast fading derived from the use of frequency 
diversity is assumed to be ideal. This is somewhat optimistic since one 
of the RF frequencies is used twice over the spatial dwell time. An 
underestimate of the power required would be to set the number of RF 
frequencies to N=5, and the spatial dwell time to 30 msec. 

For Budget #2, the required average power is 6.3 watts, and the 
required peak power is 31.5 watts. Both of these values are below the 
average TWT power required for the Ku-band communication system. Based 




1 


1 

r 


1 I ! 



[ 


15 


Table 3. Power Budget for Passive Searcli at 12 nmi Uitliout Range Designated 




Sect on 


Hudget iJl 

Oudget l!2 



Nuniher 


Case 

Case H2 



in This 


(Slow fade) 

{Fast Fade) 

Item 

Term 

Report 

Value 

(da) 

(dB) 

1 


3.1 

Required Peak S!1R 

18.2 

11.7 

2 

g2 

3.4 

G^3fl.5d0 

-77.0 

-77.0 

3 


3.3 

X- 0.0218 m 

33.24 

33.24 

4 

0 

3.2 

1 r/ 

0 

0 

5 

(4tt)^ 



33 

33 

6 

R-» 

3.6 

R “ 1 2 nmi = 22 . 2 kni 

173.9 

173.9 

7 

k 

3.7 

1.38X 10"^^ 

»228.6 

-228.6 

8 

Ts 

3,7 

ISOO^K 

31.76 

31.76 

9 

Coherent Addition 
Time 

3.5 

1,195 msec 

29.23 

29,23 

10 

T 

2.0 

2 

^3 

-3 


Losses 

3.8 

i 



11 

Transmit Losses 

3.8.1 

RFP 

3.5 

3.5 

12 

Beam Shape 

3.8.2 

Scan Alignment - 
Lateral Scan 

3.2 

3.2 

13 

Threshold 

3.8.3 

Constant FAP 

1.0 

0.5 

14 

Processor Loss 

3.8.4 ' 


1.7 

1.0 

15 

Range Gate 
Straddling Loss 

3.8.5 


0 

0 

16 

Doppler 

3.8.6 

Hi smatch 

1.0 

0.7 

17 

POI 


Postdetection Inte- 
gration BEF [1] 

1.25 

1.25 

18 

Pre-Sum Mismatch 

3.8.7 

Due to Doppler 

0.6 

0.6 

19 


3.5 

Required Peak Pov/er 

22.93 dB 

14.98 dB 


P 


{198.6 w) 

(31.5 w) 

20 


3.5 

Duty Factor = 0.2 

-7.0 

-7.0 . 

21 

p 


Required Average 

15.99 <JB 

7.99 dD 


avg 


Pov/er 

(39,7 w) 

(6.3 w) 

22 

Pp Design 


Reference [1] 

50 w 

SO v; 

23 

Pavg 


Reference [1] 

10 w 

10 W 

24 

Design Margin 



-6.0 dB 

+2.0 dB 


f 






1 


1 


1 
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on a peak power of 50 mtts and an average power of 10 v/atts, the design 
margin is 2.0 dB. Without the advantage of frequency diversity, as 
described by Budget #1, the design margin is -6.0 dB. 

The more realistic power budget is Budget where there is a 
positive design margin. It is also expected that some of the losses can 
be improved so that the design margin improves over that presented. 

5.0 ALTERNATE POWER BUDGETS FOR PASSIVE SEARCH 

In Table 4, alternate power budgets for passive search at 12 nmi 
are presented. The assumption that the target range is not designated 
is kept. In Table 4, the values which have altered from those in Table 3 
are designated with an asterisk. The new values are a result of changes 
suggested in [9]. The altered values are now considered. 

Item 3 . A new set of RF frequencies for the Ku-band radar have been 

specified, and are; 13.779 GHz 

13.831 GHz 
13.883 GHz 
13.935 GHz 
13.987 GHz . 

Using the center of these RFs, the value of the wavelength is 
X = 0.02161 m. 

Item 9 . The updated value of PRF for target ranges greater than 7.2 nmi 
is 2987 Hz and the updated value of pulse width is 66.4 ysec. Therefore, 

PRI = 335 ysec 
Maximum unambiguous range = 21.8 nmi 

Duty factor = 0.20 (unchanged) 

Sixteen pulses are still assumed to be transmitted at each RF 
during search, so that the frequency dwell time is 

Tt = (16){335 usee) = 5.36 msec 

and the coherent addition time is 

= (16)(66.4 ysec) = 1.062 msec. 

Number of Independent Observations of the Target . If the spatial 
dwell time of 36 msec is unchanged (which is assumed), then 6.72 RF 
frequencies illuminate the target during each spatial observation. 
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Table 4* Alternate Power Budget for Passive Search at 12 mni Without flange D'^signated 


Item 

Term 

Section 
Number 
in This 
Report 

Value 

Ouclget 113 
Case //i 
(Slow Fade) 
(dB) 

Budget M 
Case if2 
(Fast Fade) 
(dB) 

1 


3.1 

Required Pciik S!1R 

18.2 

11.7 

2 

g2 

3,4 

G=38.5 dR 

-77.0 

-77.0 

3 


3.3* 

0,02161 m 

33,31 

33.31 

4 

a 

1 3.2 

; ip,2 

0 

0 

5 

(411)3 



33 

33 

G 


3.6 

R= 12 nnti = 22.2 km 

173.9 

173.9 

7 

1: 

3.7 

1.38x10'^^ 

-228.6 

-228.6 

8 

Ts 

3.7 

1500“K 

31.76 

31.76 

9 

Coherent Addition 
Time 

3,5* 

1,062 msec 

29.74 

29.74 

10 


2.0 

2 

-3 

-3 


Losses 

3,8 




n 

Transmit losses 

3.8.1 

RFP 

3.5 

3.5 

12 

Beam Shape 

3.8.2 

Scan Alignment- 
Lateral Scan 

" 3,2 

3,2 

13 

Threshold 

3.8.3 

Constant FAP 

1,0 

0,5 

14 

Processor Loss 

3.8.4 


1,7 

1,0 

15 

Range Gate 
Straddling Loss 

1 

3.8.5* 


2,0 

2.0 

16 

1 

Doppler 

3.8.6 

Mismatch 

1.0 

0.7 

17 

PDI 


Postdetcction Inte- 
gration Ref. [1] 

1,25 

j 

1.25 

18 

Pre-Sum Mismatch 

3.8.7* 

■ ^ 

Due to Doppler I 

0.57 

0.57 

19 

'’p 

3.5 

Required Peak Pov/er 

25.53 dQ 
(357.3 w) 

17.53 dB 
(56.6 w) 

20 

‘•t 

3.5 

.Duty Factor- 0.2 

-7.0 

-7.0 

21 

P 

avg 


Required Average 
Power 

18.54 dB 
(71.5 tv) 

10.54 dB 
(11.3 w) 

22 

Pp Design 


Reference [1] 

50 w 

50 w 

23 

Pflvg 


Reference tU 

low 

10 w 

24 

Design Margin 



-8.5 dB 

-0.54 dB 
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Since there are 5 RFs in the frequency diversity algorithm, there are 1.72 
repeats during each spatial observation. Under the realistic assumption . 
that, on the average, the repeated frequencies are only 50% as effective 
as nev/ frequencies would be, the value of N = 6 independent observations 
of the target is used in the budgets of Table 4. 


Item 15 . With the new set of values proposed in [9], there is a range 
gate straddling loss. At 12 nmi, the pulse is received as shown in 
Figure 6, where four range gates are used in the search mode and range 
is not designated. The loss of received (coherently added) energy in 
Range Gate 2 is 


/5.45 - n 
^ 5.45 


2 

•) 


-2 dB . 


The square results from the fact that the received signal is being 
coherently {predetection) integrated. The straddling loss would be 
linear in the noncoherent case. 


12 nmi 
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\ XMTR i 


5.45 


10.9 


16.35 
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Time 
• — ^ 
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Figure 6. Range Gate Straddling Loss at 12 nmi 

Item 18 . There is an almost negligible change in the pre-sum mismatch 
due to doppler. From Figure 4, when Tp = 66.4 usec and f^ = 2987 pps, the 
pre-sum loss is 0.57 dB, 

All remaining entries in the budgets in Table 4 are carried forward 
from Table 3. For the passive search mode, it is noted that the new param- 
eters deteriorate performance with respect to the former values in Table 3. 
When frequency diversity is totally ineffective on a Swerling I target 
(Budget #3), the design margin drops from -6.0 dB to -8.5 dB. When the 
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frequency diversity is assumed to operate ideally (Budget M) and the 
various other loss improvements are used as described previously, then 
the design margin drops from +2.0 dB to -0.54 dB. 

By allovring 2 nmi (12 to 10) to initiate all track modes of oper- 
ation and get all RMS errors within specification, the detection must 
meet specification at 12 nmi. For coherent radar, it is most likely that 
this time is substantially less. For example, if these power budgets 
v/ere computed at 11 nmi instead of 12 nmi, +1.5 dB of design margin is 
gained, and if they were performed at 10.2 nmi, +2.8 dB is gained in the 
design margin. 
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1.0 INTRODUCTION 

The objective of the Ku-band radar system on the Shuttle Orbiter 
is to detect the presence of a targ'et and then obtain continuous accurate 
estimates of the various target parameters, namely, its range, range rate 
(velocity), azimuth, elevation angle, and angle rates. The target is 
usually detected when the radar is in the search mode , while the accurate 
estimates are obtained in the track mode . It is the objective of this 
appendix to clarify the overall radar block diagram and to explain its 
operation in both the search and track modes. 

The radar employs coherent processing over each RF frequency 
transmitted, and noncoherent postdetection integration (PDI). Following 
the intermediate amplifiers, in-phase and quadrature phase signals are 
sampled and digitized (A/D converted), providing target amplitude and 
phase information. For longer range search and acquisition (R> 3000 ft), 
16 doppler filters are formed, covering the doppler interval defined by 
the repetition rate. This is an ambiguous doppler interval; this ambi- 
guity is removed via differentiation of the range estimates. 

Following the filters, the target magnitude is determined, summed 
over the number of RF frequencies (noncoherent postdetection integration) 
and compared to a threshold. The threshold is set from noise measure- 
ments, so as to maintain a constant false alarm probability (CFAP). This 
can be implemented by any of several methods. 

In the search mode, only the sum channel is processed. When 
the threshold is exceeded and target detection is declared, the auxiliary 
antenna signal is processed and its magnitude is compared with that of 
the main antenna to eliminate sidelobe detected targets. The auxiliary 
antenna has a peak gain which is approximately 20 dB less than that 
of the main antenna. Also, the sidelobes of the main antenna are 
approximately 20 dB down from the main lobe. This provides roughly 
a 20 dB main lobe/guard antenna ratio to detect and eliminate sidelobe 
targets . 


f 
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In the modes for which angle tracking is required, namely, GPC 
acquisition and autotrack, the sum plus angle error channels are processed. 
The azimuth and elevation error signals are time-multiplexed, thereby 
eliminating the need for a second and third matched processing channel. 

In Section 2.0, a general block diagram is drawn and explained 
in more detail; the functions of the individual blocks are described in 
the subsequent sections. 

2.0 RADAR BLOCK DIAGRAM 

The Shuttle Ku-band radar block diagram is shown in Figure 1. The 
antenna assembly consists of a main antenna and an auxiliary antenna, 
followed by the deployed mechanism assembly (DMA) and the deployed elec- 
trical assembly (DEA), where various switching and filtering operations 
take place, as explained in Section 3.0. The antenna assembly is followed 
by several IF mixing stages to bring the carrier frequency down to a 
stable IF of fQ= 78.143 MHz. The output of the I-Q stage that follows is 
a pair of baseband signals which are in phase and in quadrature. with the 
input signal. The operation of the previous two stages is discussed in 
Section 4.0. The available information is converted by an A/D converter 
to digital data which samples both channels at a rate of f^= 430 kHz, 
which is slightly larger than the Nyquist rate of the baseband signals 
at the output of the I-Q stage. (The LPF bandwidth in both the I and Q 
channels is equal to 237 kHz.) The output of the A/D converter is a set 
of complex numbers which undergo various stages of digital processing 
before estimates of the target parameters (range, velocity, azimuth angle, 
elevation angle) are obtained. The first processing stage is a presumming 
stage, follov;ed by a set of range gates and doppler filters consecutively. 
The arrangement of the gates is varied for the two modes of operation. 

In the search mode, where range designation is not present, there are 
four range gates per pulse, while in the track mode, there are only two 
nonoverlapping gates: an early gate and a late gate. 

For long ranges, the magnitudes of the outputs of the doppler 
filters are summed over the RF frequencies and compared to a threshold 
to detect the existence of a target. For short ranges (R< 3000 ft), the 
detection takes place without using any doppler filtering. 

Once a target is detected, the radar is switched to a track mode 
and the digital signals undergo postdetection integration, followed by 
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Figure 1. Radar Block Diagram 
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an error estimation unit which is initiated by a crude estimate obtained 
by the search mode operation. The final estimates of the various param- 
eters are; 

(1) Fed as inputs to the servo motors for angle track. 

(2) Used to adjust the position of the early and late gates for 
range track. 

(3) Displayed for the astronauts for visual reading (all vari- 
ables being tracked). 

(4) Utilized to update the state of the error estimators to 
obtain new estimates. 

In the following sections, a detailed discussion of the various 
stages in the block diagram is presented, along v/ith an analysis of the 
signal timing in both the search mode and the track mode. 

3.0 ANTENNA ASSEMBLY 

The antenna assembly is illustrated in Figure 2. The main antenna 
has a gain G=38.5 dB. There is also an auxiliary antenna whose gain is 
G=20 dB, which is used to guard against target effect on the sidelobes 
of the main antenna. Three signals emerge from the antenna, namely, the 
sum signal (s) which carries the received signal, and two angle signals 
(aAZ and aEL) which supply the azimuth and elevation angle information. 

The E signal is passed through a bandpass filter (BPF), whose center fre- 
quency can be adjusted, by a centralized frequency synthesizer, to be one 
of the five frequencies (13.779, 13.831, 13.883, 13.935, 13.987 GHz). The 
angle signals aAZ and aEL are passed through a selecting switch v/hich 
selects one of the signals at a time at a rate of 93 Hz. The output of 
the switch is then passed through a phase encoder whose rate is tvnce 
that of the selecting switch. A typical sequence of signals out of the 
encoder is shown in Figure 3. The encoded angle signal is then passed 
through a matching BPF to that of the E signal. The two signals are then 
passed through the T/R switch to the receiver IF mixing stage. 

When angle track is operating, there are four antenna error dv;ell 
periods at each RF, during which the azimuth and elevation error signals 
are processed. During each error dwell period, 16 doppler filter outputs 
(via DFT or FFT) are computed. 


Main Antenna 
G=38.g dB 



Phase Switch 
0 - 180 ° 
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Figure 2. Antenna Assembly 
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Figure 3. TDMA and Phase Encoder Output of Angle Error 
Information, During Angle Track 

4.0 IF MIXING AND I-Q STAGES 

The antenna assembly is followed by tv/o stages of intermediate 
amplifiers and IF filters, where the frequency of the signal is beaten 
down to a stable 78,143 MHz, as shown in Figure 4. This, in turn, is 
followed by an in-phase and quadrature detectors (I-Q), where the signal 
is converted to two baseband signals. The two signals are converted to 
digital data using an A/D converter v/hose sampling frequency is slightly 
higher than twice the bandwidth of the signal to reduce the aliasing 
error to less than 0.5 dB. 

The A/D converter supplies 32 samples/pulse from each channel of 
the I-Q detector for each range gate period. These samples are summed 
at the PRESUM stage to give a complex number that represents every 
received pulse.* 

5.0 SIGNAL FORMAT 

Before explaining the different modes of operation, it, is important 
to discuss the signals that are being used. The signals are RF pulses 
with a duration (t) and a pulse repetition frequency (PRF) that are func- 
tions of the designated range. Frequency diversity, using five RF fre- 
quencies, is used in both the search and track modes. Table 1 lists 
these frequencies which are cycled in a given pattern which is convenient 
for the manual operation mode [1]. The various signal timings are sum- 
marized in Table 2. 

The 32 samples/pulse number applies to the search mode only, 
where the pulse width is 66.3 vsec. 
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Table 1. Frequency Synthesizer Frequencies for the 
First IF Mixing L.O. (in GHz) 




^3 

^4 


13.128 

13.180 

13.232 

13.284 

13.336 


Table 2. Dwell Times for Various Designated Ranges 
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Spatial 
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8-22 

66.3 

2987 

5.3 

31.8 

4-8 

33.2 

7177 

2.2 

13.4 

2-4 

16.6 

7177 

2.2 

13.4 

1-2 

8.3 

7177 

2.2 

13.4 

0.5-1 

4.1 

7177 

2.2 

13.4 
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6.0 THE SEARCH MODE 

In explaining the search mode operation, long range signal format 
is chosen for demonstration. All other ranges follow identical processing 
except for the short range (<0.5 nmi) which does not utilize the doppler 
filters. 

Sixteen pulses are transmitted at each RF frequency. The duration 
of each pulse is 66.3 nsec, while the time betv/een consecutive pulses is 
335 nsec, as shown in Figure 5. Four range gates (R-j, Rg, R3, R^) are 
used to cover the designated range. It is important to note that no 
doppler correction is used at the presummer in the search mode. For each 
RF frequency, 64 complex numbers are stored in a memory designated as 
M{1), as shown in Figure 6, The output of each range gate at every RF 
frequency is passed through a bank of 16 doppler filters — implemented as 
discrete Fourier transforms (DFT) or fast Fourier transforms (FFT)— and 
the magnitudes of the outputs of these filters are calculated and compared 
to precalculated thresholds. The target is detected v/hen the magnitude 
of one or two outputs of adjacent doppler filters pertaining to one or 
two range gates is exceeded. An initial estimate of the range and the 
range rate of the target are calculated from the knowledge of the doppler 
filters and range gates whose outputs have exceeded the threshold. These 
initial estimates are fed into the error estimators as initial conditions 
to start the successive estimation process in the track mode. The thresholds 
are calculated to produce an error alarm rate of one false alarm per hour. 

7. THE TRACK MODE 

7.1 Signal Format 

The signal in the track mode consists of 64 pulses for each RF 
frequency, instead of the 16 pulses in the search mode, as shown in 
Figure 7 [1,2]. For long range tracking, the frequency dwell time is 
equal to 21.44 msec. Each frequency dwell time is divided into four time 
slots by the switching combination of SWl and the phase encoder switch 
in the antenna assembly. The signals in the succeeding time slots are 
E + aAZ, E-AAZ, e + aEL, and E-AEL, respectively. The time slots are 
arranged so that the first two are used to measure the azimuth angle, 
v/hile the last two are used for the elevation angle. Each one of the 




0 nmi 



















Figure 7. Track Mode Signal Format 
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previous time slots consists of 16 pulses. The pulses are passed through 
an early gate and a late gate, which are located around a predicted esti- 
mate of the range [Rp(n)] from the processor timing unit. The width of 
each gate is equal to the pulse width. The output of each range gate is 
a complex number z{I, J, L,K), which is a function of four parameters {I, 
0,L, K): 


I denotes the RF frequency being transmitted; 1 = 1, 2,. ..,5. 

J denotes the time slot in each frequency such that 
J=1 corresponds to S + MZ 
J = 2 corresponds to E - aAZ 
J = 3 corresponds to s + AEL 
J = 4 corresponds to e-aEL. 

K denotes a particular pulse in each time slot; K=l,...,16. 

L denotes the range gate under consideration: 

L = “1 corresponds to Early Gate 
L = +l corresponds to Late Gate. 

Figure 7 illustrates the signal format in the track mode. 


3ler Filterinc 


The outputs of the range gates are stored in a memory bank desig- 
nated for illustration as M{1). For all combinations of I, J, and L, the 
series of 16 pulses (K=l,...,16) are passed through a bank of 16 doppler 
filters, as shown in Figure 8 [1,5], The doppler filters are uniformly 
spread over f^,, fc'‘"^PRF* overlapping of adjacent filters results 
in doppler filter loss, which has to be taken into consideration. The 
outputs of the doppler filters are designated as o)(l, J, L, M) , where M 
denotes a particular filter {M= 1 ,. . . ,16). I, J, and L represent the RF 
frequency, the time slot, and the range gate, respectively. Since doppler 
filtering is accomplished by digital Fourier transforms {DP"i/, m(I, 0,L,M) 
is given by 


w(l, J, L,M) 


16 

I z(I,J,L,K) e 
K=1 


(j2irKM)/16 


( 1 ) 


Implementing doppler filters using FFT might result in losses due to the 
lack of proper match between the pulse bandwidth, the PRF, and the band- 
width of the FFT. 
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Figure 8. Pictorial Illustration of Memory Bank M(l) 
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The outputs of the doppler filters are stored in a memory bank 
designated as M(2), as shown in Figure 9. Two possibilities might occur: 

(1) Two adjacent filter outputs exceed a prespecified threshold. 
All their values for various RF frequencies (I), time slots (J), and 
range gates (L) are then postintegrated. 

(2) Only one filter's output exceeds the threshold; then the out- 
puts of the filter along with one of the adjacent filters are postinte- 
grated. 

For illustration, it is supposed that the outputs of the filters 
M=Mi,M.j+l are processed. 


7.3 Postintegration 

The outputs of the two doppler filters that have exceeded the 
threshold consist of 40 memory cells from the memory H{2), as shov/n in 
Figure 9. The 40 cells correspond to five RF frequencies (1 = 1,5), two 
range gates (L = ±l), and four time slots (0 = 1,4). Each one of these 
cells carries a complex number w(I,J, L,M) , whose magnitude is found 
by multiplying it with its conjugate. The squares of these magnitudes 
are added in different manners to form the postdetection integration 
(PDI) outputs. To be able to demonstrate the logical procedure for add- 
ing these magnitudes, the memory cells under consideration are redrawn 
for a given RF frequency F^- in Figure 10. 
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Figure 10. Memory Cells Storing 6 j(I,J,L,H) for a Given RF 
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Eight sums are formed at the output of the postdetection inte- 
grator — tv/o corresponding to each of the four parameters of interest, 
namely, range, range gate, azimuth angle, and elevation angle. 

7,3.1 Range PDI Outputs 

The tv/o range outputs are denoted Sj^(L = -l) and S|^(L = +1). The 
R denotes the range and the L denotes the early or late gate. Thus, to 
obtain Sj^(L = -l), the square of all the complex numbers pertaining to 
the early gate are added over all the RF frequencies, Sj^{L = +l) is 
obtained in a similar fashion, using the late gate; hence, 

5 4 

= I I ( 2 ) 

i=l j=l ‘ ' 

5 4 

Sr(L = +1) = I I |4i,jM,MO|^+io)(i,j,l,M,+l)|2. (3) 

^ i=l j=T ' ' 


7.3.2 Range Rate (Doppler) PDI Outputs 


The doppler PDI outputs Sp(M=M-]) and Sp{M=M-j+l) are obtained 
by summing the square of the magnitudes of the numbers out of the doppler 
filters M-j and M-j+1, respectively, rather than the outputs of the early 
and late gates, as was done ror the range computations. Therefore, 


Sp(M=M^) 


i=l j-1 * ‘ 


(4) 


Sp(M=M^+l) = n‘^{UJ.-UM^+l)|^+>(i,a\l, 


(5) 


7.3.3 Angle PDI Outputs 

To obtain the two PDI outputs for the azimuth (elevation) angle 
information, the squares of the magnitudes of the first (last) two time 
slots corresponding to E + aAZ and z - M2 (z + AEL and z-aEL) are added 
over all RF frequencies to give and S^2^'^“2) [S^j_(J = 3) and 

SEi(d = 4)], respectively. The various angle PDI outputs are shown in 
equations (6) through (9). 
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(6) 

(7) 

( 8 ) 

(9) 


7.4 Logarithm Discriminant Formation 

The purpose of the logarithm discriminant is to obtain an estimate 
of the error (error signal) including sign of the various parameters under 
consideration. The idea is best explained by using the well known Taylor 
series expansion of ln[(l + x)/(l - x)]. 

In ^ = 2x+^ + ^+... ; ixl < 1 . (10) 


For small x , the following approximation is valid: 

In ^ - 2x; x« 1 . (11) 

Thus, if x represents the errors in the parameters under consideration, 
then forming the logarithm of (l+x)/(l -x), this error can be computed. 
The mathematical derivation of the variance of the logarithmic estimator 
is formed in [3], where the maximum likelihood estimator of the form 

I |U,l^ 

I IV, 1' 

is analyzed. U^. and represent a pair of voltages whose difference is 
proportional to the target parameter being measured. Since the target 
range, range rate, and orientation angles are mathematical duals of each 
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other, 1t is the objective of the following discussion to relate |U,. | 
and jV^. I , shown in (12), to the various outputs of the PDI. 

7.4.1 Range Discriminant 

In the formation of the range discriminant, a sum signal (S) and 
a difference signal (D) are formed such that 

S = Sj^(L=-l) + S^(L=+1) (13) 


D = Sj^(L=-l) - S^(L = +1) (14) 

where S^( ) is defined in equations (2) and (3). 

Since the output of the early gate is a measure of R- aR, where 
AR is an absolute error in the range, and the output of the late gate is 
a measure of R+aR, then S is a measure of the range R and D is a measure 
of the absolute range error. Then, replacing S by R and D by AR, it can 
be stated that 


Sr(L=-1) « R-AR 

(15) 

S|^(L = +1) « R + AR. 

(16) 


(A proportionality constant might be added; however, this coefficient 
will drop when taking the logarithms.) 

Letti ng 

20 « 

Sr{L = -1) 5 I lu^i^ 

^ i=l ^ 

20 5 

Sr(L=+ 1) e I |V.|S 

1=1 1 


(17) 

(18) 


it is easy to see that, using (15) through (18), 
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(19) 


where R^ is the required range error signal. The procedure is illustrated 
in Figure 11. 
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Figure 11. Range Discriminant 


7.4.2 Doppler Discriminant 

A similar argument can be stated for the doppler discriminant. 
In this case, from (4) and (5), 
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I |V,-r = S„(M=M,+1) , 
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from which 

20 2 

20 2 - Sb(H=H,) - ‘ ' 

i=l ^ 


( 20 ) 


7.4.3 Angle Discriminants 

As has already been noted, the data from the DFT representing the 
early and late range gates are summed to form an on- target gate, and then 
integrated over the transmitter frequency cycle. The sum and difference 
samples of angle error are then log-converted and subtracted to form the 
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normalized angle errors for the elevation and azimuth channels. The 
error signals are then D/A converted and sent to the servo filters and 
finally to the servo motors. 

In particular, for the azimuth angle, the estimation vectors U^. 
and V^. satisfy [from (6) and (7)]: 

20 5 

.1^ |Uit = 

20 9 

I “ s^2 • 

Similarly, for the elevation angle [from (8) and (9)]: 

20 9 
I lU.r = 

i=l ^ 

20 9 

I jv.r = = 

1=1 


from which the error signals become 


2a^AZ - 


2A AZ 


A 


_ 2 A EL 


( 21 ) 

( 22 ) 


The error signals in (19) through (22) are used to obtain the 
various required estimates pertaining to the parameters under consider- 
ation. 

Since the error signals are obtained successively as new readings 
are processed, they will be denoted by R (n), R^(n), A AZ(n), and A EL(n), 
where (n) denotes the nth error signal being processed. 

The angle error signals A^AZ and A^EL are fed into loop filters 
as shown in Figure 12 to obtain two sets of estimates: 

(1) Angle estimates aAZ and AEL, which are fed back to the 
servo motors. 

r*i A 

(2) Angle rate estimates AAZ and aEL, which are displayed for 
the astronauts for visual reading. 
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The range and range rate errov signals, however, are fed into the processor 
timing unit to obtain various required estimates, as discussed in the fol- 
lowing section. 

7.5 Processor Timing, Range and Range Rate Prediction 

The processor timing stage consists of tvro substages. The first 

is called the 0-3 tracker [4] and the second is designated as the ambiguous 

doppler resolver and corrector, as shown in Figure 13. The a-6 tracker 

is employed to give a smoothed estimate of the current range ft (n), a pre- 

dieted estimate of the next range reading R_(n+1), and a smoothed estimate 

P 

of the range rate Rg{n). This latter value is used along with the range 
rate error signal ft (n) as inputs to the ambiguous doppler resolver and 
corrector to produce a smoothed corrected doppler estimate R, (n) and a 

' ^ ” Sw 

predicted corrected doppler es^-imate R^^(n+1). The previous estimates 

pc 

are used for various purposes as summarized in Table 3. 


Table 3. Processor Timing Stage Outputs 


Estimate 

Rp(n+1) 

R,(n) 



Use 

Adjust early 

Display 

Display 

Adjust 


and late gates 



doppler 


location 



filters 


The one-step prediction and smoothing equations of the a-3 tracker, 
as well as a block diagram of a possible implementation of the tracker are 
shown in Figure 14. The equations v^ere derived in [3,4], along with the 
variances of these estimates. The initial conditions R(0) and R(0), which 
are necessary to start the a-3 tracker operation, are supplied by the 
search mode operation after the target is detected. 

The operation of the ambiguous doppler resolver and corrector is 
explained in Section 7.6. 

7.6 Doppler Ambiguity Resolution 

The doppler discriminant and resulting estimation will give a 
value between [0,PRF]. The maximum closing velocity is +148 fps, and 
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the maximum opening velocity is -75 fps. At 13.9 GHz, this corresponds 
to approximately 28 Hz/fps doppler. The maximum opening doppler is 
therefore -210P and maximum closing doppler is +4144 Hz, as shown 
in Figure 15. 



Figure 15. Resolution of Doppler ^biguity 

When the high PRF of 7177 is used, there is no doppler ambiguity 
from the doppler estimation from the a-3 tracker. For the low PRF at 
longer ranges, namely, 2987, there is a doppler ambiguity as indicated 
in Figure 15. The resolution is easily resolved from the more coarse 
estimation of range rate from the a-3 tracker by noting the following. 

(a) If the coarse ft^(n) > 0 (opening), then the fine estimation 
R-_(n) will be in Region 3', whereas the actual estimate should be in 
Region 3. In this case Rg(,(n) is adjusted by subtracting the PRF equal 
to 2987 Hz. 

(b) If the coarse Rg(n) is longer than 2987 Hz (high closing 
velocity), the* estimate R»-(n) will be in Region 2*. In this case, 

? SC 

R--(n) is adjusted by adding the PRF. 
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(c) If the coarse estimate is between 0 ami 2987 Hz, the fine 
estimate is not altered. 

It is assumed that the resolution of the ambiguity is carried 
out ideally, so that it does not contribute to the RMS tracking error. 
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APPENDIX C 

TIME MULTIPLEXED SINGLE CHANNEL ANGLE TRACKING OF PASSIVE TARGETS 

by • 

Charles L. Weber 


1.0 INTRODUCTION 

* 

Angle tracking is carried out by sequential processing of monopulse 
signals. The angle error signals are derived by sequentially processing 
the monopulse error signals. Angle tracking is potentially more critical 
than range or range rate tracking since angle monopuTse signals are 
subject to angular (glint) and amplitude (scintillation) variations due 
to the target motions. The easiest method to reduce these effects is by 
postdetection integration over the RF frequencies before the tracking 
error signals are generated. 

In the angle tracking loop analyzed in this appendix and described 
in [1], single channel processing is used for generation of the angle 
error signals. The single channel is time-divisioh-multiplexed between 
the ^ + ^"^AZ’ ^'*‘^EL’ ^ " ^EL signals. The alternative 

is a two- or three-channel monopulse system. The primary advantage of 
the three- channel monopulse is that it allows angle error measurements 
to be derived from a single pulse. To reduce the effects of angular 
glint, however, several pulses are averaged to obtain smoothed data. 

The primary disadvantage of the three- channel monopulse system is the 
requirement that the two additional channels be phase and gain balanced 
to the first channel. Alternatively, the advantage of the single channel 
monopulse system is the simplicity in implementation and the noncriticality 
of phase and gain, at the expense of a 3 to 6 dB reduction in performance. 

A disadvantage is that it may be more suspect to target scintillation. 

As mentioned previously, this can be circumvented by frequency diversity. 

In the next section, the single channel monopulse system under 
consideration will be described, followed in subsequent sections by a 
performance analysis. 

2.0 SINGLE CHANNEL AMPLITUDE MONOPULSE 

A single channel monopulse system is described, vjhere the channel 
is time shared by the various error signals. The description of the 
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various radar tracking functions has already been presented in the Appen- 
dix B. From that description a block diagram of the single channel mono- 
pulse is shown in Figure 1. Familiarity of the previous appendix is 
assumed. Only those components of the radar that are pertinent to the 
angle tracking function are shown in Figure 1. 

2.1 Ku-Band Amplitude Monopulse Antenna 

Referring to Figure 1, the Ku-band monopulse antenna outputs a 
sum channel and two angle error signals, which are designated ^^d 
for the azimuth and elevation signals, respectively. With respect 
to the amplitude monopulse feed, sum and difference patterns of a mono- 
pulse antenna network are usually determined experimentally in order to 
account for the effects of coupling between horns and of introduction 
of high-order transmission modes at the mouth of the horn assembly. 

Studies of the actual parabolic antertna patterns [2-4] have 
shown that the normalized sum (voltage) pattern can be represented very 
closely by 

Eg = cosMi. 18 A^/eg) . (la) 

where a^ is the angle off the boresight axis, azimuth or elevation, and 
0g is the 3 dB beiamwidth. It is assumed that A^ and 6g are given in the 
same units. This representation of the normalized sum pattern is excel- 
lent to well below 3 dB from the boresight axis of the main lobe and is 
therefore a convenient representation when evaluating antenna tracking 
performance. At the lower gain portions of the main lobe and over the 
si delobe, the above representation is not useful. 

The measured difference (voltage) pattern for the same type 
parabolic antenna is well represented by 

E^ = 0.707 sin (2.36 A^/9g) , (lb) 

where the same comments given above apply. Significant differences 
between these approximations and actual patterns appear only at the 
edges of the main beam, well beyond the 3 dB beamwidth. 

In the angle tracking monopulse system under consideration, as 
well as in most cases, only the variation within the 3 dB beamwidth is 
of interest, in which case, the above approximations are very satisfactory. 


"’V —|, ■ -j;- 



Figure 1, Block Diagram of Single Channel Honopulse Angle Tracking Loop 
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The actual error pattern is approximately linear to at least one-third 
of the beamwidth off the boresight. The normalized error pattern 
v/hicii approximately results from the AGC and logarithmic discriminant 
operation in the receiver system (Figure 1) is linear over more than 0.6 
of the beamwidth. The slope increases as the direction of the received 
signal approaches the half pov/er point of the sum channel beam. The 
antenna error slope is designated as and is approximately 1.67 for 
the representations given above. This can be seen by observing that 

E. A 

^ ^ tan (1.18-^) = 1.67A^/0g 

s 8 

for small A^. Within the linear region, the error voltage is 

E. 

Error Voltage = -g- = (2a) 

s 

or, equivalently, 

^e|deg " ®B|deg ^m^ Ej ’ 

The actual value of the antenna error slope, varies over 
different monopulse antenna designs. A conservative interval is 

1.0 < < 1.9 , (3) 

with typical values being 1.2 to 1.7. 

The connection between e+a and (E^, E^) is that s+A is propor- 
tional to Eg + E^ and s-A is proportional to E^-E^, where the propor- 
tionality constant is the same in both cases. 

Since the antenna gain patterns discussed above have been nor- 
malized to 1, there is no loss in generality in equating 



E+A 

— E + E. 
s A 

(4a) 

and 

E- a 

= • 

S A 

(4b) 

Hence, 

E + A 
E- A 

E + E 



_ *^S A 

Ee-E, ‘ 
S A 

(5) 


H 4" A 

This will be used in the discriminant function to relate — r to E./E„. 

E - A A s 
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3.0 Description of Single Channel Honopulse 

The amplitude monopulse feed of the Ku-band antenna system outputs 
three channels of data— the sum channel the azimuth error channel 
and the elevation error channel as shown in Figure 1. As described 
in more detail in the previous appendixes of this report, the elevation 
and azimuth error channels time-share the single signal processing channel. 
The modulation by ±1 either passes the error signal unchanged or inverts 
the RF phase. After recombining the error channel with the sum channel, 
the time-sharing is as shown in Figure 2. 



Figure 2. Timing for Single-Channel Time-Shared 
Monopulse Tracking 


During the track mode of operation, the frequency diversity con- 
sists of 5 RF's, each of which operates for 21.44 msec. During each RF 
frequency, the single channel processing is equally time-shared among 
the four signals— E + a^^, slots. 

Each of these is processed for 5.36 msec, during which time 16 pulses 
are transmitted and received. Throughout this description, the system 
parameters corresponding to the largest range Interval of interest are 
being employed. In particular [7], 
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PRI = T = pulse repetition interval = 335 ysec 
P 

PRF = 2987 

Pulse Width = t = 66.4 nsec 
Frequency Dwell Time = Xp = 21.44 msec/RF 

Tg is the dwell time of which is one of four time slots for each 

frequency dwell, and the range interval is 7.2 to 18.9 nmi. 

Referring again to Figure 1, the various stages of IF and A6C 
are simplified and consolidated into one bandpass filter (BPF) and AGC. 
The AGC is shown but not taken into account in the angle track perform- 
ance analysis, since for angle tracking purposes, the logarmithmic 
discriminant is already self-normalizing. The primary purpose of the 
AGC is to maintain satisfactory signal level in the receiver and avoid 
receiver saturation. The AGC is herein assumed to be performing its 
operation, so that all normalization can be accounted for in the 
discriminant. 

The final stage of IF converts the sum plus error narrowband 
signals to I and Q (complex) baseband v/aveforms. The video filter 
before the A/D is a 2-pole Butterv/orth with transfer function 

where f = 237 kHz. At these longer ranges, R>7.2 nmi, where the pulse 
width is 66.4 ysec, the signal can be assumed to be undistorted by H^(f). 

The sampling frequency of the A/D is 480 ksps, v/hich corresponds 
to 32 samples over each pulse width x. This is a sufficiently large 
number so that the pre-sum can be accurately approximated by analog 
integration. 

The notation developed in the previous appendixes is maintained. 
In particular, 

1 = 1,,.., 5 corresponds to RF frequency 

J=l,...,4 corresponds to time slot; 

J = 1 implies 2 + A^^ 

J = 2 implies 2 - 
J = 3 implies A^^^ 

J = 4 implies - A^j^ 

K=l,...,16 corresponds to the pulse number for a given 
I and 0. 



1 

f 


I 


7 


Nov/ consider the complex baseband v/aveforms before the pre-sum 
and restrict attention to the azimuth channel. The I and Q v/ave forms 
are 

SI(I,J,K) = {z± kl^) cos [(wj + u^j)t+ 0j] + n(t)} cos (wjt) 

and (7a) 

SQ(I,J,K) = cos C((.)j + Wj,)t+0j] + n{t)} sin (wjt) 

(7 b) 

for the I channel and Q channel, respectively. The peak received signal 
pov/er is given by 

V = • 

v/here it is assumed that z±h (which are the factors representing antenna 
processing) are normalized to 1. The average received power is 

P = d P 

avg r t pr ’ 


where is the duty factor of the radar pulse train. When J = l, the 
amplitude is proportional to Z + and when J = 2, the amplitude is 
proporational to amplitudes Aj are assumed to be Rayleigh 

random variables, since the target is assumed to be Sv/erling I. The 
random phase 0j is assumed to be uniformly distributed over (0,2 tt). The 
Rayleigh and uniform random variables Aj and 0j, respectively, are inde- 
pendent. The pulse P(t) has magnitude of unity over the transmitted 
pulse width t = 66.4 ysec. The doppler frequency shift is represented 


by 0)^. 

The receiver system noise n(t) after the bandpass filter has- 
the narrowband representation 


n(t) = >^[N^(t) cos ujt - Ng(t) sin Wjt] , (8) 


where N_ and are independent zero mean narrowband Gaussian processes 
c s 

with one-sided power spectral density of Nq watts/Hz, and one-sided noise 
bandv/idth approximately equal to f (f = 237 kHz), which is the 3 dB 
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bandwidth of The total noise power of n(t) is 2Ngf^, and the 

power in N(.(t) and Ng(t) is Ngf^, each. 

Before the A/D, the complex waveform is given by 

SI + jSQ = [{wj + (o^)t+ 6j] + n(t)j e 

= ® ^ ^ + N^(t) + jNg(t) . (9) 

3.1 Pre-Sum (Range Gate) 

After the A/D, each received pulse is pre-summed over the pulse 
width. The timing of this pre-sum is dependent on the setting of the 
range gate in the range tracking loop. The range tracking loop is 
presently assumed to be tracking ideally . Variations in tracking errors 
will be taken into account subsequently when design margins are 
considered. 

With ideal range tracking, the pre-sum output is 

rkTn+T 

PSI(I,J,K) + jPSQ(I,J,K) = ^ CSI(t)+jSQ(t)] dt. (10) 

•^kT„ 

P 

Errors in range gate tracking would be accounted for by errors 
in the limits of integration in (10), where again, the digital sum of 
the pre-sum is being approximated by analog integration. 


3.2 Doppler Filters (Range Rate Tracking) 

Adjacent doppler filters are being used for accurate range rate 
tracking, as described in Appendix B, The frequency midway between the 
correct two adjacent doppler filters in the DFT is assumed to be at the 
frequency w^ + Aw^, v/here Ato^ represents the error in doppler tracking. 
The effect of this DFT filtering can be approximated by multiplying 
SI + jSQ by exp [“j(wjj- Au^)t] before the pre-sum integration and extend- 
ing the equivalent integration of the pre-sum and DFT from f=66.4 psec 
to t^ = 5*36 msec. This corresponds to the ^Iwell time at one RF 

and also corresponds to the bandwidth of each of the doppler filters of 
approximately =187 Hz. 


3.3 


Statistical Parameters of Doppler Filter Outputs 

The output of the doppler filters can then be represented by 
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A 1 

*^UI 

h 


li" V 


-VI_ 


The equivalent noise term for Uj is given by 

and the equivalent noise term for Vj is given by 

2t^ 

Nyi = I * 

E 


( 11 ) 


(12a) 


(12b) 


For the signal part of Vj in (11), the integral actually is over 
but the result is exactly identical to the integral over (0,xg) 
since coherent processing is being carried out over the entire RF dwell 
time. This is not the case for the noise, since the noise at this point 
is wideband with respect to and therefore the noise contribution 
to Uj is independent from that of Vj. 

It Is also important to note that, since the processing time of 
the DFT is taken to be the entire elapsed time of the 16 pulses, namely, 
(16)Tp = 5.36 msec, then the corresponding signal pov/er is average 
signal power over that same period, as opposed to peak signal pov/er. 
Additional comments concerning the filter response of the DFT to a 
range gated coherent pulse train are given in Section 6.0. 

The noise terms are identically distributed and independent. 

At this point, the signal and noise have both been accounted 
for during both the S + slot for the X^h RF, which is designated 

as Uj, and during the time slot for the X^h RF, which is desig- 

nated as Vj. This is because the desired signal portions of Uj and Vj 
are correlated. The noise contributions remain uncorrelated. 

Carrying out the indicated integral, the representation in (11) 
can be expressed as 
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X A 

AZ j0. / jAu. xA /Sin (aw. Xp/2)> 

I - A,, ^ (-T-)( (A„,V2) , 


From the signal standpoint, the effect of nonideal range rate 
tracking can be seen to be an attenuation of signal strength by the 
factor 

sin (Aw.xp/2J1^ 

_ J 

in power, and a phase shift of A'jj^Tg/2. The phase shift is inconsequential 
since 6j is already uniform and independent of Aj. 

For the present, the doppler tracking loop is assumed to be work- 
ing ideally, so that Au^-O. 

In (13), 


^cUI ^ ^sUI 


(15a) 


^cVI ^ *^sVI 


(15b) 


All of these random variables are Gaussian with zero mean and 


variance 


= ECN^jj] = = E[n/,P = N„Bp (16) 


where Bp is the one-sided bandwidth of the DFT filters, namely, 

Bp=187 Hz. These random variables are all mutually independent. 

With this information, the statistical parameters of the doppler 
filter outputs can be determined for a given RF. In particular, the 
power in Uj is 




The power in Vj is 


+ 2N(jBp 


(17a) 


+ 2NpBp (17b) 





sin (AtOjj T|^/2y' 
I A t ^7^2 ) 


2 


and the correlation (unnormalized) is 



E[UjV*] 


P 


* 

VU 




E(Af ) 


sin (Ato^T^/2j” 


2 


(17c) 


where (*) means complex conjugate and d^ is the duty factor of the signal. 

In general, Pyy v/ould be complex. If there v/ere an unaccounted 
RF phase shift in the antenna, DMA or DEA, then the angle error channels 
v/ould not be exactly in phase with the sum channel when they are recom- 
bined into one channel (Figure 1). When this occurs, Pyy becomes a 
complex correlation. In this development, no RF phase difference between 
the sum and angle error channels is assumed. The effects of such errors 
on angle tracking performance is examined in [6]. 

If the doppler tracking loop is assumed to be working ideally, 
then Aw^=0. Recalling that the peak received signal power is given by 



Pp, = E(Af) 

(18) 

When, 

When AtOjj= 0, 




(19a) 


Pv - ’’pr 

(19b) 

and 

PUV ' "tt^+V(=-"Az'V- 

(19c) 


Define the normalized correlation coefficient between U and V as 


2 A 
p 



When Aojy = 0, 


2 

P 


is given by 


(20) 
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2 

P 




( 21 ) 


2 

From an average power point of view, (s±a) is normalized to 
unity, particularly for small angle errors, so that 


1 _ 1 
P “ 2NnBp " 2Nn 

1 + Q_£ 1 + Q__ 

avg r avg r E 


( 22 ) 


The signal-to-noise ratio (SNR) P^vg r Signal -energy- to- 

noise-spectral-density over the time slot t^, which is the dwell time of 
2 + Equivalently, this is the SNR in the doppler signal bandwidth. 

The significant observation is that the normalized correlation 
coefficient is dependent on SNR, as given by (22). Note that, as 

t 


SNR-»-“, then p ->• 1 

and, as 

SNR->-0, then p ->• 0 , 


(23) 


This will subsequently be related to the performance of the monopulse 
angle tracking loop. 

With these statistical parameters of Uj and Vj. now determined, 
the remainder of the angle tracking loop can be described and the per- 
formance analyzed. 


3.4 Angle Tracking Discriminant Generator 

Referring to Figure 1, the detector forms the magnitude squared; 

|Uj|^; I = l,..,,5 (24a) 

|Vj|^; 1 = 1,.. .,5 (24b) 

where I represents the Kh RF. The only correlation that exists is 
between Uj and Vj for each I, As noted in the previous section, the 
statistical parameters of the complex Gaussian random variable {UpVj} 
are independent of I. 


I I 


= {E[(ln x)‘]} 


13 



Figure 3. Standard Deviation of the Random Variable z as a 

Function of Gorrelation 
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The postdetection integration forms 

(25) 


(26) 

as shown in Figure 1. The computation techniques for generating z are 
considered in Appendix B. 

The random variable z is the maximum likelihood (H.L.) estimate [5] 
of th logarithm of the ratio of the power in Uj to the power in Vj, i.e.» 

z = (27) 

where {'') means maximum likelihood estimate. 

It has been shown in [5] that z is the H.L. estimate of In (P^/Py) 
whether or not there exists correlation between Uj and Vj. 

In order to see how z forms a measure of antenna angle error, 
note first that z is an unbiased estimate of In (Py/Py) [8], namely, 

Z In [Py/Py] + x, (28) 

where [5]: 

E(x) - 0 (29) 

2 2 

where Is the variance of the estimate of Py/Py. Equivalently 

stated, 

E(Z) = In (Py/Py) (30) 

The RMS error of z , namely a , has been derived and numerically 
presented in [5], and the performance curves are shown here in Figure 3. 

The RMS error is dependent only on the number of frequency dwells (N) 
and the normalized correlation coefficient (p) between Uj and Vj, as 
given in equations (20), (21), and (22). It can be seen that, as the 
normalized correlation coefficient increases, the RMS error decreases 
and, in the limit, as p->l, the RMS error-s-O. 
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1=1 
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XlVii 


1=1 


where, in this instance, N = S. 

Finally, the logarithm discriminant former generates 


z = = In 


"N 

I 

1=1 ^ 
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The improvement In performance due to increases in signal-to-noise 
ratio is related to RMS error through the correlation coefficient. As 
indicated in (22) and plotted in Figure 4, the normalized correlation 
coefficient increases as the signal-to-noise ratio increases. 

In Figure 5, the results of Figures 3 and 4 are combined. The 
RMS error of the maximum likelihood estimate is shown as a function of 
the = ^/(2NgBp) for various values of N. In this system, the 

pertinent value of N is 5, since there are 5 RF's in the frequency 
diversity. 

The second step is to obtain the appropriate understanding of 
In (Py/Py)* which is the quantity being estimated in the maximum likeli- 
hood sense by z. Upon substitution of (19) into (30), 


In (Py/Py) = In 




^)^Pavg r ^ *^F 


= In 


. 9 2NnBp 


z P 


avg r 


(1 . 4)2 + _i!loiL 

.2, 


z'-p 


avg r 


(31) 


Inspection of (31) indicates that the output of the discriminant 
generator is a function of the signal-to-noise ratio, as well as the 
normalized measure of the angle error A/z. 

As indicated in Section 2.0, when the monopulse angle error is 
small, so that the error voltage is in the linear region, then the 
desired error voltage is given by 


c 

Error Voltage “ "f" “ 

or, equivalently, the actual angular error is 



B m Ej 





(32b) 
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standard Deviation 
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For small values of (a/s), corresponding to small values of 
antenna angular error In (Py/Py) as described in (31) can be expressed 
as 


In (Py/Py) 


1 + 2(a/£) + (A/X)^ 

1 + (SMR)"^ 1 + (SNR)“^ 

1 2(A/S) , (A/S)^ 

1 + (SNR)“^ 1 + (SNR)*‘| 




1 


[1 + (snr)“Li 


= 4p(4) 


(33) 


where 

SNR S . (34) 

2NgBp 

S^ is already normalized to unity, and the approximation In {1+Y)“y 
for small y has been employed. 

The first moment of the output ot the log discriminant generator 
In (Py/Py) without the above approximation is shown in Figure 6. It can 
be observed that, for small values of A/s, the value of In (Py/Py) 
varies almost linearly in a/s, and that the slope is dependent on SNR. 
The more drastic variations in In (Py/Py) occur at larger values of A/S 
vjhere the reduction in loop gain is more significant. It is henceforth 
assumed that the approximation in (33) is satisfactory, and that the 
effect of SNR on In (Py/Py) is through the multiplication of (a/s) by p, 
thus corresponding to a reduction in loop gain. 


4.0 PERFORMANCE OF THE AMPLITUDE MONOPULSE ANGLE TRACKING LOOP 

The results of the previous sections can now be combined to 
provide the RMS angle error due to thermal noise and target scintilla- 
tion when the target is in the linear region of the antenna angle error 
channels. Only one of the angle errors is considered. 

A block diagram of the equivalent model of the antenna tracking 
loop is shown in Figure 7. The result is a system whose angle tracking 
performance is identical to that in Figure 1 when the antenna is in the 
linear region. In Figure 7, 





t 
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0j = relative angle of the target with respect to an arbitrary 
inertial reference 

8^ = relative angle of the antenna boresight in either the azimuth 
or elevation direction 

~ angular error in either the azimuth or elevation 
directions 

kj^/0g = scale factor necessary to convert angular error to antenna 
amplitude monopulse output voltage 
4p = scale factor in (33) necessary to give In (Py/Py) 

X = equivalent noise of the system, taking into account the 
effects of thermal noise and target scintillation as 
described in this report and in [5] 
z = In (Py/Py) + x is the output of the log discriminant generator, 
as given in (28) and diagrammed in Figure 1 
K = equivalent loop gain of the antenna tracking loop, which 
takes into account all of the various contributing gains 
F(s) = open loop antenna angle control system and inertial stabili- 
zation transfer function. A discussion of such transfer 
functions is given in [6]. 


Equivalent 
Noi se 



Figure 7. Block Diagram of the Equivalent Mode", of the Antenna 

Tracking Loop 


I 
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Since a transient analysis and a study of the dynamic character- 
istics of the antenna tracking loop are not a part of this development, 
the effects of the loop gain and open loop transfer functions are taken 
into account through the noise loop bandwidth (Hz). 

Referring to Figure 7, the closed loop transfer function from 
to is equal to 


H(j2uf) 



Keq F ( j2rf) 


1 + 





(35) 


and the corresponding noise bandwidth is defined as 



(36) 


The dependence of on p has been explicitly indicated to 
emphasize the dependence of noise bandwidth on the signal -to-noise 
ratio by way of the normalized correlation coefficient p. For inter- 
mediate and high gain tracking loops, the dependence of Bj^j on p is 
approximately given by 

Bj^(p) ~ pBj^j(o-l) (37) 

where p = l corresponds to the infinite SNR case. This is the case since 
loop gain varies linearly with SNR. In addition, a high values of loop 
gain, the noise bandv/idth varies approximately linearly in loop gain 
and therefore SNR. This is indeed the case for all loop transfer func- 
tions that are anticipated for F(s). 

The transfer function from the equivalent noise input x to the 
angle tracking error is given by 

so that the variance of the angle tracking error is given by 

\ "x CN-O) . (39) 
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where 

2 2 

a (N,p) = 0 is the variance of the maximum likelihood estimate 
of In (Py/Py)) which is also the variance per sample 
of the noise in the equivalent block diagram of the 
angle tracking loop. The dependence on the correla- 
tion coefficient p and the number of RF frequencies M 
is emphasized in (39) 

Tg' = time per sample of the random sequence x by the D/A 
converter and succeeding boxcar detector in Figure 1. 

In this system, 

Tg = Tp(16){4)(N), where (40) 

N=5 is the number of RF frequencies. The noise samples 
X are assumed to be statistically independent from 
frequency to frequency 

16 = number of pulses during v;hich s + & is observed 
(see Figure 2) 

4 = factor corresponding to the number of different 
angular error measurements per frequency dwell time, 

namely, t ^ ” ^AZ* ^ ^ ^EL 

T = pulse repetition interval = 335 usee. 

r 

Upon substitution of these values into (40), 

t' = 320 T_ = 107.2 msec. (41) 

5 P 

Upon substitution, the RMS angle tracking error is then given by 



(42) 


or, equivalently, 

( 43 ) 

where 6g = 3 dB antenna beamwidth in degrees. 


deg 


4k. 


1 + 


m 


^P \-l 
avg r l 

2NoBp) 
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5.0 PERFORMANCE COMPUTATIONS 


For the Ku-band Integrated Radar and Comiiiunl cation Equipment for 
the Space Shuttle Orbiter Vehicle proposal [1], the radar parameters 
needed to compute the RMS angle error 1n (43) are given as 
eB=l-6» 

= 1 (m, P. 3-47) 

= 1 Hz {[1], p. 3-47) 

Tg = 107.2 msec [eq. (41)] . 

For these values of the system parameters, the performance of 
the angle tracking loop Is shown in Figure 8. The RMS angle tracking 
error is plotted against the signal-to-noise ratio j,/(2NgBp), 
which is the SNR at the output of the doppler filter (see Figure 1). 

The specification value for the Ku-band radar system Is shown in Table 1 
as 

'^SPEC ~ 0*153 deg . (44) 


This is shown as the dashed line in Figure 8. It is seen that the 
specification value is met for 


p 

avg r 
2NqBf 


> 


0 dB . 


(45) 


Table 1. Specification Radar Parameter Allowable Measurement Errors 



Raindom 3a 

Bias 3a 

Angle Rate Error 
Including Target 

Effects 



LOS 

Pitch 


0.8 rarad 
0.458 deg 

+2“ 

LOS 

Azimuth 


0.8 mrad 
0.458 deg 

±2° 

Angle Rate Error Not 
Including Target Effects 



LOS 
Pi tch 


0.14 mrad/sec 
0.008 deg/sec 

±0.14 mrad/sec 

LOS 

Azimuth 


0.T4 mrad/sec 
0.008 deg/sec 

±0.14 mrad/sec 




^ 11 1 j ! f ; I 
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In order to obtain angle track performance as a function of 
range, the following radar equation is used to determine P 


Pay, r , _ 

2(4Tr)^L[kT5]R^ 


2NqBf 


avg r/t2NoBp) 
(46) 


where all of the parameters are defined in Appendix A of this report 

The values used in the computation correspond to 

those in Table 3, 

Budget #2, of Appendix A, namely, 


G = peak antenna gain 

= 38.5 dB 

O' = radar cross-section 


\ = wavelength 

= 0.0218 m 

d^ = duty factor 

= 0.2 

P = peak transmitted power 

r 

= 50 y/atts 

Tg = dwell time corresponding to one , 


computation of the DFT = Bp 

= 5.36 msec 

k = Boltzmann's constant 

= 1.38X 10"^^ w/Hz. 

Tg = system noise temperature 

= 1500°K 

R = range, in meters 


L = losses, itemized as follows for angle 

track; 

see Appendix- A for more details 


Transmit 

3.5 dB 

Scan alignment 

0 

Lateral scan 

0 

Threshold 

0.5 

Processor loss 

1.0 

Range gate straddling 

0 

Doppler mismatch 

0 

PDI 

1.25 

Pre-sum 

0.6 

Signal loss due to coupling 

2.0 

(see Section 7.0) 


Noise increase due to combining 

2.1 

of £ and A channels 


(see Section 7.0) 


Total Losses 

10.95 dB 


With these parameters, the RMS angle tracking error is shown in 
Figure 9 versus range in nautical miles. For all ranges below 18 nmi, 
the RMS angle error is below the specification value, v/here the effects 


'VI 
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of both thermal noise and target scintmation have been taken into 
account . The performance curves do not take into account any degrada- 
tion in the angle tracking loops due to errors in the range tracking and 
doppler tracking loops. With the SMR design margin for angle tracking 
of passive targets at approximately 10 dB, there is room to allow degra- 
dations for these errors. 

In summary, the angle tracking specification for passive targets 
is met with ample design margin (=10 dB) by the Ku-band system. 

The other angle tracking specification is the estimation of 
angular rate. Note that the angular rate specification does not require 
that target effects be taken into account. Therefore, a nonfluctuating 
target can be assumed, thereby greatly increasing the average signal-to- 
noise ratio. The performance of angular rate estimation is not taken 
into account in this report. 

6.0 COMMENTS CONCERNING THE FILTER RESPONSE OF THE 

DISCRETE FOURIER TRANSFORM 

In this section, certain additional comments are made concerning 
the response of a DFT to a coherent pulse train. The particular wave- 
form under discussion is shovm in Figure 10. It is a complex baseband 
signal which is assumed to be ideally range gated as shov/n. It corre- 
sponds to the signal I + jQ in Figure 1. The signal is sampled and pre- 
summed to form the complex numbers s-j,...,Sj^ where N = 16. These are 
the inputs to the DFT. 

The time of pre-sum is equal to the pulse width, which is equal 
to the width of the range gate. During the remainder of the time, the 
signal does not exist, although implicitly it does, since coherence is 
maintained from pulse to pulse. 

In the frequency domain, the picture is as shown in Figure 11. 

In the upper part of the figure is the spectrum of the coherent pulse 
train. The bandwidth to the first crossover is I/t. The "needles" of 
the spectrum are spaced by the PRF. The one-sided width of each "needle" 

is (NT_)" , where N is the number of pulses observed. 

P _1 

The number of "needles" within the bandwidth 1/t is , where 

d^=x/Tp is the duty factor. The shift of the entire spectrum due to 
a doppler frequency shift is also indicated. 
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Range 



Figure 10. Coherent Pulse Train with Pre-Sum 

In the lower part of Figure 11 is the location of each of the 
filters in the DFT filter bank. In this instance, the fact that digital 
filters have periodic frequency domain responses is an advantage and 
all "needles" in the signal spectrum will fall v/ithin one of the filters. 

The bandwidth of each filter is (NTp)"^. If the pulse width 
remains constant and Tp varies, then 

(a) The bandwidth of the signal, 1/t, remains unchanged. 

(b) The. bandwidth of each needle varies. 

(c) The bandwidth of each filter varies exactly as that of 

the needle in the signal spectrum. 

Since the filter bandwidth is matched to the bandwidth of each 
needle in the signal spectrum, the output signal-to-noise ratio at the 
end of the observation period is approximately optimized. This is what 
is done, since the samples are made at the end of the observation interval. 

Whether in a search or track mode, a coherent radar can extract 
the maximum energy from the signal by use of the optimal receiver shown 
in Figure 12. All of the necessary ideal assumptions have been made, 
namely, doppler is assumed known, the range to the target is assumed 
known, and ideal coherence is maintained. The digital signal processor, 
which is a slightly suboptimal implementation of the above,. is shown in 
Figure 13. The only difference is the reversal of order of the ■’'•.tegrator 
over the range gate, and the mixing by the doppler frequency. In most 
coherent radars, the resulting loss is negligible (in the hundredths of 
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Doppler Ideal Coherent Addition Detector 

Assumed Range (Predetection 

Known Gate Integration) 


Figure 12. Block Diagram of Optional Coherent Radar Receiver 



Doppler Coherent Addition 
Assumed (Predetection 

Known Integration) 


Figure 13. Digital Signal Processor Implementation of the Coherent 

Radar Receiver 
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a dB). In this radar, it is shown in Appendix A that the loss is approxi- 
mately 0.6 dB. The purpose for making the change is that the doppler 
mixing and predetection integration are then most conveniently carried 
out by the DFT. The digital implementation is thus essentially the 
same as the optimal implementation. The time domain representation and 
the frequency domain representation have the same appearance for both 
implementations. 

7.0 ANGLE TRACKING COUPLING LOSSES 

In the power budget of Section 5.0, tv/o coupling losses are noted. 
These are determined as follows. In Figure 14, the power divide and 
power combiner for the sum and angle error channels are shov/n. In the 
power divide, the signal and noise powers decrease by the coupling 
factor A = -2 dB. The overall power in the e + a or e-a combined signal 
is therefore reduced by 2 dB, since the signal power in the angle error 
channel is very small. This is the signal loss due to coupling. 


E Channel 


+ Noise 




E Channel 
Comm Only 


A=-2 dB 


A Channel 


+ Noise 




E + A or E - A Channel 


Figure 14. Block Diagram of Coupling Losses 


For the noise, the thermal noise spectral density level in the 
sum and angle error channels are the same. These two noises are also 
statistically independent and their spectral density levels can there- 
fore be added. Therefore, the overall noise level is 


"ol = *"o| ■'"ol = ■'-“"o- 

“ E+4 E A 


r 


which is an increase in noise spectral density of 2.1 dB. This is the 
noise increase due to the thermal noise in the difference channel being 
combined with that in the sum channel. 
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RADAR PERFORMANCE IN PRESENCE OF MAIN LOBE GROUND CLUHER 

by 

Sergei Udalov 

1.0 MOTIVATION FOR GROUND RETURN CONSIDERATIONS 

The possibility of the ground return interfering with the Shuttle 
rendezvous radar operation arises when the target is in such a position 
that the radar antenna points towards the earth at an angle close to 
vertical. Figure 1 depicts such geometry. 

In addition to the geometry shown, i.e., the target being in the 
line-of-sight (LOS) to ground, the conditions which results in a poten- 
tially troublesome ground clutter are: 

(1) The ranging waveform repeating at a period shorter than 
the round trip time to ground, thus precluding the rejection of the 
ground clutter in real time. 

(2) The radar operating at the same frequency for a time inter- 
val longer than the round trip time to ground, thus precluding clutter 
rejection in the frequency domain. 

(3) The clutter return being of the same order of magnitude as 
(or possibly higher than) the signal, thus precluding clutter rejection 
by magnitude discrimination. 

If the system interpulse period exceeds the radar time to ground, 
receiver gating can be used to prevent clutter from interfering with 
radar operation. Also, if the operating frequency of the radar is changed 
so that the ground return arrives when the receiver has been already 
switched to the next frequency, the clutter interference is eliminated. 
Consequently, if either one of these clutter suppression techniques is 
available to the radar system, the clutter magnitude is irrelevant, 
except for the receiver front end saturation consideration. 

However, if the radar system is designed to operate at a rela- 
tively high PRF* and if the signal processing is such that the radar is 
required to dwell at each of the operating frequencies for a time which 


iC 

Such is the case of the Ku-band rendezvous radar system, which 
is limited to peak transmitter power of 50 watts. 


Shuttle 

Orbiter 



■Igure 1. Rendezvous Geometry* With Ground Clutter Return in Main Lobe 


Not to scale 



3 


exceeds the round trip time to the ground, the clutter will enter the 
radar receiver at ambiguous ranges (fold-overs). Consequently, it 
behooves us to examine the factors which may result in condition (3). 

We must also examine the possibility of alleviating the effect of ground 
clutter on system performance by exploring the possibility of adapting 
the aforementioned clutter rejection techniques to the system under 
consideration. 

2.0 GENERAL CONSIDERATIONS 

Before proceeding with the detailed analysis of the ground return 
and its estimated impact on rendezvous system performance, it is proper 
to present a few general considerations. These will then detenaine the 
specifics of the subsequent analysis. 

Consider first the general relationship betvreen the orbit alti- 
tude and the target detection ranges. The range of orbital altitudes, h, 
at which the rendezvous may take place is from 100 nmi (185.3 km) to 
500 nmi (927 km), with the lowest altitude presenting a potentially most 
troublesome rendezvous situation [1], In comparison, the range, R, for 
passive target detection and acquisition may vary from about 0.5 nmi 
(0.927 km) to 12 nmi (22.2 km). Hov/ever, the range advantage in favor 
of the passive target may be outv^eighed by the excessive cross-section 
area of the ground return. This is particularly true for the passive 
targets at the extreme of the acquisition range (-12 nmi) and for antenna 
line-of-sight (LOS) close to the downward vertical. It must also be 

pointed out that, for a near vertical incidence, the clutter power falls 

-2 -4 

off as h compared to R for the passive target return. Therefore, 
for a rendezvous at any nominal altitude, the target- to-cl utter pov/er 
ratio depends primarily on the target range and is relatively insensitive 
to altitude changes resulting from the rendezvous manevuers. For example, 
lowering the orbit altitude from 110 nmi to 100 nmi, as a result of 
rendezvous "from above," will change the clutter return by 20 log (110/100) 
or 0.83 dB, assuming a vertical descent on the target. In comparison, 
the change in target signal strength resulting from reducing the range 
to the target from 12 nmi to 0.5 nmi is 40 log (12/0.5) or 55.2 dB! 

This simple example points to the fact that, at ranges an order 
of magnitude shorter than the maximum detection range (12 nmi specified 
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for 1 m target), the signal advantage in favor of the desired passive 
target outweighs the possible deleterious effects of ground clutter. On 
the other hand, the potential degradation of system performance must be 
estimated quantitatively, as is done by the use of the time domain analysis 
in Section 3.0. 

Figure 2 shov/s, qualitatively, the time domain relationship 
between the various power levels which influence the performance of the 
rendezvous radar in the presence of the clutter. As shown in this figure, 
the pertinent power levels are; 

= signal (target) return power 

Pq = ground clutter return power 

P|^ = system noise in detect (or track) bandv/idth. 

The ratios of importance are then: 

Pj/Pc = si gnal-to-cl utter pov/er ratio 

Pc/Pn = clutter-to-noise power ratio 

Pj/Pj, = signal-to-noise power ratio. 

It is intuitively obvious that, if the rendezvous geometry* is 
such that clutter return is completely submerged in the 

thermal noise and, thus, the system performance is determined, such as 
in a "clear" mode, fay the Pg/P,^ ratio. When Pq/P^" the clutter may, 
depending on its nature, begin to degrade the target detection perform- 
ance and possibly target tracking performance of the system, as well. 

For P(;/Pn^^> radar performance may be affected consider- 
ably, but it can be regained by lov/ering the transmitter pov/er. The 
reduction of the transmitter power "buries" the clutter in the thermal 
noise, but it also reduces the received signal pov/er. Therefore, the 
transmitter power reduction should be carried out only if the available 
Pg/Pn IS high enough so that, v/hen the received signal drops in strength 
(relative to P^) along with the clutter, the resulting Pg/Pp is suffi- 
cient to perform the appropriate radar function. 

Finally, if Pg/Pj.<U, where U is some minimum ratio required to 
perform a designated function, the radar becomes inoperative and "fixes" 
have to be employed. Some of the potentially usable fixes are described 
in Section 4.0. 


Rel ati ve 
Pov/er 

(not to scale) 



= target return povjer 
Pj, = ground clutter return pov/er 
P^ = system noise in detect (or track) bandwidth 


Figure 2 . The Degree of Clutter Interference in Time Domain Depends on the Relative 

Magnitudes of P5/P0 and Pq/Pj, Ratios 


The acceptable value of U depends on the mode in which the radar 
is operating. We will consider radar operation with clutter in terms 
of two requirements: (1) a general requirement stated in the original 

specification [1], and (2) a specific requirement [2] pertaining to the 
rendezvous for the Skylab reboost mission. 

2.1 The General Requirement 

The general requirement for radar operation in the presence of 
the ground clutter is defined by paragraph 3. 2. 1.3. 5.1 of [1]: 

"3.2.1 .3.5.1 Ground Clutter . The radar shall meet the perform- 
ance requirements of this specification in both the detection and 
tracking phases and in both the passive and active modes of opera- 
tion in the presence of transmitted power reflected from the surface 
of the earth. Detection and tracking shall not be degraded by 
earth- reflected signals (clutter) for any antenna pointing angle 
which results in a total signal-to-clutter ratio of 20 decibels or 
greater. The clutter may result, from main-beam and sidelobe 
reflections or sidelobe reflections only. The reference signal 
level shall be equivalent to that provided by the reflection from 
a 1 -square-meter target. Clutter signal levels shall be computed 
in accordance with the clutter model given in Appendix VI. ... " 

From the general requirement, it is evident that as long as the 
Pg/Pc ratio is kept at 20 dB or above, the system operation must be com- 
patible with the nominally specified operation "in the clear." Note 
that paragraph 3. 2. 1.3. 5.1 does not demand meeting system performance 
requirements for incidence angles close to vertical; yet there is no 
relief in either the range or the cross-section specification as a 
tradeoff for the capability of operating with the targets being vn'thin 
a near-vertical main lobe. The situation is different v/ith the Skylab 
reboost mission requirements. 

2.2 Skylab Reboost Mission Rendezvous Requirement (Preliminary) 

In comparison, the preliminary requirements [2] for the radar 
rendezvous for the Skylab reboost mission are as follows: 

Mandatory - Detect, acquire, and track the passive target at a 
minimum range of 0,5 nmi (0.927 km) with the ground 
clutter in the main lobe directed downward, including 
vertical direction. 
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Highly Desirable - Detect, acquire, and track the passive target at a 

range of up to 3 nmi (5.56 km) with the ground clutter 
in the main lobe directed downward, including vertical 
direction . 

The tradeoffs in this case are between the requirement to operate 

with the ground clutter in the downward- directed main lobe and the range 

capability of the system. Also, because the nominal cross-section of the 
2 2 

Skylab is 10 m (versus 1 m general specification), another tradeoff 
is offered for the requirement to operate in the presence of the main 
nadir clutter. 

The analysis which follows considers the radar performance in 
view of both the general system and the Skylab reboost mission require- 
ments. 


3.0 GROUND RETURN CALCULATIONS 

3.1 Transient Response of Main Lobe Nadir Ground Return 

For a relatively flat earth surface, i.e., surface void of 
topological irregularities, the strongest main lobe ground clutter will 
be observed at an incidence angle close to vertical. The primary reason 
for this is that the ground back-scatter coefficient is largest at 
vertical incidence. Also, because the reflective properties of the 
incrementally-illuminated areas of the ground change as the radar signal 
propagates along the ground, the power of the ground clutter, as seen 
by the radar, also changes with time. 

The general nature of a time-varying ground return can be under- 
stood in considering the simplified case of vertical incidence illustrated 
in Figure 3. From this figure it is evident that, as the leading edge 
of an RF signal hits the ground, progressively larger areas are illumi- 
nated and thus, theoretically, the total return increases as the time 
progresses. For most cases of practical interest, this is not the case, 
however. The main factors contributing to the "diminishing return" are: 

(1) the average reflectivity of the incremental areas decreases with 
the incidence angle and (2) the antenna beam shape attenuates returns 
at 0j^:sL0g, where 6g is the nominal antenna beamwidth. The latter condition 
occurs when the time required for the incident RF wavefront to spread 


1 t: 
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across the ground area, defined by 0g» is short compared with the duration 
of the radar pulse. This condition is generally referred to as "beam- 
width limited" operation. 

To obtain a highly accurate representation of the transient 
behavior of the main lobe ground clutter, one must first obtain an 
"impulse response" of the ground return by examining the temporal pattern 
of the power reflected off the incremental areas aa{ 0). Then, to obtain 
the shape of the ground return for any arbitrary v/aveform, one must per- 
form the convolution of the transmitted waveform v/ith the impulse response 
of the ground. The complexity of obtaining the impulse response by the 
aforementioned method is evident by considering the expression for the 
power returned from an incremental area aA(9); 


AP 


rec 


P. G^(0) 

■ ' ' o ' 4 — 

(4Tr)^d^(0) ° 


( 1 ) 


where “ power received from an incremental surface area ,AA(e) 

P^ = transmitted power 

G(0) = antenna gain in the direction of the incremental 
area 

0 = angle from LOS to incremental area 

A = operating wavelength 

d(e) = distance between radar antenna and incremental area 
AA(e) 

p = surface roughness 
A' = wavelength-dependent losses 

Oq(o,p,a') = average ground reflectivity (referred to as baek- 
scatter coefficient 6 in [1]) 

aA( 0) = incremental area. 

As can be seen from (1), the evaluation of the total return from 
an extended surface is rather complex and therefore, wherever applicable, 
various simplifications should be used. One such simplification is that, 
for a relatively narrow beam (i.e., 0g<3°), Pq may be assumed to be 
nearly constant over the illuminated area. This is particularly true 
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for the back- scatter coefficient of the earth, as represented by the 
RI/SD model described in [1] and in Addendum B. 

Another simplifying assumption is that, within the antenna beam- 
width, the differences between the ranges to the individual reflectors, 
i.e., dr, are insignificant compared to the total height, h, above the 
ground. Specifically, for a near-vertical return, the range R(e) to 
the reflecting areas lying along the perimeter defined by the antenna 
beamwidth is 


R{0) 


h 

cos.Og/2 ’ 


( 2 ) 


where h is the altitude and is the 3 dB antenna beamwidth, 
loss ratio Ln is then 



The range 


{3) 


For 0g= 1.6°, 


^(e) = (cos 1.6°/2)^ = (0.9999)^ 

= 0.9996 . 


(4) 


Note that the value shown in (4) is independent of altitude h 
and, obviously, it causes only a negligible reduction in the power levels 
contributed by the ground areas lying along the periphery defined by 0g 
intersecting the ground surface. 

Based on these two assumptions, a simplified analytical model for 
the transient behavior of the near-vertical return can be derived, as 
explained in Addendum A. Equation (19) of that addendum is of particular 
importance because it represents the transient response of the ground 
return to a unit step. Thus, using (19) of Addendum A, we have 

Ug'(t) = A(6g.h)b(1 . (5) 

A(0g,h) = " effective ground area, a function of h and 0g 

b = weighting coefficient, function of beam shape approxi- 
mation used 

h = radar altitude 


n 


Gg = 3 dB antenna beamwidth In radians 

t« = system "time constant," function of h. On. and antenna 
^ beam shape approximation. 

Note that both b and t^ are functions of the beam shape approxi- 
mation used. For a Gaussian approximation, v/hich we use here for our 
analysis, we have 

b = 0.720 

and tp = 0.180 0g^(h/c) , • ’ (6) 

where Og is in radians. 

The time varying power level of the vertical ground return can 
be computed from a range equation into which all the pertinent parameters 
have been substituted. Thus, 


Pc(t) 


2 2 


Oq A{eg,h)b(l-e~^^S 
(4ir)^ h^ L.,. 


(7) 


v/here = peak transmitted pov/er 

G = maximum gain of radar antenna 
A = system wavelength 
= transmitter losses 
Og = ground reflectivity coefficient 

and h, A(9g,h), b, and tg are as defined previously. The maximum expected 
vertical ground return power can be found from (7) by letting t»tg. 

The time behavior of P^Ct) can also be found from (7) by examining P^Ct) 
for all values of t . 

Let us compute the maximum values of P^Ct) for vertical return 
and for typical orbit altitudes, such as 100 nmi, 200 nmi, and 500 nmi. 

We can also then compute the temporal behavior of to determine the 

effects of various pulse widths on system performance. Table 1 shows 
the numbers pertinent to such computations. 

The power budgets for computing the value of a*'® shown 

in Table 2. 
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Table 1. Effective Ground Area and System Time Constant as Functions 
of Orbit Altitude (Antenna beamwidtli eg = 1.6° or 0.028 
radians. Vertical return) 


Orbit Altitude, h 
nmi ( km) 

A(0g jh) 

^0^ 

ysec 

40 log (h) 
dB-m^ 

n,2 

m 

dB-m^ 

100 (185.3) 

2.10x10^ 

73.2 

0.087 

210.7 

200 (370.6) 

8.40x10^ 

79.2 

0.173 

222.8 

500 (927) 

5.25X10® 

87.2 

0.434 

238.7 


The values of Pgtt) shown are the asymptotic values, because 

they represent the value of Pg(t) for t»tQ. Figure 4 shows how these 

asymptotic values are approached as the function of time. The responses 

* 

shov/n are the result of transmitting either the classic unit step excita- 
tion pulse, or a pulse whose rise time is short compared to t^, but 
whose duration exceeds tg by at least an order of magnitude. 

The most important feature shown by the ground response curves 
of Figure 4 is that, for values of t>0.7 nsec, all of the responses 
are within about 1 dB of their asymptotic value. Considering the pulse 
width used in the radar system being built by the Ku-band equipment sub- 
contractor (Hughes Aircraft Company), this means that all pulses equal 
to or longer than 2.07 ysec v/ill result in the full magnitude of the 
clutter return for all orbit altitudes of interest. 

The importance of transient analysis of ground return is also 
indicated by the set of curves given in Figure 5. The curves shown 
there illustrate the effect of orbit altitude on the rise time of the 
main beam vertical ground return. 

It can be seen that the rise time of the return is similar to 
that of a single-pole RC low-pass circuit responding to either a step 
or a long pulse input. We also know that the time constant, tg, which 
determines the “speed" of the response, is a function of altitude and 
the antenna beamwidth. 

Again, as in Figure 4, one can see that for t>0.7 to l.Q ysec, 
the ground return power reaches a level within 5» of the maximum v/ith 
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Table 2. Power Budget for Computing Maximum Value pQ(t) 


Term 

Value 

dB (+) 

dB (-) 

Pt 

50 w 

47 dBm 

— 


(2) (38.5 dB) 

77 

— 


0. 0218 m 
(2.18 cm) 


33.2 

l-t 

3.5 dB 


3.5 


2000 

— 

30.0 

b 

0.720 

— 

1,4 

*^0 

0.80 


1.0 


Function of h 

lOlogAe(h) 

— 

h^ 

Variable 

— 

40 log (h) 



= 124-MOlog Ag(h) 

■69,1 - 401og (h) 



= 54.9 dBm + 10 log Ag(h) 

- 40 log (h) 


Substituting the values of Table 1 into the equation in Table 2, one 
obtains: 

Altitude 'C' 


100 nmi 
200 nmi 
500 nmi 


Pp(t) 


max 

-82.6 dBm 
-88.6 dBm 
-96.6 dBm 





Ground Clutter Power (dBm) 



Time (ysec) 



Figure 4. 


Main Beam Ground Return Clutter Power Level as a Function of 

Time (Nadir Return) 




Received Pov^er Relative to Maximum Pp(t) 
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Figure 5. Effect of the Orbit Altitude on the Rise Time of Main Lobe Nadir 

Ground Return 


1 1 
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the eKceptlon of the h = 500 nmi curve which has a relatively slow 
response. 

Another interesting fact revealed by the curves of Figures 4 
and 5 is that, for t less than about 0.1 vsec, the ground return does 
not reach its full amplitude. Thus, for the radar mode using a pulse 
0.122 usee wide, the estimated clutter level maxima v/ould be; 


Altitude 

P ' (0.122 ysec) 

dB Below Maximum 

100 nmi 

-83,9 dBm 

-1.3 

200 nmi 

-91.7 dBm 

-3.1 

500 nmi 

-102.8 dBm 

- 6.2 


The fact that the clutter power drops as the transmitted pulse 
is narrov/ed can be used for improving the target- to-cl utter power ratio, 
Ps/Pg- A close examination of data in Figures 4 and 5 Indicates, however, 
that for the system parameters under consideration to provide worthwhile 
improvement, say 6 dB or better for the h=100 nmi case, the pulse v/idth 
must be about 0.025 ysec. Such a pulse width may be usable only at the 
short ranges, say less than 3000 feet, at which the P 5 /PQ ratio is 
already sufficient due to range advantage. 

Thus, one concludes that the prime purpose of the time analysis 
of the nadir return is to indicate at which pulse vn‘dth the system oper- 
ates in the "beamwidth limited" mode. As shown in Figures 4, 5, and 6 ., 
this mode takes place at x^ 2.02 usee and thus clutter return for greater 
values of x is the same as for x~ 2.02 ysec. 


Effect of Target Range on Si gnal-to-Cl utter Ratio 


Having estimated the absolute values of the main lobe ground 
clutter for the worst case of nadir return, one can proceed to calculate 
the signal-to-clutter ratio Pj/Pq a function of range. The first 
necessary step for this is, of course, the calculation of absolute level 
of power return from a passive target of specified cross-section. Table 3 
shows the power budget for this calculation. 

To make the data available to us complete, v/e also need to know 
the absolute noise powers in the system bandwidths associated v/ith the 
various functions of the rendezvous radar. The noise levels can be 
obtained by multiplying the system noise density, v/hich includes system 



Main Lobe Clutter Power (dBm) 



Figure 6. Main Lobe Clutter Power Approaches Asjmiptotic Value as Pulse 
Width is Increased to About 2 vsec 
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noise temperature for the radar mode, with the appropriate system band- 
width. The latter can be assumed to be approximately equal to the inverse 
of the pulse width utilized during any particular mode. Table 4 shows 
the noise pov/br levels in dBm as a function of the transmitted pulse 
width [3]. 

Note that, for all pulse v/idths longer than 0.122 nsec, the noise 
pov^er is referred to the output of the pre-sum circuit, which is essen- 
tially a time domain matched filter for pulses longer »han 2.07 ysec. 

Figure 7 shov/s the power received from a passive target relative 
to range, ground clutter and system noise. The data presented in chis 
figure is based on solutions provided in Tables 3, 2 and 4, respectively. 
From Figure 7, it is evident that the specification requirement [1] for 

p 

detecting a 1 m passive target at 12 nmi ( 22.2 km) cannot be met regard- 
less of the orbit altitude if the target is along the downv/ard vertical 

direction from the Shuttle* On the other hand, Figure 7 shows that the 

2 

requirement for the Skylab reboost mission, v/hich calls for 10 m target 
detection at 0.5 nmi (0.927 km), can be met even at the worst case alti- 
tude of 100 nmi. The margin with which this requirement can be met will 
depend, of course, on the detection method used, but from data in Figure 7 
it is evident that P 5 /PQ of more than 20 dB is available at 0.5 nmi range 
for a of 10 m^. 

From Figure 7, it is also evident that the main lobe nadir clutter 
exceeds the level of system noise, even for the narrow pulse operation 
with T = 0.122 usee. A 12 dB reduction in transmitter power v/ill bring 
the 100 nmi altitude clutter approximately to the noise level of the 
0.122 usee pulse operation. Similarly, a 24 dB power-down will bring 
the 100 nmi clutter level to about the 2.07 usec pulse operation level. 
But, surprisingly enough, the operation vrith the 0.122 usec pulse will 
then be thermal noise limited and not clutter limited! Figure 8 shows 
this phenomenon. 

The explanation for data shown in Figure 3 is obtained from 
Figure 7. Specifically, at a range of about 0.9 km, where the 0.322 usec 
pulse mode takes over, Ps/P^ about 20.5 clB and is .about 35,5 dB 

This statement does not pertain to cases v/here fixes such as 
frequency agility can be used to combat clutter. 




1 

r 


1 


Table 3. Power Budget for Computing the Absolute Level of Passive 

Target Return 


Term 

Value 

dB (+) 

dB (-) 

Pt 

50 w 

47 dBm 

— 


(2)(38.5 dB) 

77 




0.0218 m 


33.2 

Lt 

3.5 dB 

— 

3.5 

(4ir}^ 

2000 

— 

30.0 


1 

(specification) 

10 log a 

— 


Range 

(variable) 


40 log R 
(in meters) 


Ps(R) 

= 124 +10 log a 
= 57.3 dBm - 10 log o 

-66.7 - 40 log R 
- 40 log R 


Table 4. Noise 

Power as a 

Function 

of Transmitted 

Pulse Width [3] 

Transmi tted 
Pulse Width 
T (jisec) 

Pre-sum 

Ratio 

(N) 


Noise Power (dBm) 

Video 

After 

Pre-sum 

After 
Ooppl er 
Filter 

0.122 

N/A 

-96 

N/A 

N/A 

2.07 

1 

-108 

-108 

-120 

4.15 

2 

-108 

-109 

-121 

8.3 

4 

-108 

-112 

-124 

16.6 

8 

-108 

-115 

-127 

33.2 

16 

-108 

-118 

-130 

66.4 

32 

-108 

-121 

-133 




(-123 

(-135 




in manual) 

in manual) 


power (dBm) 


Figure 7. Power Received From a Passive Target Relative to 
Range, Ground Clutter, and System Noise (All 
power refers to pre-sum output) 


/- R= 0.5 nmi (0.927 km) 
\r See Note 2 


a = 10 m 


NOTE 1: Specification requirement for 
passive target detection 

NOTE 2: Potential Skylab reboost mission 
mandatory requirement 

NOTE 3: Potential Skylab reboost goal 
requirement 


-t t 

a = 1 m 


Pg (Passive Target Return) 


Main Lobe Nadir Clutter 


^R=3 nmi (5.6 km) 
r See Note 3 


h = 1 00 nmi 


-82.6 


System Noise 


■T= 0.122 nsec 



h = 200 nmi 


h = 500 nmi 


T = 2.07 vsec \ 
r /*T = 4.15ysec 


rx= \ 
r 8.3 usec 


-88.7 


-96.6 


R=12 nmi (22.2 km) 
See Note 1 



Passive Target Range 




Si gnal-to-Cl utter {P5/PQ) and Signal-to-Ncise (P^/P^) Ratios (dB) 


\ \ \ 
\A 

a=lOm^ \ 

a = 1 \ ^ 


V^n“\ \ 

\ 

For 24 dB power-dovrti \ 
and T = 0.122 psec, \ 

Ps/Pc " Ps/Pn 


h= 100 nmi, Main Lobe Nadir Clutter 


For: 

1) Full transmit power 
Pc/Pp 2) 12 dB power-down 
^ ^ 3) 24 dB power-down 

Note: Ps/Pq^Ps^^o all t 
except 0.122 psec 


cf = 10 m 


50m 100m 


Passive Target Range 


Figure 8. Si gnal-to-Cl utter and Signal-to-Noise Ratios Versus Range 


? 

for a 0=1 m target. The clutter signal for h=100 nmi is about 13.4 dB 
above system noise for t= 0.122 usee. Therefore, reducing the transmitter 
power by 24 dB places the clutter 10.6 dB below system noise. But the 
24 dB power-down also lowers the signal-to-noise ratio by the same amount, 
thus reducing this ratio from 35.5 dB to only 11.5 dB, a ratio which is 
far below the si gnal-to-cl utter ratio of 20.5 dB available at this range, 
i.e., approximately 0.9 km. 

Although the results of a 24 dB power-down will not affect the 

tracking function, because of an increased SNR available owing to tracker 

integration, the effect of power reduction must be examined in detail, 

2 

particularly if one is considering acquiring a 10 m target at about 
0.9 km, a task required by the Skylab reboost mission. 

Another issue which must be addressed when considering the effect 
of clutter is that of providing a 20 dB si gnal-to-cl utter ratio as called 
out by the specification [1]. From Figure 8, it is evident that, for an 
orbit altitude of 100 nmi and a 1 m^ target, the 20 dB Pj/Pq ratio vn'll 
be at about 1 km for the case of main lobe nadir clutter. Consultation 
of the clutter versus incidence angle curves of Addendum C permit one 
to construct a set of Pg/PQ>.20 dB range contours as a function of 
pulse width and incidence angle. These contours are shown in Figure 9. 
From this figure, it is evident that the longer pulses provide better 
clutter immunity at angles of incidence other than vertical. This is 
evidenced by the fact that, for any particular angle of incidence other 
than vertical, the ranges included within the Pg/PQi20 dB increases 
with the pulse length. This is opposite to the aforementioned short 
pulse at vertical incidence. Note also that, for vertical incidence, 
the P3/PQ is the same for all pulse widths.* 

The leveling of the contours at the longer ranges is due to 
thermal noise limitation. In other words, as the incidence angle 
increases, the Pg/P^<20 dB contour extends further into the range 
until the condition is reached where p2/pQ^20 dB=Pg/Pj^. 

The ranges at which the thermal noise limitation takes over can 
be obtained from data in Tables 3 and 4, and also from Figure 7. For 
example, the Pg/Pj^=20 dB for a 2.07 wsec pulse operation is reached 
at about 4.3 km; for 8.3 psec, it is reached at about 5 l:m, etc. 

For t 2^2 vsec, the system operates in the beamvndth- limited mode 
.and, hence, the short pulse (t«1 psec argument) improvement is not 
observed. 


h = 1 00 nml 



Incidence Angle (deg) 


Figure 9. Contours of 20 dB Si gnal-to-Cl utter (Ps/P^) 
Ratio Versus Angle of Incidence and Range 




It must be emphasized, however, that thermal noise contour limits 
shown in Figure 9 are determined by the pre-sum output bandwidth, which 
is approximately equal to the inverse of the pulse v;idth used. But, as 
shown in Table 4, the doppler filter output noise is 12 dB lower than 
the pre-sum output noise. Consequently, from the standpoint of thermal 
noise limitation, the ranges to the noise contours can be doubled. But, 
because the effect of doppler filters on the clutter level is difficult 
to predict, the plots of Figure 9 are presented in terms of pre-sum 
output ratios only. Once the effect of doppler filtering on clutter 
is established, one can work out the Pg/Pj,i20 dB contours, referred to 
as doppler filter output rather than pre-sum output. Considering- the 
complexity of the doppler structure of the ground return as the function 
of the incidence angle, this task may require analytical modeling and 
numerical evaluation. 


4.0 FREQUENCY AGILITY EFFECTS 


As indicated in the previous section, the requirement of the 

2 

system to detect and acquire aim target at 12 nmi cannot be met at 
angles close to vertical and at the lowest orbit altitude of 100 nmi. 

The question to ask is then; What techniques can be used to alleviate 
the ground clutter problem? 

One of the methods which suggests itself is the use of frequency 
agility. This method is of particular interest because the frequency 
agility is already employed by the radar system to combat the effects 
of target return fluctuation. Let us thus consider the numbers involved. 

The round trip time, which is the radar time to the nadir ground 
return is 


t 




Q 

where h is the orbit altitude in meters and C’'3xl0 meters/second. 
Figure 10 shows the values of ‘tpyipCGR) for the range of orbit altitudes 
from 100 nmi to 500 nmi. Now the radar operation, as it is presently 
defined by the contractor, requires 16 independent pulses at each fre- 
quency to form a 16- frequencies Fourier transformation. An additional 
pulse may be required for initialization, thus making it a 17- pulse 
group at each of five diversity frequencies. For the acquisition PRF’s 


. (GR), Minimum Radar Time for Nadir Return (msec) 


6.2 msec 


5.69 msec 



17-Pulse Group / I ( 

at 2987 kps / j ^ 

/ ' i 
/ 1 

/ 1 

I" 

2.37 msec 

/ 1 1 

/ 1 

/ 1 
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j 

/ 

I 

460 nmi j 

1.24 msec y j 

- 1 

)-■ i 

1 
1 
1 
1 

1 1 
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1 
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Orbit Altitude (nmi) 

Figure 10. Minimum Radar Time to Nadir Return Versus Orbit Altitude 
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of 2987 pps and 7177 pps, the durations of the 17-pulse groups are thus 

1^^(2987) = = 5,69 msec (9a) 

7^^(7177) = = 2.37 msec. (9b) 

These delay times transform into 460 nmi and 192 nmi clutter-free alti- 
tudes, respectively, provided the operating frequencies are also changed 
at the same intervals. 

Consequently, if we v/ere to construct the maps of altitude versus 
angle-off-nadir areas where one can expect the cl utter- free operation, 
one vrauld correlate the data of Addendum C v/ith the clutter- free altitudes 
indicated by (9a) and (9b). Such maps are shovm in Figure 11 for the 
tv/o PRF's under consideration. The corresponding time-frequency patterns 
are also shown. The significance of the maps is that, within the shaded 
boundaries, the ground clutter is above thermal noise at the output of 
the pre-sum device and, consequently, it is a potential threat. It must 
be noted that, outside the angle barriers (vertical boundaries) the 
clutter is below the thermal noise, yet it is present at the input to 
the system. However, above the orbit altitude limits (horizontal boun- 
daries), the clutter is virtually eliminated (except for some very weak 
multiple time around) at the input to the receiver due to frequency hop- 
ping. To distinguish the nature of clutter suppression, the symbols 
consisting of two orthogonal arrows have been attached to the appropriate 
boundaries. 

« 

Suppose now that one were to consider a time-frequency pattern 
which would change the frequency at the interval of i.e. , 

every 1.24 msec. The area of clutter interference would then be as 
shown in Figure 12. From this figure, it is evident that such operation 
would clear all altitudes up to 5G0 nmi from primary clutter interference. 

At 500 nmi, the ground return will be falling with the operating fre- 
quency of the system, but it will be attenuated relative to clutter at 
100 nmi by 

or 14 dB. The 14 dB attenuation indicates that system range performance 
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can be at least doubled at the 500 nml altitude. It can be further 
improved by employing alternating frequency time patterns to avoid the 
return at 500 nmi. For example, a simple technique would be to intro- 
duce a fixed frequency offset to each alternate time-frequency pattern 
of five frequencies. 

Consider nov/ the effect of such "fast" hopping on the number of 
pulses available for integration at each frequency. Thus, 

N(7177) = (1.24x 10'^ sec)(7177) = 8.9 (11a) 

or about 9 pulses; 

N(2987) = (1 .24x 10"^ sec)(2987) = 3.7 (11b) 

or about 4 pulses . 

From (11a) and (11b), we conclude that "fast" frequency hopping will 
result in an insufficient amount of pulses to be coherently integrated 
at each frequency. In principle, however, one should consider the pos- 
sibility of storing the individual small subgroups of pulses until a 
sufficient number are available to perform a 16-point Fourier transform. 

If permitted by antenna dwell times and target return coherency, such 
a method would provide a very promising solution to the clutter return 
problem at all altitudes. 
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ADDENDUM A 

TRANSIENT RESPONSE OF THE NEAR-VERTICAL GROUND RETURN 

The analysis presented below takes into account the effect of the 
time dependence, as well as the radar pulse width dependence, of the 
ground return. It is also shown that the antenna beamwidth is a major 
factor which influences the target- to-cl utter power ratio. 

We begin by referring to Figure 1, v/hich shows the relationship 
between the quasi-stationary and the time-variable parameters which 
determine the transient behavior of the near- vertical ground return. 
Specifically, we assume that the altitude h varies slowly compared to 
other time-dependent parameters, and thus is quasi-stationary. We also 
assume that the average ground reflectivity Oq is virtually uniform 
within the limits defined by the antenna bearawidths under consideration* 
and that the illuminated ground patch is level. 

Thus, as the RF wavefront hits the ground at nadir point a, the 
boundary of the reflecting area begins to spread and, at some time t 
after nadir contact, it spreads to a circle containing point _b. Also, 
as shown in Figure 1, the boundary of the illuminated area is defined 
by the time-varying angle e(t). From the geometry shown, we have 

“ hi-Mt) 

where 

4v(t) = . (2) 

since we are considering only the monostatic radar return. 

Substituting (2) into (1), we have 

cos e(t) = . — L-r , (3) 

where v/e define 

T = = nadir return delay . (4) 


* 

See Addendum 
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Note that, for the case of h= 185.3 km (100 nmi), T is 


T • “ 1235.3 psec. 


( 6 ) 


Now, since we are considering cases where t«T and e«l radian, v/e can 
expand both sides of (3) into their respective series. Thus, (3) becomes 


1 _ e^(t) + 6^(t) _ 
' 21 4! 


I T ^ 2PT^ 


( 6 ) 


Retaining only the first lowest order terms of both sides, we have 

eV) = 2|. (7) 

Equation (7) establishes a relatively simple relationship between e(t) 
and t , i.e., the time after start of radar return. 

Now we will derive the expression for the time-varying area of 
the ground return. From Figure 1, it is evident that 

A(t) = irh^ tan^e(t) (8) 

= Trh^0^(t) for0(t)«l radian. 

Substituting (7) into (8), we obtain the expression for A(t) explicitly 
as the function of t : 


A(t) 


2tt h^ 
T 


(9) 


Equation (9) shows that the illuminated area varies directly as a function 
of t , provided that we do not consider the effect of the antenna direc- 
tivity pattern G(e). For a Gaussian beam shape approximation, the two- 
way normalized directivity (gain) pattern is given by 

G^(0) = exp (-5.552 0^/0g^) , (10) 

where 0g is the one-v/ay 3 dB beamwidth and e is the direction of signal 
arrival with respect to beam center. But, from (7), we know that e is 
a time-varying function 0(t). thus, we can rewrite (TO) as a time- 
varying weighing function: 

w[e(t)] = exp [-5.552 e^(t)/0g^] . 


( 11 ) 
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The weighing factor expressed by (11) can be applied only to the differ- 
ential areas defined by an expanding ring having a time "thickness" At, 

In other words, (11) can be applied only to the "impulse response" of 
the ground. To obtain the expression for such a response, we use (9). 

Thus, the differentia! increase in area AA due to time increase At is 
computed by subtracting A(t) from A(t+At); 

2 2 

A = A(t+At)-A(t) = -^(t+At-t) = -^At. (12) 

The magnitude of the impulse response is obtained by dividing AA by At 
and letting At-->-0. Consequently, 




lim 
At ■+ 0 


M = 

At 


lim 
At + 0 


A(t+At) - A(t) 
At 


2tt h^ 

T * 


03 ) 


Because this is a definition of the derivative, we thus obtain: 

«s(t) = t ^ ' (H) 

for the unweighted impulse response of the ground. Note that within the 
assumptions used, i.e., e(t) « 1 radian and 8g « 1 radian, this response 
is constant, which is a characteristic of a perfect integrator. This is 
also evident from (9), which tells us that the reflecting area increases 
linearly with time. Since such a ramp-like response is obtained from an 
integrator when a unit step function is applied to its input, we could 
have obtained (14) by direct differentiation of (9). 

We can now proceed to apply the antenna pattern v/eighting func- 
tion (11) to the unweighted impulse response (14); 

2 

6^(t) = ^ exp [-5.552 , (15) 

v/here the prime indicates the weighted response. 

Substituting (7) into (15): 

2 

Sa^t) = -^exp C-(5.552)(2) t/T] 

2 

= exp [-t/tQ(0g,h) , ( 16 ) 
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where 


to(h.6g) = 


8g^T _ e|(Hh) 
(2){5.552) (2)(5.552)c 


= 0.180 0n (h/c) * 


Note that (16) and (17) can be used to obtain the ground return due to 
any envelope shape of the transmitted radar signal. For this, the 
envelope of the transmitted radar signal ^-'.s to be convolved with (16). 

But since we are primarily interested in ground response due to 
rectangular pulses, we can integrate (16) to obtain the v/eighted ground 
return waveform due to a unit step. Once v/e know the response to the 
unit step, the response to a pulse of any length can be obtained by the 
conventional method of using positive and delayed negative step functions. 

The weighted step response of the ground is the integral of (16). 

Thus, 




Jo 

2trh2r. -t'/tr 


h ^ / ^ g"^/^0 

T ^ ^0 


+ tft) 


v/here T and t^ have been defined by (4) and (17), respectively. The 
physical significance of (18) can be understood better if we replace the 
tg factor outside the parentheses by its definition given in (17): 


Ug'tt) = 


(2)(5.552) 


1 ■ 

4 1.388 




0 - . 


The term irh 0g/4 is the area of the ground within the cone 
defined by the 3 dB beamwidth of the antenna. Furthermore, the factor 
1.388 is the weighting factor for this area due to the Gaussian beam 
shape approximation.* Therefore, we can write (19) in terms of its 


(4) (1.388). 


In (19), we deliberately factored the coefficient 5.552 as 
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coefficients as 


Ug(t) = A(0B,h)b(l-e'’^'^S (20) 

t,2o2 

V n ®n 

where A(0g,h) = — — ~ = ground area, function of h and 0g 

b = weighting coefficient, function of beam shape approx- 
imation used (b= 1/1.388 for a Gaussian beam) 

tp = system “time constant,” function of h and 6g. 

Note that (20) is rather general and thus can be used along with (16) 
for various combinations of parameters, provided the assumptions used 
for expansion (6) are valid. 

Figure 2 shows the vertical ground return response to a unit 
step and to a pulse of finite duration. The latter can be computed 
from the following: 


fp(t) = A(03,h)b(l -e"^''^°) ; 0<t<.x > (21) 

fg(t) = A(0g,h)b(l-e"'^^O)e"^^'*''^^^O; t>x (22) 


where t is the radar pulse duration in time units compatible v/ith units 
of tp. 

It is interesting to note that the transient response of the 
near-vertical ground clutter to a Gaussian beam antenna illumination 
has the same form as the response of a single-pole low-pass RC filter. 

To compute the maximum target- to-ground clutter ratio, Pg/pQ> 
for any pulse width, one estimates the highest value of A' (^^jh) reached 
at x/tp and substitutes it into the appropriate range equations. Table 1 
shows the expression for P5/Pj. as a function of x/tp and of all other 
pertinent parameters. Note that substituting (21) and (22) into the 
equation of Table 1, the transient behavior of P5/PQ can be obtained 
for all values of t/tg. 


Fraction of Ground Return Pi 
Returned Frotn Area A(eD»h 


(1 - e 0) = Unit Step Response 



Figure 2, Normalized Near-Vertical Ground Return Response to Unit Step and 

to Pulse of Finite Duration 



Table 1. Maximum Target**to-Ground Clutter Pov/er Ratio for a Boresight 
Target and Near- Vertical Ground Return 


(t) = 


2 2 

Pt*^ " °t 

k3„4, 


(4Ti)^h^L 




{ 411 ) L Pt^ ^ Ag{6g jh jt) b Oq R Ag(0g»hjT) bog 


(4) 


4o. h' 
t 


T/tf 


,4.2, 




R'*iih^ 0 j b(l-e ■ “)og nbR’e^n-e ' “) Oj, 


where 

p^ = peak transmitted power 

G = bore axis antenna gain 

X = wavelength 

R = target range (meters) 

h = altitude (meters) 

= target cross-section (1 m ) 

L = system losses 

Ag(0g,h,T) = ground return area 

b = antenna beam shape factor (b=0.72 for Gaussian 
approximation) 

0g = antenna beamwidth (in radians) 

T = pulse width 

tg = 0.180 0g (h/c) for Gaussian beam shape 
Og = ground reflectivity 
Pg = target power 
Pg = ground clutter power 
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ADDENDUM B 

GROUND BACK-SCATTER COEFFICIENT MODEL FOR EARTH'S SURFACE 

According to the RI/SD specification [1]» the values of the back- 
scatter coefficient of the earth's surface v/ill be those associated with 
a smooth sea, i.e., approximately Seastate 1. The value of the back- 
scatter coefficient (s numeric, referred to as 6 in the RI/SD specifi- 
cation), according to the specification, ranges from 0,8 for vertical 

^ O 

incidence (specular reflection) to 10“ (diffuse reflection) for incidence 
angles greater than 20 degrees. Furthermore, for clutter estimation pur- 
poses, the Oq is assumed to be a linear function of the angle off the 
vertical in accordance with the following: 

a °0 

0 ° 0.8 

20 ® 10 “^ 

40® 10"^ 

60® 10“^ 

or greater 

v/here Oq is the back- scatter coefficient and e is the angle off nadir. 
Figure 1 shows the Oq value as a function of e for the ranges of e 
specified in the table. 

Of particular interest is the behavior of Qq close to nadir and 
v/ithin the angles commensurate with the radar antenna beamwidth. Figure 2 
shows this behavior,. From this figure, it is evident that, for calcula- 
tions of the near-vertical ground return, the value of the ground reflec- 
tion coefficient may be considered to be constant at 0.8 numeric. Note 
that, across the 6 dB tv/o-way antenna pattern, the maximum deviation 
from the 0.8 value is only 0.155 dB. 
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Ground Back-Scatter Coefficient (o«) 



Figure 1. Rockwell Model for Grou 
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(a) Ground Reflectivity Versus e 

0 (deg) 



Figure 2. Ground Reflectivity (a) Changes Very Slowly Compared to 
Two-Way Antenna Pattern Attenuation (b) 
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ADDENDUM C 

MAIN LOBE GROUND CLUTTER VERSUS INCIDENCE ANGLE 
AND PULSE WIDTH 

This addendum presents data on the behavior of the main lobe 
clutter at incidence angles other than vertical and also as the function 
of pulse width t. The data summarized below is based on the technical 
presentation [3^ dealing with this subject. It is included in this 
report to provide completeness to the discussion of the clutter effects. 

The video signal processor model assumed for the clutter and 
thermal noise calculations prior to doppler processing is shown by the 
the block diagram of Figure 1. 



Transfer Functions 

Video filter: |H,(f)l^ = 1 

’ lt(f/fc)’ 

f^= 237 kHz fpr t= 2.07 to 66.4 ysec 
MHz for x = 0.122 psec 

, iSin^TrNT^f 
Pre-sum: lH 2 (f)l = N- 7 — ^ 

{^Tgf) 

N= number of pulses summed 
Tg= sampling interval 
l/Tg = 483 kHz 

Figure 1. Block Diagram and Transfer Function of the Video Signal 

Processor 
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The clutter data shown in this addendum is referred to the output 
of the pre-sum unit and so is the system thermal noise floor. The latter 
is determined by the bandv/idth approximately equal to the inverse of the 
pulse v/idth and by the system noise temperature in the radar mode. 

table 1 presents the relationship between the pulse vndth used 
and the corresponding thermal noise power at the output of the video 
processor. 

Table 1. Thermal Noise Power Versus Pulse Vlidth 


Transmitted 
Pulse Vlidth 
X (nsec) 

Pre-sum 

Ratio 

Noise 

Power (dBm) 

Video 

After Pre-sum 

0.122 

N/A 

-96 

N/A 

2.07 

T 

-108 ■ 

-108 

4.15 

2 

, -108 

-109 

8.3 

4 

-108 

-112 

16.6 

8 

-108 

-115 

33.2 

16 

-108 

-118 

66.4 

32 

-108 

-121 

(-123 in manual) 


Note that, for the short pulse operation, i.e. , t= 0.122 sec, 

the pre-summer is bypassed. Figures 2 through 8 show the variation of 

main lobe clutter povjer P^ as a function of the incidence angle 0 (angle 

off nadir) and pulse v/idth x. Figure 2 is for video filter output Pq(6) 

and P„. Figures 3 through 8 refer Pf.(0) and P„ to the pre-sum output, 
n u n 
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Figure 8. Ground Clutter Versus 0 for t=64.4 psec 
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RANGING WAVEFORMS FOR SHORT RANGE OPERATION 
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APPENDIX E 

RANGING WAVEFORMS FOR SHORT RANGE OPERATION 

by 

Waddah Alem 

1.0 INTRODUCTION 

Several candidates for short range {R<3000 ft) radar waveforms 
in the Shuttle Ku-band radar are introduced. The basic requirement in 
choosing any one of these v/aveforms is its high range resolution accuracy. 
The range rate accuracies can be achieved by the use of a bank of doppler 
filters or via the a-p estimator at these short ranges. The ranging wave- 
forms considered here are (1) linear FM signal, (2) short pulse signal, 

(3) tone ranging signal, and (4) pulse coded signal. The resolution . 
capabilities, as v/ell as the ambiguity functions, of these waveforms 
are presented and compared. 


2.0 AMBIGUITY FUNCTIONS AND ERROR VARIANCES 

The -ambiguity function of a radar signal is defined in [1-4] as 

x(t,<|)) = [ u(t) u*{t+T) exp [-j 2it (|) t] dt , (1) 

■•—CO 


where u(t)‘is the complex envelope representation of the transmitted signal 

s(t) exp [j0{t)] (Ref. [1]), T is the time delay between actual reception 

and anticipated reception, (j> is the doppler frequency displacement caused 

by the target radial motion, and (*) denotes the complex conjugate. The 

properties of the ambiguity functions are v/ell documented in [1,4], The 

2 2 

minimum error variance of the target range 0 ^ and the target velocity 
are given in the literature [1,4] as 


4(2E/Nq) 


4(2E/NQ)a2 


( 2 ) 


( 3 ) 


where p is the rms bandwidth and a is the rms time duration*, that is. 
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( 1 / 
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2 



^ 2 E J 

ru(f)u*(f) df 

-CO 

(4) 

and 

a = ^2^^" 

“ 2E J 

fW 9 

c u(t) u*(t) dt . 

—CO 

(5) 


E denotes the received signal energy, «q is the carrier frequency, U{f) 
is the Fourier transform of u(t), c is the velocity of light in free space 
and \ is the wavelength. It is shown in [1] that the system that leads 
to the realization of these error variances is a maximum likelihood receiver 
composed of a bank of matched filters so that only range ambiguity is the 
main concern in designing the radar signal. 

As seen from (2) and (3), the values of a and e provide an indi- 
cation of the measurement accuracies in range and range rate, respectively. 


3.0 LINEAR FM SIGNALS 


Linear FM is a means of achieving pulse compression. The complex 
baseband representation can be written as 


u(t) 


= I'ect 2“ ’ 


( 6 ) 


v/here T is the pulse duration, k is a constant determining the slope of 
the FM. The absolute value of k is given by 

B = g HZ. 


where B is the swept frequency band. The ambiguity function of the linear 
FM signal introduced in ( 6 ) is found in [4] to be 


X[t,^] 


ejir<|.T(i 


hi) sin [uT(l--i^)U-t--I^T)] ^ 

T T(l- 27 T) 


t| ^T 


0 *, otherwise 


(7) 


Figure 1 illustrates the ambiguity function of a linear FM pulse 
for TB = 5. It should be noted that more elaborate schemes exist (e.g., 
VFM) with improved ambiguity function at the expense of an increase in 
complexity in both the transmitter and the receiver. 
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In an actual implementation, a signal 

-r 

s(t) = cos((rtQt+-~) ; |tl = (8) 

is transmitted and then passed after reception through a matched filter 
of the form 

h(t) = A cos 1^1 ‘'I' • 

If the doppler frequency shift is w^=2‘ir(|), then the output of the filter 
will be [5]; 

fib j + Rt “I 

sin — (T- |xl) 0 ). 

So(t,td^) = A cos ; KHT. (9) 

(ij^ + Rt 

Since at close range the effect of doppler using doppler filters is essen- 
tially eliminated, the maximum of Sg(x,0) peaRs at t = 0. 

Figure 2 compares the response of the output of a matched filter 
for a rectangular RF pulse and a matched filter for a linear FM rectangular 
pulse. The compression effect is readily obvious from the figure.* 

The only major modification in the Ku-band blocR diagram is the 
substitution of the existing pre-summer v/ith a weighted pre-summing. The 
sample weights are chosen so as to simulate the effect of linear FM matched 
filtering. 

As was indicated in Section 2.0, the rms bandwidth (b) and the 
rms time duration (a) are indications of range accuracy and doppler accu- 
racy, respectively. For the linear FM signal given in (8), it can be 
found that the rms bandwidth is given by [6]: 



Pulse compression means more energy in a certain time interval 
which results in more accurate measurement capabilities. 
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where S(t) and C(t) are Fresnel integrals defined as 


S(z) = f sin (f t^) dt 

Jo ^ 

C(z) ^ f cos t^) dt . 

Jo 

2 

For large time-bandwidth product kT»l, the value of e is equal to 



This results in a minimum mean square error of 
^ k(2E/Ny)^^^ ’ 


( 11 ) 


( 12 ) 


Similarly, for large kT, the minimum mean square velocity error is found 
to be 


= -L( k .1/2 

4 4ttTE/NqJ 


To compare the performance of a linear FH pulse with a rectangular 
pulse with the same duration T and bandwidth k/2iT, the relation of the 
rms time errors is computed: 

"Rectangular” pulse rms time-error „ u /kTJ/2 
Linear FM pulse rms time-error ^ ' 

For kT/2Tr=400, the theoretical time error ratio (improvement) is 2S.65. 
When considering a linear FM signal or any other pulse compression as a 
candidate for close range Shuttle Ku-band radar, the effect of range 
sidelobes has to be taken into consideration, as well as means of reducing 
them by various weighting functions [1]. Another effect that should be 
taken into consideration is the transmitter tube phase ripple (b^) on the 
sidelobes. For the TWT used in the Shuttle, 


= 1(^1 0 
3 W J 


where AV/V is the voltage ripple ratio and is the electrical length 
of the device. 


X 


4.0 


SHORT PULSE SIGNAL 


Since the resolution Improves with decreasing the pulse duration, 
a short pulse should be considered as a candidate at short ranges (R<2000). 
The main advantage of a short pulse is its simplicity when compared with 
any other waveform format. 

The ambiguity function can be derived from (7) by letting k=0. 
sin [^T(l - -i^)] 


xC'Cj'J') = 




ttT 


|t|<T 


0 *, otherwi se 


(H) 


The ambiguity function given in (14) is shov^n in Figure 3. 

The output of a matched filter to the RF pulse with zero doppler 
was shov/n in Figure 3. Direct application of (4) results in This 

is due to the infinite rise time of the pulse. In actual implementation, 
the rise time and fall time of the pulse cannot be Infinite due to the 
limited bandwidth; hence, if the bandwidth of the pulse is restricted to 
be B Hz, while time width is kept approximately equal to T, it is shown 
in [6] that 


_ _L ir B T - si n If B t 

^ ~ y2 Si (ir B t) + (cos ‘nBT-l)/TrBT 


(15) 


where 


Si(z) = p 
^0 


sin t 


dt . 


For large time bandwidth product (BT>>1), the theoretical rms 
range error can be given as 


(i 


4 'B E/Nq 


-) 


1/2 


(16) 


Applying (5) to a perfect squarevmve results in a range rate error of 






2hT(2E/Nq) 


1/2 * 


For short pulses, the range accuracy improves at the expense of the 
doppler resolution. However, since in the range rate, resolution is 


I 



9 


obtained using a bank of 16 doppler filters and, since the doppler vari- 
ations are small at close range, a short pulse should be considered as a 
candidate. 

The theoretical rms range error for a pulse with duration T=0.1 
nsec and a time bandwidth product BT=100 as a function of signal-to- 
noise ratio per pulse is shown in Table 1. 

Table 1. Theoretical RMS Range Error for T=0.1 nsec, TB = 100 


E/Nq (dB) 

-20 

-10 

0 

10 

20 

\ (m) 

30 

9.49 

3 

0.95 

0.30 


5.0 FIXED TONE RANGING SIGNAL 

The fixed tone ranging technique consists of transmitting a 
properly weighted series of N single tone signals so that the lowest 
tone, which can arbitrarily be called satisfies the ambiguity require- 
ment while the accuracy is maintained by the highest tone, 

The accuracy of fixed tone measurements is derived in [7]. If 
the average power of the nth tone signal is then the square of the 
rms bandwidth is found to be 



2 

n 




(17) 


where Pj is the total power in all the N tones. 

If the suboptimum receiver shov^m in Figure 4 is used, then assum- 
ing that the period of the lowest tone signal is equal to the maximum 
expected target delay and that the ratio between tv/o consecutive fre- 
quencies is a constant equal to M (M- > 1), it is shown [7] that 
the variance of the range error is given by 


-L 

4Tr 


^*^ax 

PtN^^O 


“i 


ll/2 


M' 


N-1 


(18) 
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where is the loop bandwidth of the highest frequency tone and is 
the maximum expected range. 

As is clear from equation (18), the variance of the range is 
increased v/ith the increase in the number of tones (N) which results in 
more complicated' receivers. The variance also decreases with the increase 
of the ratio between tv;o consecutive frequencies at the expense of more 
bandwi dth . 

The decision algorithm uses all the outputs of the N channels to 
estimate the range of the target. This is done by tracking the phase of 
the highest frequency so as to locate the peaks of this sinusoidal wave- 
form within a period of the lowest tone signal. The tracking of the phase 
of the second highest tone helps eliminate half of the ambiguity of the 
first measurement by reducing the number of possible ambiguous high fre- 
quency peaks by a factor of two. The process is repeated N-1 steps until 
the ambiguity is completely resolved. 

This method is very accurate at the expense of complex hardware 
implementation. An additional factor that has to be considered if this 
method is chosen for further study for close range tracking is the neces- 
sary acquisition time required by the phase lock loops (PLL) before an 
estimate can be derived. This might be a major factor in determining 
the eligibility of this ranging waveform. 

6.0 PULSE CODED SIGNAL 

Phase coding techniques can be used for pulse compression purposes 
which, in turn, is a means of improving measurement accuracies. Two can- 
didates are considered, namely, maximal length sequences (pseudonoise) and 
Barker sequences (perfect sequences). The main purpose of this section 
is to investigate the possibilities of using these codes and to point out 
the various problems that have to undergo further study before deciding 
if a given code can be used as a short ranging waveform. 

6.1 Pseudonoise (PN) Sequences 

The generation and properties of PN sequences have been discussed 
in detail in the literature [8,10]. Their application in ranging and 
radar measurements is found in [1,1,7,11]. The most important features 
about PN sequences from the ranging point of view is that they are easy 
to generate using linear shift registers with proper tap connections [10] 


1 


i 


I 






12 


and their correlation function resembles that of random walk signals 
for large sequence periods. The period (N) of a PN sequence is given 
by 

N = 2 " - 1 , (19) 

where n is the length of the linear shift register. 

The autocorrelation function of a PN sequence, averaged over a 
whole period of the sequence, is shown in Figure 5. The idea behind any 
PN tracking system is to transmit a PN sequence and then correlate the 
received target return v;ith a shifted version of the transmitted sequence 
until the peak in the correlation function is detected. The amount of 
local PN generator shift is a measure of the range. If the chip time 
is A and the amount of shift is x chips, then the range R is given by 

R = 

R 2 • 

Since the autocorrelation function of PN sequences is periodic, it is 
important to have an initial estimate of the range so that the correct 
peak is detected. 



Figure 5. PN Sequence Autocorrelation Function 

The ranging capabilities of PN sequences can best be understood 
by the following analogy. If a CW system is paired to an unmarked yard- 
stick used in measuring distances, then a PN system is paired to a marked 
yardstick whose marks are the code bits. 
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The squared magnitude of the ambiguity function of a PN sequence 
with a period N is given in [11] as 

= |l aKafe^.„exp (20) 

k 

= for m=Q (mod N), ij)=0 
= 1 for mj*Q (mod N), <|i=0 
= 0 for m = 0 (mod N) , ^ 0 
= N + 1 , otherwi se . 

The ambiguity function in (20) peaks at m=0 (mod N) and ^=0, which cor- 
. responds to accurate range measurement and no doppler. 

A possible i deal i sti c block diagram of the tracking system using 
a PN waveform is shown in Figure 6. In this block diagram, the received 
signal is mixed down to baseband, then correlated v/ith a delayed version 
of the same (local) PN code. The delay is adjusted by stepping the local 
code until a threshold is exceeded, which means that the two codes are 
brought to within a fraction of a chip time. The switch S is then closed 
and a tracking loop is used to bring the two codes to within the required 
range accuracy [12-14]. The range estimate is obtained by recording the 
amount of delay that the local code has to undergo before it is in com- 
plete synchronization with the received code. The thresholds are chosen 
so as to satisfy certain prespecified probabilities of error along with 
the examination period [15], The accuracy in range is related to the mean 
square error of the tracking system. Expressions for the. tracking errors 
for various systems are given in [12-14]. It is to be noted that, if the 
chosen examination period is not equal to the code period, the partial 
autocorrelation function of the code has to be used in the analysis [15]. 

Two closely related factors have to be considered when analyzing 
the system in Figure 6: 

(1) An initial range estimate (designated range) has to be 
known to set the initial time delay. 

(2) The acquisition time has to be calculated to see if it can 
be tolerated at close range measurements. This, in turn, is a function 
of the accuracy of the initial estimate, because as the latter is improved, 
the ansount of possible shifts between the two sequences becomes less and, 
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hence, faster acquisition can be expected. 

The block diagram in Figure 6 is idealistic . In any actual imple- 
mentation, the correlation process could be done at the early and late 
range gates to generate a pair of signals which are fed back after pro- 
cessing to step the local PN generator. Doppler filters have to be used 
to eliminate the effect of target range rate. Also, the pre-summing has 
to be modified by proper weights so as to match the PN signal. 

As in any other pulse compression system, sidelobe effect has to 
be analyzed and accounted for. A recursive least square estimate of 
channel response has been analyzed in [16] to eliminate the sidelobes 
effect. 


6.2 Barker Sequences 

The Barker sequences are those whose ambiguity function in the 
absence of any doppler is characterized by 

(N ; k = 0 

x{k,0) = \ (21) 

(±1 ,0 5 k 7^ 0 


There are only 9 Barker sequences known to exist. These sequences are 
found in the literature [1,4, 10,11] and are listed in Table 2 for 
convenience. 


JL 

2 

2 

3 

4 

4 

5 
7 

11 

13 


Table 2, Barker Sequences 
W X(k,0) ; k= 0 , 1 , . ..,(N-1 


+ + 

- ^ 

+ + ^ 

+ f -f - + 

+ + ^ ^ ^ 

+ + 4 * — “ -^ 4 * 


2 + 

2 - . 

3 0 - 

4 " 0 

4 + 0- 

5 0 + 0 + 

7 0-0-0- 

n 0 - 0 - Q - 0 - 0 - 

130 + 0 + 0 + 0 + 0 + 0 + 
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Barker sequences are called perfect sequences because of the 
definition of their ambiguity function in absence of doppler. The prop- 
erty that the magnitude of the si delobes is always less than or equal to 
unity does not apply when doppler is present [17,18]. 

Code 13, v^hich is applied in practical systems, has the auto- 
correlation shown in Figure 7 and the ambiguity function shown in 
Figure 8 [4]. 



Figure 7. Autocorrelation Function of Barker Code of Length 13 

The ambiguity function stresses the fact that the use of phase coding is 
useless in the presence of excess doppler. However, since in close range 
tracking, doppler is small and almost negligible, phase coding remains 
a modulation candidate. 

The actual implementation of a Barker code in phase coding a 
radar signal for tracking purposes is similar to that of PN codes where 
the received code is correlated with its shifted replica until the peak 
in the autocorrelation function is detected. The amount by which the 
local code is shifted is the measure of target range. 

7.0 CONCLUSION 

In order to choose the best ranging waveform for close range 
tracking, a rigorous analysis of all the candidates is required. This 
analysis should take all actual implementation losses into consideration, 
as well as the resulting complexity of the receiver, before any final 
selection is made. 
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Figure 8. x|{t»i{')| for the 13-Element Barker Code 
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MAXIMUM LIKELIHOOD ESTIMATION OF THE RATIO OF POWERS 
IN TWO NARROWBAND PROCESSES 

by 

Charles L. Weber 
ABSTRACT 

The maximum likelihood estimate of the ratio of the power in 
two narrowband processes is developed. Its RMS error is determined. 

I 

The estimate of this ratio is used in all radar tracking loops 
of the Space Shuttle Ku-band integrated rendezvous radar/communication 
system. 



2 


1.0 INTRODUCTION 

In this report, v/e consider the maximum likelihood (M.L.) estimate 
of the ratio of the power in two narrowband processes. These processes 
are assumed to be Gaussian and, in general, will be correlated. The in- 
phasp and quadrature-phase components of each process are sampled N times, 
where the time between samples is assumed to be large enough so that the 
samples are statistically independent. 

The estimate of this ratio and its performance play a fundamental 
role for all types of tracking in the radar described in [1], 

■2.0 THE MAXIMUM LIKELIHOOD ESTIMATE OF POWER RATIO 

Consider the narrowband Gaussian process u(t), which we describe 
as 

u(t) = u^(t) cos tOpt + Ug(tj sin UQt , (!) 

v/here u (t) and u (t) are baseband, zero mean, stationary, and Gaussian 

c s 

random processes. 

The in-phase and quadrature-phase samples of u(t) are represented 
as the real and imaginary components of a sequence of complex Gaussian 
random variables, namely, 

U^ = u^^. + , i = l,...,N . (2) 

2 

The expected value of [U^. j is denoted by 

Pg = E[|U.|^] = for alii, (3) 

which is twice the average power in u(t), u (t), and u (t). Equivalently, 

E[u^(t)] = E[u^^(t)] = E[Uj^(t)] = (4) 

The same notation can be established for the other narrowband Gaussian 
process, 

v(tj = v^{t) cos upt + Vg(t) sin Wgt , 
with average power Py/2. 

With this notation, we can restate the estimation problem as 
fol lows; 


Let ^ and V be N-dlmensional complex Gaussian vectors whose N 
complex components V. (i = l,...,N) have magnitudes and phases equal 
to the magnitudes and phases of a sequence of statistically independent 
samples of the narrowband processes and v(t), respectively. 

Given these observations, the maximum likelihood estimate of Py/Py 
and In (Py/Py) desired, as well as the RMS error of the estimates. 

In Addendum A, we show that the M.L. estimate of In (Py/Py) is 
given by 


where 


and 


In 


[pyp^i = 


In 


[pyp„] 



( 5 ) 




N 

I 

i=l 


( 6 ) 




(7) 


The estimator in (5) says that the H.L. estimate of In (Py/Py) is 
equal to the natural logarithm of the M.L. estimate of Py/Py. In Adden- 
dum A, this is shown to be the case for any strictly monotonically 
increasing function of the estimate. 


3.0 STATISTICAL CHARACTERISTICS OF THE OBSERVATIONS 


At time t., when the i;yi^ sample of u^(tj, u^Ct), v^,(t), and v^Ct) 
are observed, the covariance matrix of the resulting four-dimensional 
Gaussian vector will be of the form 



Py/Z 

a 

a 

pyz 

1 c 

i -d 
1 
1 

d 

c 

G 

-d 

-p- 

1 

b 

d 

c 

1 b 

t 

Pv/2 


( 8 ) 


In Addendum B, the statistical characteristics of two narrowband 
Gaussian processes are considered. In particular, it is shown that the 


4 


narrowband processes u{t) and v{t) that are observed in each of the track- 
ing modes of the radar are "spherically synunetric." The physical meaning 
of a spherically symmetric random process is also described in Addendum B. 
Tne result is that a=b=0 in (8). In addition, the conditions under 
which d = 0 are also enumeratjd, and it is noted that these conditions 
are met by the narrowband processes u(tj and v(t). 

The covariance matrix in (8) therefore simplifies to 


P„/2 

0 

! c 
1 

0 

0 

Py/2 

! 0 
I 

c 

c 

0 

m J- p 

1 Pv/2 

0 

0 

c 

1 0 
1 

Py/2 


(9) 


With this covariance matrix, all in-phase components are statisti- 
cally independent of all quadrature components. As a result, each complex 
sample of u(t) and v{t) can be interpreted as two consecutive real samples. 

Stated equivalently, the observations can be thought of as 2N inde- 
pendent samples of the real baseband Gaussian processes u (t) and v (t), 

Vi W 

Where these zero mean processes for each sample have covariance matrix 


M, 


Pu/2 c 


c Py/2 


, i = l,...,2N 


( 10 ) 


The maximum likelihood estimate of Py/Py is now given by 


/ \ A ? 


2N g 
i=l 

2N 5 

^ 4 

CK 


01 ) 


The performance of this estimator, as well as z = ln Py/Py is 
developed in the next section. ' 


1 




r 
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4.0 PERFORMANCE OF M.L. ESTIMATE OF POWER RATIO 

In order to develop most 
the following notation is used. 


In order to develop most directly the pdf of defined in (11), 


The normalized correlation between and v^^. is given by 


2c 


p = 


^/V 


p < 1 . 


(12) ^ 


The inverse of the covariance matrix Mg is denoted as 


‘ 


'^ll ''^12 


^12 Wg2 


(13) 


With this notation, the pdf of x^ in (11) Is given by ([2,p. 50] or [5]): 

a 


where 


nfy ) = 2 [det (Wg)]”^^ ^ [w^^ Xg^ + Wgg] 

B (f ,f)C(w^l Xg^ + Wgg)^ - (Zx^w^g)^]^'^'^^^/^ 


n = total number of independent samples = 2N 


. B(a,p) = Beta Function = 

r(a) = Gamma Function = (a-1)! if a is an integer 

det Wg = determinant of Wg . 

Substituting (12), (13) and (15) through (18) into (14), the pdf 
of can be written as 


(14) 

(15) 

(16) 

(17) 

(18) 






1 


L°u (I-P^) ay2(l-P^)J 


[(N-I)IT 

{2N-1)! 


1 


2 Xg P 


I^Olf(l“P^) ay^(l-p^) 


CTyPy (1 - P ), 


where 




2-|(2N+l)/2 

(19) 

( 20 ) 


> 


6 


The estimate of Py/Py is 


& 


i> 


^ ^ A 2 

V’’\i - h = \ 


(21) 


Carrying out this transformation, 


a cl 


d z = 2 x, dx. 
a a a 


The Jacobian of the transformation is 

1 


iji = 




and the pdf of z^, after some algebraic manipulation, is given by 

a 


p(/a) = 


(2N-1)! 

[4] 

" (1 - 


L'^vJ 

\ 0,, / 


z^ by setting 

a 


C(N - 1)1]^ jjz^ + [a^/a^] (1 - 2p^)) + 4 (ay/oy)^ p^O - 

(22) 

The form of this pdf can be substantially simplified if we normalize 

(23) 


A \ 

h - ^ ■ 


Then,’ 


p(x) = A 


(x^-1) x'^"'' 


[(x + 1)^ - 


(24) 


where 


A = (l-P^)'^(2N-l)l/[(N-l)l]^ 

The performance of the estimate of 


(25) 


In Py/Py = In Py/Py = 2 = IP Z^ ^ lO (Py/Py) + lO X (26) 


The pdf of p(x) in (24) agrees exactly with that developed in 
[l,p. C-9]. 
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is evaluated by determining the RHS error. 

The first moment is initially evaluated and is given by 

E = iri (Pu/Py) ' 

By making the change of variable x=e'^, 

Efonx)''! = aT (In x) . ^(x±iu'!:l'te . 


= 2A 


r 


cosh (y/2) 


[2 cosh y + 2 - 4p^3^’^2 


dy 


( 27 ) 


(28) 


When k=l in (28), the integrand is an odd function over even limits and 
therefore vanishes. Hence, 


E (tTyp;;) = In (pyp^,) 

and In P^j/Py is an unbiased estimate of In (Py/P^) 
The variance of error is 

Var ^ 0 ^ = E{[lnx]^}. 

Letting k = 2 in (28),* 


a/ = E[(lnx)^] = 4 a| 


cosh (y/2) 


= 2 A 


i 


“ dy 

° [2 cosh y + 2 - 4p^]'^'^2 

n + e“^1 e-^y 


(29) 


(30) 


^ [e“^^ + (2- 4p^) e"-^' + 1]'^’^ 2 


y dy . (31 ) 


This integral has been numerically computed. The results are 
illustrated in Figure 1. The computed values of the rms error* 
agree with the results shown in [T, p. B-4]. 


* 


Equation (31) agrees with the performance result in [l,p. C-10]. 




■[ 


.. f 

J 


I 





p = Normalized Correlation 


Figure 1. Standard Deviation of the Random Variable z as a 

Function of Correlation 
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ADDENDUM A 

MAXIMUM LIKELIHOOD ESTIMATE OF POWER RATIO 
MAXIMUM LIKELIHOOD ESTIMATE 

The maximum likelihood estimate of the ratio of power in two complex 
Gaussian vectors is developed. If U_ and V_ are random vectors and a is a 
parameter of the probability density function (pdf), p(U,Vja), then the 
M.L. estimate of a is defined’as that a which is 

max p(U>V;«) - P(MJ«) • (A- 1 ) 

a 

If the pdf is unimodal, then a is that a which is the solution to 

0 p (U,V;a) 

= 0 (A-2) 

3o 

or equivalently, that a which is the solution to 

3 In p(U,V;a) 

— = .0 . (A-3) 

3 a 

If there is more than one unknown parameter, whether or not it is 
to be estimated, the M.L. estimate is given by the following. Suppose 
a and 3 are unknown parameters and a is to be estimated. The M.L. esti- 
mate for a and 3 is formed by simultaneously solving [6-7] the following: 

3 In p(ii,V;ct,3) ' 

0 

'3a 

a In p(li,l;a,3) 

— =0 (A-4) 

3 3 

for a and 3 . The solutions are the M.L. estimates for a and 3 . If only 
one of the parameters is to be estimated, the estimate does not change. 

We also note that: The maximum likelihood estimate of any strictly 
monotonical ly increasing function of a parameter is equal to that same 
function of the M.L. estimate of the parameter . 


n 


For example, tf we have and we desire the M.L. estimate 


Y = In a , 


(A-5) 


then we form 


3 p(U,e^) 


(A-6) 


and solve for y- By the chain scale for differentiation, however. 


3 p(U,e'') _ 3 p(U,,a) 3 ^ 

3 Y ” 3 a 3 Y * 


(A-7) 


Since 


= e*''^ > 0 for all y s 


the only way (A-7) can equal zero is for 


a p(,U,a) 


= 0 . 


(A-8) 


(A-9) 


since e^ never vanishes. 

Therefore, when determining the M.L. estimate of In (Py/Py)s we 
need only consider the M.L. estimate of Py/Py 

THE MAXIMUM LIKELIHOOD ESTIMATE OF POWER RATIOS 

The M.L. estimate of power ratios is developed in steps. Consider 
first real uncor related Gaussian samples. The pdf p(U,V) is given by 


P(M) = -7^ 


2 exp 

J “v °v 


(A-10) 


The M.L. estimate of 


A Pui^ 


(A-11) 


is desired. The second parameter in p(]J,^) is chosen to be Oy, so that 
p(0,^) can be written as 


t 


1 


1 


r -iN 


P(U,V) = 


2'n‘ a, 


exp 


^|jl!J||^ + R||V||^] . 


The M.L. estimates are solutions to 

s In p(U,V) 


3 R 


3 In p(U,V) 
3 cr.. 


= .0 


= 0 , 


12 

(A^12) 


(A-13) 


Without difficulty, the M.L. estimate for R is shown to be given 


by 


Also, 


R = 


vr = 


lu 


i2 • 


(A-14) 


U.I 

ml 


(A-15) 


since the square root is a strictly monotonically increasing function. 
In the case of power ratios, it is noted that 


•{A-16) 


when the random vectors are uncorrelated. We shall see that this is also 
the case when u(t^.) and v(t|) are correlated. 

CORRELATED SAMPLES 

Correlated Gaussian samples are now considered. Since the extension 
from one sample of u and v to N samples is direct, we restrict attention 
to one sample of u and y. The joint pdf is given by 




•T ■ 


p(u,v) = 


1 


Zv < f .. a .. n/i~-*p 


exp 


u V 


-1 


20 -P^) 


2 u V p 

o a 
u V 


Where p is the normalized correlation coefficient between u and v , 
defined by 

/ EiMl . 

%°v 

With the definition of R given in (A-11), p(u,v) can be v/ritten 




-1 


20 -p") 


2u V p 


There are now three parameters, so that the M.L. estimates are given bj 
the solutions of 


a In p(u,v) 
a R 


= 0 


a In p(u,v) _ 


a p 


= 0 


8 In p(u,vl _ 

a (% ) 


Performing the indicated operations, we obtain 


R = 


1 


n 2, 2 

(1 - P ) ^2 


[u'^ “ u V p M“] , 


1 

2 


a-p' 


VP 

R /r 


and 


O-P 


u^ 2^ 

2/\T-^'^ 


u V 

VT 


T+p' 


Lp (1 - P )J 


(A-17) 

(A-18) 

as 

(A-19) 

(A-20) 
(A-21 ) 
(A-22) 

(A-23) 

(A-24) 

(A-25) 


from (A-20), (A-21), and (A-22), respectively. 
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From (A-23), (A-24), and (A-25), the M.L. estimate of R is 


or equivalently, 

2 2 

In R = In [uVv ] . (A-26) 

The key point is that when the random variables become correlated, 

2 2 

the M.L. estimate for does not change. 

The extension of these results to complex correlated Gaussian vectors 
is straightforward, although somewhat more cumbersome, with the result as 
shov/n in (5). With the simpler development that is shown here, the M.L. 
estimate in (5) is the expected estimate. 


£ 






9 ) 


» 
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ADDENDUM B 


STATISTICAL CONSIDERATIONS OF TWO NARROWBAND 
GAUSSIAN PROCESSES 

Preliminary statistical considerations of the tv/o narrov/band Gaussian 
processes are presented which are necessary before the performance of the 
maximum likelihood estimate can be determined. 

The processes used in this estimation are marginally and jointly . 
"spherically symmetric." Spherically symmetric means that an arbitrary 
phase shift of the carrier frequency v/ill not change the statistical char- 
acteristics of the in-phase and quadrature-phase components of either process 

For spherically symmetric complex Gaussian vector processes, we can 
invoke Grettenberg's Theorem, [3 - 4], an extension of which is quoted for 
this application. 

Grettenberg's Theorem . Let W(t) be a complex (column) vector Gaus- 
sian process with mean zero. Then a necessary and sufficient condition 
that 

E[W(t)/(s)] = 0 (B-1) 

for all s,t is that w{t) and e'^®w(t) be identically distributed for all 
real e (i.e., spherically s^mimetric). The T in (B-1) means "transpose." 

Mutually spherically symmetric for vector complex Gaussian processes 
implies that the information is in the phase differences, and an arbitrary 
but same phase shift in all components of the vector process produces no 
change in the joint statistics. 

As pointed out above, the processes u(t) and v(t) are such that 
they are jointly spherically symmetric. It has already been assumed that 
the successive samples of {u(t),v(t)} are sufficiently separated in time 
so as to be independent. The utility of Grettenberg's Theorem here is 
when s= t, for which 

E[u^(t)] = 0 (B-2) 

E[v^(t)] = 0 (B-3)' 

= 0 . 


E[u(t)v(t)] 


(B-4) 
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From {B-2)» 

ECu^(t)] = E(u^.) “ + 2 j'e(Uj,^. Ug.) = 0. (B-5) 

Hence, 

E(u^ 2,) = E(u/p (B-6) 

and 

E(u^,.Uj,.) = 0 . ■ (B-7) 

Therefore, and are statistically independent and identically 
distributed. 

A similar statement can be made about v(t) from (B-3). 

From (B-4), 

0 = Etu{t)v(t)] = E[“crci - "si''cP ■ 

Hence, 

where c and d are real . 

The cross-correlation of u(t) and v{t) is defined as 

^ E[u*(t)v(t)] 

= E [u^, v^,. + V 3 ,.] + J E Cu„- V 3 , - u^j . (B- 11 ) 

For spherically symmetric processes, 

PqV " 

and is real . 

T 

The covariance matrix of [u -,u v ., v •} can now be written 

W I O I >■ 1 M I 

as 


1 ) 




3) 




0 


e 


-T- 


\ 


T7 


M4 « 


V" 

0 

c 

d 

0 

Pu/2 

“d 

c 

c 

-d 

Pv/2 

0 

d 

c 

0 

f 

V2 


(B-13) 


Next consider u(t) and v(t) expressed 1n terms of envelope and 
phase variations: 


U(t) = R^(t) COS [uigt - 8^{t)] 

v(t) = R 2 (t) cos [Upt - 02 (t)l 


(B-14) 


Then, 


E [Uj, v^3 = E [R.J Rg cos 0^ cos 0g] = c 

E [lig v^l = E [R-j Rg sin 0-j sin 0g3 " c 

E [u^ Vg] * E [R^ Rg COS 0^ sin Og] = d 

ECdgV^] - ECR^Rg sin 0-j cos 0g] » -d 


(B-15) 


Let <j» .= 02" ®1 * where ^ represents the "phase difference information" 
betv/een u(t) and v{t). Then, 


E(u,Vc) - iE 


R| Rg (cos i> + cos (i) + 20^)^! = 

Rg (cos <{> - cos (<{> + 20^ )j 


= c 


(B-16) 


= d 


ECUjV^) = ^E 


Rg (sin , - sin + 20^)U = -d . (B-17) 


For spherically symmetric processes, 0-j is uniformly distributed 
over (“It, it) . When this is incorporated into (8-16) and (B-17), 
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ECu^v^] = -E[UgVg] = i E ER^ Rg cos <].] = c (B-18) 

E(u^,Vg) = E(u 3V^] = l.ECR^RgSin^] = d. (B-19) 

Finally, if the phase difference (f has a symmetric pdf about 0, 
then E (sin i|)) = 0, and d = 0. In the applications anticipated with this 
estimator, the statistics of the estimator are independent of phase. In 
addition, the assumption that <}> has a symmetric pdf is a very reasonable 
one for most applications. This assumption is therefore made and, hence, 
d = 0. 

The covariance matrix in (B-13) now simplifies to 


r pu /2 

0 

c 

0 

0 

py? 

0 

C 

c 

0 

Py/2 

0 

0 

c 

0 



(B-20) 


Inspection of the covariance matrix in (B-20) shows that the only 
dependence is between the real part of u(t) and the real part of v(t), and 
the imaginary part of u(t) and the imaginary part of v(t). These correla- 
tions are also identical. Since everything else is uncorrelated (and 
there'fore independent, since all processes are Gaussian), the observations 
could equally well be modeled as 2N samples of two real baseband processes 
u^(t) and v^(t). 

When the pdf of the estimator is determined, v/e assume that 2H 
samples of u_(t) and v_(t) are observed. The resulting pdf of the esti- 
mator is Identical to that when. N complex samples are assumed to be 
observed from u(t) and v{t) with the covariance matrix given by (B-20). 


